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Preface 



In this work we study basic properties of unstable particles and scalar hadronic resonances, respectively, 
within simple quantum mechanical and quantum field theoretical (effective) models. The term 'particle' 
is usually assigned to entities, described by physical theories, that are able to propagate over sufficiently 
large time scales (e.g. from a source to a detector) and hence could be identified in experiments - one 
especially should be able to measure some of their distinct properties like spin or charge. Nevertheless, it 
is well known that there exists a huge amount of unstable particles to which it seems difficult to allocate 
such definite values for their mass and decay width. In fact, for extremely short-lived members of that 
species, so called resonances, the theoretical description turns out to be highly complicated and requires 
some very interesting concepts of complex analysis. 

In the first chapter, we start with the basic ideas of quantum field theory. In particular, we intro- 
duce the Feynman propagator for unstable scalar resonances and motivate the idea that this kind of 
correlation function should possess complex poles which parameterize the mass and decay width of the 
considered particle. We also briefly discuss the problematic scalar sector in particle physics, emphasizing 
that hadronic loop contributions, given by strongly coupled hadronic intermediate states, dominate its 
dynamics. After that, the second chapter is dedicated to the method of analytic continuation of complex 
functions through branch cuts. As will be seen in the upcoming sections, this method is crucial in order 
to describe physics of scalar resonances because the relevant functions to be investigated (namely, the 
Feynman propagator of interacting quantum field theories) will also have branch cuts in the complex 
energy plane due to the already mentioned loop contributions. As is consensus among the physical 
community, the understanding of the physical behaviour of resonances requires a deeper insight of what 
is going on beyond the branch cut. This will lead us to the idea of a Riemann surface, a one-dimensional 
complex manifold on which the Feynman propagator is defined. 

We then apply these concepts to a simple non-relativistic Lee model in the third chapter and 
demonstrate the physical implications, i.e., the motion of the propagator poles and the behaviour of 
the spectral function. Besides that, we investigate the time evolution of a particle described by such a 
model. All this will serve as a detailed preparation in order to encounter the rich phenomena occuring 
on the Riemann surface in quantum field theory. In the last chapter, we finally concentrate on a simple 
quantum field theoretical model which describes the decay of a scalar state into two (pseudo) scalar 
ones. It is investigated how the motion of the propagator poles is influenced by loop contributions 
of the two (pseudo) scalar particles. We perform a numerical study for a hadronic system involving a 
scalar seed state (alias the <r-meson) that couples to pions. The unexpected emergence of a putative 
stable state below the two-pion threshold is investigated and it is clarified under which conditions such 
a stable state appears. 

*Note that this thesis was revised on the 21st May 2013.* 
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1 Introduction 



At the beginning of the 20th century most physicists around the world were faced with all 
the crucial changes that certainly unexpected ideas concerning the microscopic structure 
of nature (later on incorporated in the theory of quantum mechanics) resulted in. In 
retrospect, this time marks a turning point not only in the thousands- year-long history 
of science, but also in the mere way of how human beings look at the world surrounding 
them, leaving behind all coming generations with a mixture of amusement and curiosity 
about the universe. It is pure irony: while a huge number of our ancestors believed they 
were so close of obtaining a deep and conclusive understanding of the world, it seems 
nowadays that this kind of search for knowledge may never reach a final end. We have 
arranged with this. In fact, there are less people trying to reach for the answers to all 
things. Thinking about the universe has changed. And so we started a new venture at 
the beginning of the 21st centurey since it is up to us clarifying what our ancestors have 
left behind. 

1.1 Historical remarks 

Besides philosophical and fundamental challenges after finding the appropriate mathe- 
matical formalism, (non-relativistic) quantum mechanics faced a huge problem in estab- 
lishing a theory of nuclear forces. In 1935, YukawsQ used field theoretical methods to 
derive the nucleon-nucleon force as an interaction through one-pion exchange |1|. Al- 
though this description finally turned out to be not the correct way, it gave motivation 
for a lot of new approaches and theories in particle physics during the next decades. We 
will not try to review all those ideas, failures and milestones. This short introduction 
shall only sketch the basic concepts for working with one special kind of framework called 
effective (quantum) field theory (EFT) when dealing with resonances. 

Physicists are somehow 'cursed': they believe that everything should be made as 
simple as possible, but no simpler]^] Especially the basic first lines in books of particle 
physics are based on this principle, where we try to build up all matter from some very 
few and hopefully simple blocks of matter called elementary particles. Such efforts were 
first not successful after experiments discovered more and more heavy particles, named 
hadrons, in the early 60s whose existence was not covered by the theoretical models 
constructed before. It was realized soon after that most of the new particles are very 
short-lived states, so called resonances, during scattering processes and do not manifest 

1 In 1949, H. Yukawa was honoured with the Nobel Prize in Physics for his work on nuclear forces, in 

which he predicted the existence of mesons. 
2 This quote is often attributed to A. Einstein. 
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themselves as propagating entities. They were usually identified from the decay products 
they disintegrated to, hence it became dramatically clear that they could not be taken 
as elementary in the above sense. 

A theory to describe all the new hadrons and their interactions in a unified framework 
was developed from this time on, where the foundations were laid by Gell-ManrQ [2], 
Ne'eman and Zweig [3j. One of the main questions to be answered in this context was 
if a classification scheme for the new and already known particles could be found. Gell- 
Mann and Zweig proposed a solution using group theoretical methods, namely they 
treated all the different hadronic states as manifestations of multiplets within the SU (3) 
(flavour) group. As a consequence, the physical interpretation of quarks, elementary 
particles that are the building blocks of the observed hadrons, became the major idea to 
turn ones attention to. Quantum chromodynamics ( QCD) is the theory of the interaction 
between such quarks. Without going into details, the problem with QCD is the fact that 
this theory is non-perturbative in the low-energy regime, for example characteristic in 
nuclear physics. Despite all the (technical) methods physicists and mathematicians have 
developed in the past to find solutions of non-pertubative problems, it was up to now not 
possible to solve QCD, although enormous efforts were made in recent years. By applying 
the fundamental theory on a discretized space-time lattice, performing all calculations by 
using a huge amount of computational power and finally letting the lattice spacing turn 
to zero, lattice QCD has become a well-established non-perturbative approach of solving 
QCD and though lots of poblems have not been solved yet, e.g. concerning the numerical 
implementation, this field of research is under continuous growth. Other approaches have 
been found by using holographic models and gauge/gravity correspondence, for instance 
to extract meson masses with good accuracy [4j[5]. 

A substantial progress in hadron physics was nevertheless made when the concept of 
an effective field theory was applied to the low-energy regime. Weinberg has pointed 
out the general ideas in Ref. |6j, i.e., the key point is to identify the appropriate de- 
grees of freedom and to write down the most general Lagrangian consistent with the 
assumed symmetry principles. By this, it is not necessary anymore to solve the un- 
derlying fundamental theory due to the fact that within the new framework the set of 
degrees of freedom ('the basis') are not quarks (and gluons) but composite particles, 
namely hadrons. We then only need to deal with effective fields that are obtained in 
the particular energy region of interest. Good introductions to EFT can be found in 
Refs. |7f{9]. 

Whenever studying effective field theories one needs to understand the 'replaced' the- 
ory, namely QCD, and its symmetries. The effective Lagrangian will have the same 
symmetries and some of them may be broken - so clarifying the properties of the QCD- 
Lagrangian would give constraints on the effective Lagrangian. In this work we will not 
focus on specifiying an effective model for low-energy QCD. Instead we want to investi- 
gate the influence of loop contributions on the mass and decay width of a scalar particle 

3 In 1969, M. Gell-Mann was honoured with the Nobel Prize in Physics for his work on the classification 

of elementary particles and their interactions. 
4 In 1979, S. Weinberg, S. L. Glashow and A. Salam were honoured with the Nobel Prize in Physics for 

their work on the theory of the unified weak and electromagnetic interaction. 
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within a simple toy model containing mesonic scalar degrees of freedom and interactions. 
It is therefore not necessary to summarize the properties of the QCD-Lagrangian as usu- 



ally done, e.g. see Refs. 10,11 , or to go into details of chiral symmetry. We rather refer 



the reader to short summaries on the topic, provided for example in Ref. 12 



1.2 Scalar quantum field theory 

At least in modern literature it seems 'more elegant' to treat quantum field theory within 
the formalism of Feynman's path integrals. Yet we will choose the canonical formalism 
here for pedagogical reasons and follow the presentations of Refs. [13, 14 during the rest 



of this section, but use the conventions as denoted in |Appendix B within this work 



1.2.1 The scalar Feynman propagator 

The free Feynman propagator for a scalar field S in position space is defined as a par- 
ticular expectation value of two field operators: 

A| ee (x -y) = -i(Q\T{S(x)S(y)}\0), (1.1) 

where T is the time-ordering operator. In order to give an explicit form of this function 
we need to consider the action of the field operators on the vacuum state. To this end, let 
us first look at the case x° > y° by writing Ag ^ free (x — y). By introducing superscripts 
at the field S denoting those parts that contain a creation or annihilation operator, i.e., 

sl+) s = / ^ ' (L2) 

with p 11 on-shell, we can immediately drop most of the emerging terms, 

A { - )hcc (x-y) = -i(0|^5^+^5^+5(-)5(+)+^5^|0) 

^3 i 

since the creation and annihilation operators combine in such a way that their scalar 
product vanishes. In exactly the same manner we find for the other case x° < y°: 

A ( + )free (x -y) = -i [ j^ x ~ y) • (1.5) 

Both results are strictly valid only for their specific time order, but can be combined 
when using Heaviside step functions: 

Af%x-y) = e(x°-y )A { s )bcc (x-y) + e(y -x°)A^ irc \x-y) 

^° _ q fl\ p -W-{x-y) i Q/„0 _ ~0\ p lp-[x-y) 



[ [9(x° - y°) e-^-v) + 9(y° - x°) M 

J (2tt) 3 2E p I y y 1 yy 1 



(1.6) 



3 



This result can be further simplified. First we split the integrand in its time and spa- 
tial coordinates and then perform the analytic continuation into the complex p°-plane. 
Realizing that the remaining expression is nothing else than the summation of its two 
residues, we can rewrite it as a contour integral according to the residue theorem: 



1 

2vr7 



dp 



-ip (x°-y ) 



(P° ~ E p )(p° + E p ) 



1 

2E P 

1 

2£ D 



,-ip (x°-y ) | 
Jp (x°-y°)\ 



x° > y° 
x° < y° 



(1.7) 



The simple poles are located at p° = ±E P . The contour path of integration C has to 
be closed either in the upper half plane (x° > y°) or in the lower half plane (x° < y ) 
in order to make the contribution from the semicircle tend to zero. Together with the 
spatial part the scalar Feynman propagator can be written as 



a free 
^5 



(x-y) 



d 3 p 
(2^r)3 

d A p 
(2vr) 4 



drP 

r e -ip-{x-y) 



1 



C ^ 

ip-(x-y)_ 



Mi 



P 



M^ + ie 



(1. 



where in the last step we have evaluated the p°-integration by shifting the two poles as 
shown in Fig. |1.1| {p^ is no longer on-shell). This is the famous Feynman prescription in 
which the negative pole moves off the real p°-axis into the upper half plane, while the 
positive one moves into the lower half plane. The propagator in momentum space now 
simply becomes: 



A free / 
S 



It should be stressed that in the rest frame the bare mass Mq is the total energy of the 
particle corresponding to the field S. If we switch to the complex p 2 -plane, then there 
will be a single pole just below the real axis at p 2 = Mq — ie. Such poles describe stable 
states. 



Figure 1.1: Contour path of integration C (black, dashed) enclosing the two simple poles 
of the scalar propagator in the complex p°-plane. 



With the result just obtained we may show that the scalar Feynman propagator as 
ffven in Eq. (1.9) is the Green's function of the Klein-Gordon equation which fully 
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determines the evolution of the field S: 



(n x + M 2 )A^(x-y) = J ^ (V + M 2 )e-^-,) 



p 2 — Mq + it 



i4 



(2vr) 4 

= -^(x-y) . (1.10) 
1.2.2 Kallen-Lehmann spectral representation 

The free Feynman propagator for a scalar field 5, defined as the two-point correlation 



function in position space given by Eq. (1.1), plays an important role in quantum field 



theory, e.g. the LSZ reduction formula for the S-matrix mainly requires these kind of 
correlation functions. Additionally, the propagator can be interpreted as the amplitude 
for a particle to propagate from the space-time point y to x. Note that there is no 
problem to calculate the amplitude in a closed form when the particle is free, while this 
statement is not true anymore for interacting fields. 

One can express the full interacting propagator as a similar product of field operators, 
but where the scalar product is taken between the new ground state |0) of the interacting 
theory: 

A 5 (x -y) = -i(n\T{S(x)S(y)}\n) . (1.11) 
Let us write down the general completeness relation as 

l = J>)(a|, (1.12) 

a 

where we do not know the explicit form of the basis. Yet the basis will surely contain the 
single-particle states and all many-particle states, both types collected in |A P ), that shall 
be relativistically normalized momentum eigenstates. One can choose the set of three- 
momentum eigenstates to be the basis, since the three-momentum operator P commutes 
with the full Hamiltonian. As a direct consequence of postulating Lorentz invariance of 
the theory, any eigenstate at rest |Ao) of the Hamiltonian can be boosted in such a way 
that every possible value of the three-momentum is achieved. These new states are of 
course also eigenstates of the momentum operator and the Hamiltonian. By assuming 
the many-particle states to satisfy the energy dispersion relation with m\ as the total 
energy of a state at rest, this is also true for the boosted many-particle states in |A P ). 
The single-particle states have invariant mass p 2 = m 2 , the total energy of an unboosted 
state where the particle is at rest. The above completeness relation is then given by 

»=|fi}<«l+E/(03^IV<Apl, (1.13) 
summing over all unboosted momentum states. Inserting the completeness relation 



between the field operators in Eq. (1.11) and considering x° > y° for the moment, we 
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are left with 

A ( s\x-y) = -iJ2 J { fj r2E \ x) ^\S(x)\\ P )(X P \S(y)\n) . (1.14) 



Here, we have dropped the ground state term since it usually vanishes |l5l. The remain- 
ing matrix elements can be simplified using the transformation relations of scalar field 
operators under the translation group, 



S( x ) = e ty - x S(0)e- l ^ x , (1.15) 

and the translation invariance of the ground state. Additionally, remembering that 
p° = E p and realizing that S(x) is a Lorentz scalar and that the ground state obeys 
Lorentz invariance, the propagator for x° > y° reduces to 

A«(* -„) = -<£ J 2E~JX) ^ P ' M l^l^(0)|Ao)| 2 • (1.16) 
An explicit calculation makes this clear: 

(n\s(x)\x p ) = <n|5(o)|A p ) e-^\ p0=Ep 

= (niu-'usm-'uiXp) e-^\ p0=Ep 

= (fi|C/- 1 C/5(0)L>- 1 |Ao) e-^\ p0=Ep 

= (O|5(0)|A ) e-^\ p0=Ep . (1.17) 

Here, U is some unitary operator that implements a Lorentz boost from the three- 
momentum p to 0. If we now define the spectral density function p(s 2 ) by the expression 

p(s 2 ) = Tr^^s 2 ~ ml)\(Q\S(0)\X )\ 2 , (1.18) 

A 

the full interacting propagator for the field S can be rewritten without neglecting any 
terms: 

As(x-y) = e(x -y )A { -\x-y) + @(y -x°)A { +\x-y) 



_ 7 .0^ / d 'P 1 p -ip-(x-y) 



|(Q|5(0)|A )| 



2 



(2vr)3 2E P (X) 
J {2,-Kp p z — mj K + ie 



^e^lfd^ ■ (1.19) 

(2vr) 4 ttJo p 2 -s 2 + ie V ; 
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In the last two steps an integration over p° was introduced. This yields a simple expres- 
sion for the propagator in momentum space: 

a S (p 2 )=- r°d, 2 2 p{ i\. . (1.20) 

n Jo P - s + ie 

Obviously, the full propagator is a superposition of free single-particle propagators, each 
with different mass, weighted with the spectral density function p{s 2 ). This function 
is positive semi-definite and Lorentz invariant, it furthermore vanishes if the running 
argument s 2 does not lie within the future light cone. The integral representation from 
above is referred to as the Kallen-Lehmann spectral representation, independently found 



in the early 50s by Kallen 16 and Lehmann |17| . If we are able to identify the spectral 
function of a theory, we can directly calculate its propagator. The spectral function is 
normalized to one, a direct consequence of the equal-time commutation relations for the 
interacting field operators, 

[S(x ,x),S{y Q ,y)]\ Xo=yo =^ 3 )(x-y) , (1.21) 
[S(xo,x),S(y ,y)}\ Xo=yo = [S(x ,x.),S(y ,y)]\ xo=yQ = . (1.22) 

See Ref. |14| for a detailed derivation. It can be shown that 

I POO 

i<5 (3) (x - y) = itf( 3 )(x - y) - / ds 2 p(s 2 ) , (1.23) 

n Jo 

which requires the normalization condition 

- / ds 2 p(s 2 ) = 1 . (1.24) 

7T Jo 

For typical theories all accessible intermediate many-particle states will appear in the 
spectral function as a continuous spectrum above some threshold (we assume a mass 
gap that separates single-particle states from the many-particle states). It is commonly 
accepted to interpret 'bumps' within this continuum as unstable particles or resonances, 
depending on the dimension of their broadness. The state with the lowest total energy 
s 2 = m 2 is a single-particle state. If it is prepared in its rest frame the value of m is the 
physical mass of the particle, the mass that enters the energy dispersion relation. In this 
work we will go into the details on how resonances can be related with analytic properties 
of its propagator. A first hint is given when we try to perform the integration indicated 



in Eq. ( 1.20 ). If there are really single-particle states in the spectrum separated from the 
continuous distribution, then we expect their propagators to show up since the spectral 
function should give zero until the threshold value 4m 2 : 

^ Jo P - s + ie 

= 2 Z 2 ■ + d ^ 2 Pa ' tX 2 2) - » (1-25) 

pz _ m i _|_ ie J Am 2 p z — s A + ie 
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where a.t. means 'above threshold'. The spectral function takes the form: 

p(s 2 ) = Zit5(s 2 - m 2 ) + p a . t Xs 2 ) . (1.26) 

Here, the quantity Z stands for the wave function renormalization constant, the prob- 
ability for the operator S(0) to create a given state from the ground state. Hence, 
single-particle states are manifested as simple poles in the complex p 2 -plane on which 
the propagator Ag(p 2 ) is defined, while many-particle intermediate states make the 
propagator to possess a branch cut starting at the threshold value. Such a singularity 
provides a purely imaginary discontinuity along the positive real axis, dictated by the 
spectral function above threshold: 

1 f°° dg2 Pa.t.(s 2 ) ^- P 2 
Ji m 2 P 2 ~ s 2 + ie 



After we made use of the Sokhotski-Plemelj theorem, it is obvious that the additional 
second term only appears when Rep 2 > 4m 2 . We will clarify the meaning of such a 
behaviour in the next chapter from a mathematical point of view and apply the conse- 
quences to particle physics. 



^ J Arrfi—'p* 



dp 



-p[°° d„^ ( " + p2) 

Jirrfi—p 2 



P- 



iPa.t.(P 



1.27) 



1.3 Resonances 

1.3.1 Decay width and Breit-Wigner parameterization 

In the first section we have found an analytic expression for the full interacting prop- 



agator, see Eq. (1.25), by using the Kallen-Lehman spectral representation. On the 



other hand there exists a direct connection between the propagator and the spectral 



function. By looking at the discontinuous behaviour shown in Eq. (1.27) we realize that 
the principal value (of the propagator) is purely real. Thus, we may take the imaginary 
part of the propagator right above the real axis and obtain a very important relation: 



lmA s (p 2 + ie) = -— 2 1 -Pa.t.iP 2 ) , (1-28) 



=>- -ImA 5 (p 2 + ie) = Ztt5( P 2 - m 2 ) + p a .tXp 2 ) 

= P(P 2 ) , (1-29) 

where we have applied the Lorentzian representation of the delta distribution function. 
The infinitesimal number ie inside the argument of the propagator may look superfluous 
but in fact we will see later that it is crucial. So, if we are able to construct a propagator 
by specifying the interaction terms which lead to a continuous spectrum above some 
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threshold value, then we can also write down the spectral function. This will be used in 
chapter 4 of this work to study resonances in a simple hadronic toy model. 

We now have to make some remarks on resonances. In the framework of quantum 
field theory the term 'particle' is assigned to excitations of the fields (like the scalar 
field S) that are able to propagate over sufficiently large time scales (e.g. from a source 
to a detector) and hence could be identified in experiments. Those entities would have 
distinct measurable properties like spin or mass, and consequently should satisfy the 
energy dispersion relation. It then depends on the definition how to fix terminology: 
a stable particle should be a particle that propagates over an infinite amount of time 
without decay, except when interactions with other particles occur. This leads us to 
unstable particles - it is not difficult to fill dozens of pages about this topic. In order to 
satisfy expectations for an introductory chapter let us proceed in huge steps and simply 
write down results of less illuminating derivations by referring to the literature. 

For example, we know that charged pions as part of the particle shower caused by 
cosmic radiation have a mean life time of about 10 -8 s and can be described nearly as 
stable as long as the propagation and interaction time is much smaller than the mean 
life time (including relativistic time dilation effects). Nevertheless, when considering 
time scales of some seconds those particles decay into other particles, namely muons 
and neutrinos or photons. Unstable particles also build up chemical elements (and their 
compounds) that emit radiation in form of decay products of their nuclei - in fact we 
are surrounded by usually invisible effects of decaying unstable particles. 

During the past century physicists have found some 'particles' - and constructed 
theories with them - that cannot be even named as unstable, because they would be 
extremely short-lived unstable particles with mean life times on the order of 10 -22 s. It 
therefore makes more sense to treat them like excitations emerging when investigating 
nuclear matter or when performing high-energy collision experiments. Travelling at the 
speed of light, these particles could only overcome a distance of about 10~ 14 m before 
decaying. Such resonant states can nevertheless be interpreted as fluctuations of a field, 
too, and so we may use the words 'unstable particle' and 'resonance' equivalently. A 
very helpful formalism to describe an unstable particle is to make connection between 
its mean life time r and its spectral function (or its propagator). Treating the particle 
decay as a Poisson process, one usually defines 

r = r- 1 , (1.30) 

where V is called the decay width of the resonance associated with a specific set of final 
states, namely its decay products. This directly yields an exponential decay law for the 



survival probability pit) of the particle in its rest frame [15, 18 : 

pit) = e~ rt . (1.31) 

The particle decay width is a measurable object and theoretical determinations can 
surely be treated as predictions. We will calculate decay widths in the framework of 
non-relativistic quantum mechanics and effective quantum field theory; this is done in 
both cases with perturbation theory. However, it seems after all very problematic to 
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define a propagator for a resonance if the decay width is comparable to the assigned 
particle mass M. Such a statement is sometimes made clear by a rough idea: Suppose 
we investigate a scattering process between two particles, then a resonant state would be 
the manifestation of the very short-lived unstable state created and somehow exchanged 
during the scattering process. In the end, of course, the resonance disappears and may 
decay into the former two incoming particles. During the short time of its existence, 
according to the uncertainty principle the unstable particle could travel for a time interval 
about of 1/M. But if T ~ M the resonance decays already after 1/r and so cannot 
propagate over the full interaction distance - the main problem contains the vague 
definition of the unstable particle's mass. 

For the moment, let us assume T <C M, so that there is no problem by defining a 



propagator A$(p 2 ) in momentum space (where Z = in Eq. (1.25)). Since we deal with 
a particle that can decay, i.e., disintegrate into other particles, we are clearly faced with 
a problem of interacting quantum field theory. The complex pole of the free propagator 
situated at p 2 = Mq — ie can no longer be assigned to the now unstable particle, but it 
is also not obvious what exactly happens with this pole when the interaction is turned 
on. The only thing we can say for now is that there appears a branch cut above some 
threshold value and, respectively, the spectral function shows a continuous spectrum with 
maybe a narrow bump roughly at the mass M of the unstable particle. The mass M 
differs from the bare mass parameter Mq due to mesonic quantum fluctuations which can 
be included in a Dyson resummation scheme. The scalar propagator for the unstable 
particle will then be dressed by the self-energy contribution H(p 2 ) (see chapter 4 for 
details) that is specified by the interaction term(s) appearing in the Lagrangian of the 
theory: 

W) = p^rnpj ' < L32 > 

In general, the self-energy will turn out to be a complex-valued function with non-zero 
imaginary part, so we can drop the infinitesimal number ie. Whatever we may think to 
be 'the mass' of the unstable particle S, a reasonable definition is the real part of the 
denominator of the dressed propagator: 

M 2 = Ml + Ren(M 2 ) . (1.33) 

However, if the imaginary part of II(p 2 ) contributes to the denominator, the pole of the 
propagator is displaced from the real p 2 -axis, determining the mass and width of the 
resonance through 

r pU = M-t^. (1.34) 



This follows directly from the dressed propagator (1.32) and the well-known optical 
theorem, 

r = _^£l m n(M 2 ) , (1.35) 

which relates the decay width of the resonance to the imaginary part of the self-energy 
taken at the value of the particle's mass [131. Z po \ e is some renormalization constant. 
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Expecting ImII(M 2 ) to be small and moreover neglecting the full energy dependence 
of the self-energy, i.e., approximating it by its value on the real axis at the mass of the 
resonance, the optical theorem can be used to write down the propagator in the vicinity 
of the new pole as 



A s (p 2 ) ~ 



Z 



pole 



p 2 _ M 2 _ iZpolc i m n(M 2 ) 

•^polc 

p 2 — M 2 + iMT + 
= Zp °' e - 2 , (1.36) 

where in the second step T 2 /4 was added to the denominator. It should be stressed that 
this is only a good approximation for our assumption r <C M. In fact, the search for 
such poles will be the main part of this work, though we will realize soon that for broad 
resonances the full analytic structure and energy dependence of the self-energy has to 
be taken into account. This will cause the resonance pole to vanish from the complex 
p 2 -plane. 

The corresponding spectral function as the negative imaginary part of the propagator 
right above the real axis becomes 



p(p 2 ) ~ 



Z pole Imn(M 2 



(p 2 - M 2 ) 2 + (Z pole Imn(M 2 )) 2 

Mr (1 37) 

(p 2 — m 2 ) 2 + M 2 r 2 ' K ' 



known as the relativistic Breit-Wigner distribution for a resonance with mass M and 
decay width T. Note that the presented parameterization is strictly valid only for T <C 
M, i.e., long-lived unstable particles and sharp resonances, respectively. For example, 
positively charged pions with dominant leptonic decay channel tt + — > [i^v^ and a mass 
of about 140 MeV possess a mean life time of about 10~ 8 s, which gives a ratio T/M ~ 
1CP 16 , while for neutral pions with 7r° — > 77, a mean life time of about 10 -17 s and a 
mass of about 135 MeV we obtain T/M ~ 1CT 8 . The decay process of such particles can 



be illustrated like in Fig. 1.2 



On the contrary, very short-lived particles cannot be directly observed, yet it is possible 
to establish their existence from a scattering process A + B — > C + D of two incoming 
(stable) particles A and B, and a set of outgoing particles C + D, where the subset C 
has an intermediate resonance S such that S — > C without detection. From literature we 
know that the general expression for the decay width (see chapter 4) has nearly the same 
formal structure as the differential cross section da [13||19| . After performing a scattering 
experiment and measuring the invariant mass distribution of the particles in the subset 
C produced during the interaction, it is possible to find a peak in the differential cross 
section located around a value y/s = M, the mass of the resonance S. If all our above 
simplifications can be justified, namely if T is not too large, the obtained curve in the 
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other particles 




Production of S 



Decay of S 
and detection 



Figure 1.2: Unstable particle production and detection: The production and decay of the 
particle S occurs in a localized time interval such that, comparing to its mean 
life time, there is enough time for S to propagate. 



region of yfs = M will be narrow and can be approximated with good accuracy by a 
Breit-Wigner distribution with a full width at half maximum of about T, the decay 
width of the resonance S (this clarifies why one usually speaks of 'broad' and 'narrow' 
resonances). The differential cross section for the total process A + B — > C + D is then 



proportional to the spectral function from Eq. (1.37), 



MT 

MA + B^C + D)~ (p2 _ M2)2 + M2r2 , (1-38) 

so measurements of the cross section can look for resonances predicted by theory [15| . 
Note that we have only treated a non-interfering production cross section of a single 
resonant state with two incoming (stable) particles, while usually background reactions 
and other multi-channel effects distort the pure contribution from the resonance S, such 
that it is harder to observe if there really is something or not. 

1.3.2 The problematic scalar sector 

The original quark model based on SU (3) flavour symmetry is not the full story of par- 
ticle physics. A very famous sector of particles is the scalar mesonic sector. The scalar 
mesons have bothered the physical community for a long time because their identifica- 
tion and clear evidence are far from having a simple answer. Especially some of those 
mesons possess large decay widths, several decay channels and a huge background |20|. 
In this work we will focus on the lightest scalar meson, namely the a- meson or /o(500) 
state, which cannot be parameterized by ordinary methods, for example by a simple 
Breit-Wigner distribution. During the past decades a large number of theoretical ap- 
proaches were invented to handle related problems and to extract particle information 
from experimental data, while some effective models for QCD involving chiral symmetry 
treated the cr-meson as existent Various reasons for this assumption can be cited, for 



3 The sigma model was introduced by Gell-Mann and Levy 21 after they referred to a field corre 



sponding to a spinless meson that was introduced by Schwinger [22] . In 1965, J. Schwinger, R. P. 
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example it is often argued that the a-meson is the putative chiral partner of the pion in 
linear sigma models, since the vacuum expectation value of this state is used to imple- 
ment spontaneous chiral symmetry breaking. Furthermore, 7T7r-scattering data seems to 



require the existence of the a- meson in order to fit the data accurately, see Refs. 11 20 
and references therein for further information. The state /o(500) is often interpreted as 
a so-called tetraquark state [23T[25] , a hypothetical mesonic structure first suggested by 



Jaffe 26 



In the 70s, the Particle Data Group (PDG) removed the c-meson from its listings only 
to put it back in the 90s in consequence of new theoretical efforts made by Tornqvist 
and Roos [27] when locating the pole position of /q(500) after fitting scattering data. It 
is also possible to determine the <r-pole by using Roy equations with crossing symmetry, 
analyticity and unitarity - many works going in that direction were published in the last 
years and have found similar results. Only to mention one, Caprini et al. have shown in 
Ref. [28] that the 7T7r-scattering amplitude contains a pole with the quantum numbers 
of /o(500) and calculated its mass and decay width within small uncertainties: 

^s~T e = (441^ 6 - i272j* 2>5 ) MeV . (1.39) 

The PDG recently revised its values for the mass and width estimates on the cr-meson 
[20] , giving a range 

= (400-550) - i(200-350) MeV , (1.40) 
which is a much stronger constraint in comparison with the older window from 2010, 

= (400-1200) - i(250-500) MeV . (1.41) 

Both specifications as well as the pole of Caprini et al. should make clear that the 
/o(500) state is a very broad resonance with a decay width T ~ M. The assumptions 
included in the Breit-Wigner parameterization derived in the last subsection are clearly 
invalid. 

As already mentioned, the problematic scalar sector is so far not fully understood 
and thus it is not suprising that different interpretations and models for scalar mesons 
were and are in use. In 1995, Tornqvist studied the scalar sector by including hadronic 
loop contributions [29] , though a similar idea was already mentioned by Hohler 1 30 1 in 
the late 50s. Following the pedagogical introduction of Boglione and Pennington [3l] , 
one can think of scalar states as belonging to ideally mixed quark multiplets. In the 
simplest case of a single resonance we can assign a free propagator to the mesonic state, 
which becomes dressed if the coupling to a decay channel is non-zero. The corresponding 
field S cannot be considered anymore as an asymptotic state of the Lagrangian, since 
the loop contributions given by strongly coupled hadronic intermediate states dominate 
its dynamics. Because of this strong coupling to intermediate states, the scalar sector 
escapes from the general approach of the naive quark model. Such kind of mechanism 
is called dynamical generation, though this term is not used consistently by all authors 



Feynman and S. Tomonaga were honoured with the Nobel Prize in Physics for their work on quantum 
elekctrodynamics (QED). 
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working in the field, see Ref. [32] for an overview. The mass and width of the resonance 
are now determined by the position of the complex pole of the full interacting propa- 
gator, a procedure first proposed by Peierls [33] a long time ago, making the quantum 
theoretical treatment of unstable particles to become an object of much interest (see 
Refs. 



[30 , 34 - 36] for only some articles published after Peierl's work) . This is usually 
called a pole approach. The power of such a definition of mass and width (of resonances) 
lies in obtaining gauge and field-redefinition invariant values in order to qualify physi- 
cal observables 37 . It is well-known that the Breit-Wigner parameterization does not 
fulfill these properties [38[j42]. It should also be noted that the vast majority of pole 



approaches first model the scattering amplitude by fitting experimental data and then 
extract the values of M and T from the propagator pole. In this work we will only focus 
on the pole position without looking at data, yet there is no problem of extending our 
work in that direction. For a short introduction on how to use the formalism regarding 
scattering amplitudes we highly recommend Refs. 29 31 43 44] (and references therein). 
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2 On multi-valued complex functions 



2.1 General problem 

In ordinary mathematical analysis, the definition of a function is not very sophisticated 
at first glance. It is said to be a mapping of elements x from a domain X to elements y of 
another domain Y, where the latter is called image of the specified function / [45] . One 
can also say that the function / is a relation between X and Y, a set of ordered pairs 
(x,y) such that / is a subset of the cartesian product X x Y [46]. Now, the difference 
between a function and an arbitray mapping or relation is the fact that the first one, 
given domains X and Y, is a so-called unique mapping or relation between X and Y: 
by looking at all ordered pairs of elements, for every x in X there is a unique y in Y for 
which the pair (x, y) is in /. Strictly speaking, if an element x of X is given and it is 
not possible to point at one and only one element y of Y such that the pair (x, y) is in 
/, then / is not a function. The same statement in the context of a mapping becomes: 
each input value x needs to be associated with exactly one output y. 

We already know some examples of relations of which we usually speak of functions 
but where we actually need modifications to be in line with the upper definition. The 
root function 



is well-defined for all positive real numbers including zero, and it is known to be single- 
valued because the mapping is indeed unique. Every ordinary function is single-valued 
(in fact, this was the key point in our definition). Nevertheless, we can expand the 
image of the root function onto all real numbers M. Although the new relation is not 
anymore a function in the spirit of our above definition, it is commonly handled like 
that and simply called a multi-valued function. The main reason for such an ambivalent 
treatment is grounded on the possibility to make / single-valued whenever it is needed 
(so to fix the domain Y in a way that the mapping / becomes unique). The full root 
function can be split into two single well-defined functions and the only difficulty may 
lie in the necessity to invent a notation for distinguishing between those two. In our case 
it is convenient to write 



One can now perform all calculations and mathematical manipulations with respect 
to the rich analytic structure of the full function and make it single-valued right at the 
moment when it is needed to be. In physical practice, for example when solving quadratic 




(2.1) 




(2.2) 
(2.3) 
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equations, negative results are often arising as contributions from the negative branch of 
the root (but are not desired, e.g. because of boundary conditions). They can be wiped 
out by setting iy^ — > +^Jx. The root is then purely single-valued. 

All this is quite intuitive by considering only real numbers. If we now extend our 
domains to complex numbers, a new kind of multi- valued character of functions comes 
into play. To begin with, let us study the complex root function 

/:€->• C, z ^+V^=y/z = w . (2.4) 

We will encode points in the complex plane by using polar coordinates z = pe 11 ^. The 
root function then can be written as 

f(z) = Vt = y /pe i %, 0€(-ir,7r]. (2.5) 

While there is no problem to give the root of a negative real number, the function / as 
denoted in the last equation is not single-valued and therefore not well-defined for all 
z € C. We can convince ourselves of this fact by looking at the complex z- and zu-planes 



(see Fig. 2.1). Imagine we take a path C in the z-plane starting from the point z = —p 
and walk counterclockwise back to that point. In the w-plane this would correspond to 
a path starting from —i^fp and ending at iyfp, which is a semicircle on the right side of 
the imaginary axis. So, coming back to the starting point in the z-plane does not give us 
the same value for yfz = w. The complex root function is obviously not single-valued: 
we must take the same path as before in the z- plane, denoted as C , in order to arrive 
at the same point in the u>-plane. In addition to this geometric view, we can investigate 
the behaviour of / by approaching the negative real axis from the upper and lower side 
of the complex z-plane: 

lim f{—p + ie) = y/pe 1 ? 

= iyfp, (2.6) 

lim f(-p - ie) = Jpe~ % ^ 
t— >o+ 

= -VP- ( 2 -7) 

By coming from different directions the limits of the root function on the negative real 
axis are in fact not the same. It should be clear now that there is something strange 
going on in this region and one may ask how we could manipulate / to get rid of its 
multi-valued character. Usually, this is done by taking a 'mathematical scissor' and 
making a cut along the negative real axis of the z-plane, starting from complex infinity 
on the left and proceeding down to the origin. It is important to understand the cut 
only separating the upper and lower half planes on the left - in our definition of the 
root function all the points on the respective axis are placed above the cut. If we now 
take the same paths C and C in this new sliced z-plane, there will be no problem by 
approaching the negative real axis from above because there is an edge that cannot be 



overcome. With this procedure we can rewrite our definition (2.4) in the following way: 
/:C^C\{weC:Ret«<OV(Ren) = OAlniw<0)}, z i-> = = w . (2.8) 

This is called the principal branch of the complex root function. 
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Figure 2.1: Multi-valued character of /(z) = s/z with paths C (black, dashed) and C (red, 
dashed) in the complex z- and w-planes. 



2.2 Branch points and branch cuts 

The example just illustrated is not as special as it might seem. There exists a lot of 
common functions which have a similar analytic structure as the complex root, or an 
even more complicated one. For another famous example let us look at the complex 
logarithm 

/ : C -)• C\{w G C : -7T < Imio < vr}, z^lnz = w . (2.9) 

This function has also a cut along the negative real axis (this time including the origin) . 
We should remark as a general statement that the cut is not a unique entity but can 
be any arbitrary line (or line segment) connecting two so-called branch points [47| . In 
the case of the root function and the logarithm the branch points are situated at z = 
and z = oo, while the cut is taken to be a straight line connecting them. Nevertheless, 
we could also have chosen the positive real axis along which to be cutted along (or any 



17 



other, even curved line) and in practice one should take the one that makes calculations 
as simple as possible^] In this manner all problems of considering multi-valued functions 
simply reduce to the study of their branch points. 

It was shown in the previous section that by circling around the branch point at the 
origin the value of the root function is different from the value with which we have 



started (see Fig. 2.1). This can be taken as the definition of a branch point. There 
would have been no difference if we took any other contour path in the z-plane than a 
simple circle; as was already stated, only the branch points are the independent objects 
which cannot be moved to another position. Consider the complex logarithm / from Eq. 



(2.9) in the whole complex plane now. By taking again a circle z = pe 1 ^ with (ft = —it. ..it 



around the branch point at z = 0, we find: 



lim f(—p + ie) = In pe l 



= In p + in 

= In p - iir + i2tt , (2.10) 
lim f(—p — ie) = \npe~ m 

e-S>0+ 

= hip -iir. (2.11) 

In contrast to the root function, the multi-valued character of the logarithm appears 
as an additional imaginary part of 2tt. It can also be shown very easily that it is not 
possible to come back to the value of In z at the starting point z = —p — ie by crossing 
the negative real axis, since every additional circle around the branch point adds another 
2tt. As a consequence, the cut leads to different branches of the logarithm, exactly as 
in the case of the root function (where only two branches exist). So it is natural to call 
such objects branch cuts from now on. The case without imaginary shift is taken to be 
the principal branch of the logarithm and one can define the kth branch by 

\n(z;k) = ln\z\ + i(arg z + 2A;vr) , kZ, (2.12) 

where the domain for every single branch is assigned in the same way as was given in 



the naive definition (2.9). 

We motivated branch points and branch cuts in order to obtain a single- valued complex 
function out from a multi- valued one. The key point consisted in the 'exclusion' of the 
cut region from the domain of the function in such a way that the cut acts like an edge 



in the complex plane. Note that there is no exclusion at all, see also Fig. 2.2 for a better 
understanding of this subtle but important detail. A huge amount of literature about 
mathematics for physicists is dealing with this kind of description of branch cuts. As will 
be seen in the next section, this is not enough to describe physics of scalar resonances, 
because the upcoming relevant functions to be analyzed in this field (especially the 
Feynman propagator of interacting quantum field theories) will also have branch cuts in 



1 In the mathematical literature and computational software, branch cuts are chosen by convention to 
give simple analytic properties. In this work we adopt the cut structure of multi- valued functions as 
is used by Mathematica. 
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Figure 2.2: Understanding the branch cut: In this picture the points that are on the real 
axis are placed above the cut (indicated by a solid line) , so to be situated right 
above the axis. Branch cuts of a specific function can be also defined the other 
way around. 



the complex energy plane. In order to understand the physical behaviour of resonances 
we need to get a deeper insight of what is going on at the cut. In fact, we need to go 
beyond it. 

2.3 Introducing Riemann sheets 
2.3.1 Motivation and definition 

We have seen how a multi- valued complex function can be made well-defined by cutting 
the complex z-plane along a line (or line segment) that effectively behaves like an edge for 
every path taken from one side to the other. Within this framework it is not possible to 
cross the branch cut - so we could in principle grab both separate layers of the complex 
plane and move them in the direction of the imaginary axis to visualize the discontinuity 
of the function / in the tw-plane. In case of the root function, for someone coming from 
below the cut this deformation would look like a wall, while it would look like an abyss 



for someone standing on the other side (see Fig. 2.3). By approaching the branch cut 
the difference of w between both sides is only in the imaginary part. 

We now want to abandon our above handling of branch cuts and introduce not only 
another, but also a more general interpretation of multi-valued complex functions. In 
fact, we already have worked out all necessary aspects of this new way of thinking on the 
last pages, we just need to merge them differently from how we did before. Remembering 



Fig. 2.1, we notice that the two sketched z-planes are actually the same, although the 
value of it; in both planes is not. It seems that the problem of obtaining a single- valued 
function is only a question of choosing one of those z-planes. This was also suggested 
when we defined the kth branch of the complex logarithm: every single value of k 
corresponds to one specific z-plane where the branch lives in. The same could have been 
done in the case of the root function, yielding only two distinct complex planes. So, 
let us exploit the idea of having really two z-planes. Our complex function / is then a 
mapping from those two planes onto a single -u;-plane and the distinction between the 
two planes is made by giving the fc-value. It is fully natural to take this value as an 
additional coordinate - by giving all three coordinates, namely the real- and imaginary 
part of a given point in the complex z-plane and the k- value of that plane, we can assign 
exactly one value w to every given point z. The function appears to be single-valued. 
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Figure 2.3: Discontinuity of the root function f(z) = y/z = w. 



But what has happened to the branch cut? In our new interpretation the cut can be 
understood as the connection between the two different z-planes because now there is 
no discontinuity at all when slipping from one plane to the other. We can attach both 
planes along the cut such that circling around the origin will make us leave the first 
plane after a full polar angle of 2tt, walking on the second plane and letting us arrive 
at the starting point in the first plane after a (global) turn of 4tt. The structure of the 
two combined z-planes appears as a very simple closed spiral stairway. To make this 



aspect clearer, we visualize the structure in Fig. 2.4 by using a color encoding scheme 
for complex numbers. 

In this new framework any crossing of the branch cut is not special anymore since 
the function considered passes continuously from the first plane to the second one or 
vice versa from the second to the first one, depending on the path taken. We call every 
individual z-plane a Riemann sheet of the Riemann surface. For the root function the 



latter is the mentioned closed spiral starway structure shown in Fig. 2.4 As a result, a 
new definition of a single-valued complex function / is needed here: 

f:X^C,f(z) = w, (2.13) 

where X is a Riemann surface, a one-dimensional complex manifold. We do not want 
to get into the mathematical details of such manifolds, but shall mention the most 
important theorems required in this work. For a deeper mathematical introduction see 



for example the standard textbook by Forster 48 , whereas a very good presentation 
can be found in Ref. [49] covering a wide range of the practical aspects for calculations. 
For just a quick but adequate look we highly recommend Ref. |50|. First of all, it can be 
proven that Stokes' theorem is still valid on a Riemann surface and this leads directly 
to a generalization of the residue theorem. Contour integration along a path, however, 
needs to be defined by using differential forms, but this will not affect us here. Apart 
from that, the very important identity theorem also holds true and therefore we can 
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Figure 2.4: Riemann surface of the complex root function. Each complex value w is repre- 
sented as a particular color: the arg of the complex number is encoded as the 
hue of the color, the modulus as its saturation (the colored background graph- 
ics on the left as well as the figure on the right were created by Jan Homann 
from the University of Pennsylvania) . 



apply analytic continuation not only from the real to the complex, but also from one 
Riemann sheet to the other. Since this technique will be crucial, let us review the basic 
ideas and apply them to Riemann surfaces. 

2.3.2 Analytic continuation 

Suppose we have two holomorphic functions f\ and defined each on domains £7i , O2 C 
C such that both domains have a non-empty intersection S7 = n O2 with fi(z) = f2(z) 
for every z £ f2. Then both functions are analytic in their domains and can be expressed 
in terms of a power series, 

00 00 

fl(z)=J2^n(z-Z 1 ) n , f2{z) = Y,K{Z - Z 2 ) n , (2.14) 
n=0 n=0 
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where z\ 6 $7i and zi 6 $72- Either of the two expressions is certainly valid in the 
intersection region Q and is consequently the series representation of one and the same 
analytic function / around the points Z\ and Z2, respectively 51 . Thus, the function ji 



is called the analytic continuation of f\ onto Q2 and vice versa. 

We can immediately apply this procedure to the complex root function f(z) = \fz 
around the point z\ = 1 on an open disc with radius one. The latter limitation reflects 
the radius of convergence of the power series of /, so the root function seems to be only 
analytic on the disc. Analytic continuation now provides us with a powerful feature. We 
expand / around the point Z2 = i in another open disc, but with the same radius as 



before, and hence there arises an overlapping region, see Fig. 2.5 The expanded function 
is again holomorphic in the second disc and, moreover, it is still the root function! In 
fact, we could take the point Z2 in such a way that it lies in the first disc and where a 
point z = i is included in the new disc resulting from the series expansion around that 
Z2- Either way, we realize that both resulting functions agree in the intersection region 
and are consequently the same in both discs. 
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Figure 2.5: Analytic continuation of the complex root function by expanding it in a power 
series in two different discs (gray) and realizing that both representations equal 
each other in the intersection region (white). 



It is possible to extend the natural domain (i.e., a disc with finite radius of convergence) 
of a single series representation of a function without passing through any singularity 
by adopting the outlined method. For instance, if we add a third disc with center at 
23 = — 1 the root function is extended from the second disc to negative real numbers. 
Note that there is no branch cut preventing us to go to and beyond the negative real axis. 
Instead we would clearly rediscover its necessity when continuing further back to our 
starting point at z = 1, because we would end up with the negative root function. This 
is a very non-trivial statement, though. Remember, the most important result of the 
previous subsection was the introduction of a general procedure to make a multi- valued 
complex function well-defined by introducing a Riemann surface. Usually, physical and 
mathematical literature provides us only with complex functions or (in) definite integrals 
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on the first Riemann sheet. We are finally able to pass through the branch cut by apply- 
ing analytic continuation down into the second sheet where the analytic structure can be 
very different than in the first one. This will be shown in the next chapters when we will 
consider decay processes in quantum mechanics and quantum field theory. In order to 
study these processes, we need to know the explicit form of the propagator on the second 
Riemann sheet. However, we have already seen two examples, the root function and the 
logarithm, where the function values are known on the complete Riemann surface. For 
the first one it is 

f(z) = y^, (2-15) 
fn(z) = -f(z) . (2.16) 

We have denoted the second Riemann sheet with a roman number. The function on 
the first sheet will never have such a number in this work. We got the above result 
by approaching the branch cut on the negative real axis from two different directions, 



compare Eq. (2.6) and (2.7). The observed discontinuity is of course only the difference 
between two sheets and not an intrinsic property of the root function. To gain an 
appropriate continuation of / across the cut into the second Riemann sheet, we first 
need to calculate that mentioned discontinuity at z = —p: 

Disc f{-p) = lim f(-p + ie)-f{-p-ie) 

= *v^-(-*Vp) 

= liyfp. (2.17) 

The analytic continuation of the root function down into the second sheet is then per- 
formed by accepting the requirement that the value fa(—p — ie) (i.e., the function just 
below the cut in the second sheet) equals the value /(— p + ie) (i.e., the function just 
above the cut in the first sheet) along the whole negative real axis: 

lim fu{-p - ie) = lim f(-p + ie) 

e^0+ 6-5.0+ 

= lim, f(-p - ie) + 2i^p 

£—5-0+ 

= -iy/p + 2iy/p 

= ijp. (2.18) 

The analytic extension of this result into the lower half plane is after all what we were 
looking for, 

f u ( z ) = , (2.19) 

and exactly the second branch of the complex root in the context of our first interpre- 
tation of multi-valued functions. We should also notice the very useful identity between 
the imaginary part of /, taken right above the real axis, and its discontinuity: 

Disc f(-p) = 2i lim Im f(-p + ie) . (2.20) 
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One may ask immediately why the extension (2.19) of the purely imaginary result from 
Eq. (2.18) is so simple. In fact, we should be more precise in this point: The identity 



theorem states that if two given functions / and g, both holomorphic on a domain Q, 
equal each other on some line segment U lying in 0, then they equal each other on the 
whole domain Q [52] . Here, we see at first glance what is really the difference between 
real and complex analysis: a holomorphic function defined on Q is completely determined 
by its values on a line segment Pc!J. In this manner, there is no problem of continuing 



Eq. (2.18) into the complex plane. 



For completeness, we quickly calculate the complex logarithm in the second Riemann 



sheet and we should of course obtain the first branch as was given in Eq. (2.11) and in 
(2.12) for k = 1, respectively^] 



Disc/(-p) 



lim 

e->0+ L' 



f(-p + ie)-f(-p-ie) 

In p + in — (In p — in) 
2ni , 



lim fn(-p - ie) = lim f(-p + ie) 



lim /(— p — ie) + 2ni , 

e->0+ 



and so we indeed find for every z G C\{z G M. : z < 0}: 

fll(z) = In z + 2ni . 



(2.21) 
(2.22) 

(2.23) 



2.4 Dispersion relations 

By using branch cuts instead of Riemann sheets one can perform contour integration in 
the complex plane without any problems. This is clear: in that case we simply restrict 
ourselves to a single (sliced) complex plane and the only new thing to deal with is 
the integration along the cut. In the following section we shall first derive the useful 



identity (2.20) between a multi- valued complex function and its discontinuity at the 



branch cut, then we will make use of this relation to obtain the root function only from 
its cut structure. This can be done by performing a single integration with the help of 
Cauchy's integral formula. All this will be important in the next chapters to varify that 
every calculated expression for the propagator indeed is the right one. 

The most important complex functions in this work, like the mentioned self-energy as 
well as the scalar propagator, have in common the powerful property 



f(z) = f*(z*) 



(2.24) 



2 The Riemann surface of the complex logarithm is constructed by attaching the single complex z-planes 
in a slightly different way than in the case of the root function. It can be visualized by an infinite 
spiral stairway. 
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Figure 2.6: Contour path of integration C (red, dashed) enclosing the finite branch cut 
when blowing up the contour C (black, dashed). 



which leads to the useful identity from Eq. (2.20), since 



lim 

e->0+ 



f(x + ie) - f(x-ie) 



f(x + ie)-f*(x + ie) 



lim 

2i lim Im f(x + ie) 

e-S>0+ 



(2.25) 



From this we also obtain that the function / is either purely real on the real axis or has 
a branch cut with the discontinuity just calculated 47 . We now assume the function 



to be holomorphic except at the cut and to vanish faster than 0(l/|,z|' ? ) (with r\ > 0). 
If the cut has finite length then due to Cauchy's integral formula the function can be 
expressed as the sum of two integrals, namely 



/(*) 



1 

2vrl 



c Z,-z Jc< £ 



f(0 \ 



z I 



(2.26) 



where the contour paths of integration are taken as illustrated in Fig. 2.6 This is true 
when we blow up the contour, because in this case the first integral does not contribute 
at all. If we additionally shrink the contour C in such a way that only the paths below 
and above the real axis survive, the function can be represented in the limit of e — > + 
by its imaginary part right above the branch cut: 



/(*) 



1 



lim / dx 



f{x + ie) f{x - ie) 



2iri <e->o+ J a \x + ie — z x — ie — z 



1 

27Ti 

1 

7T 



6 , 2i lim (E ^ >0 + Im/(x + ie) 
ax 



dx 



x — z 

lim e _>. + — Im f(x + ie) 

z — x 



(2.27) 



where a and b mark the branch points on the real axis. The integral representation for 
the Feynman propagator follows directly from this calculation. There is no difficulty 
at all to implement additional simple poles of / on the real axis, since they can be 
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taken into account by deforming the integration contour to make use of the residue 
theorem. Physically, the poles will appear as free propagators with some renormalization 
constants. 

So, by restricting the integration to be performed only in the first Riemann sheet we 
need to be aware of taking the correct limit of the function /, depending on which side 
of the cut the contour path actually lies. Suppose we know that there is a branch cut on 
a part of the real axis and there shall exist an analytic expression for the discontinuity 
across this cut. We now ask for the complex function / with exactly this discontinuity 



and which also fullfils the relation (2.24). From our above considerations it is obvious 



that the whole function is determined only by the discontinuity and can be calculated 



by evaluating the dispersion integral (2.27). As an instructive example we take again 



the complex root function f(z) = \fz. Since it does not even decrease for \z\ — > oo, we 
need to modify the dispersion integral by using a slightly changed function: 

g(z) = ^ . (2.28) 

This new function has the same branch cut structure as / (note that there is no simple 
pole at z = 0). The discontinuity of the pure root function at the cut was found in the 
last section to be simply two times itself. For the new function g this means 

Disc g(-p) = -2i^l . (2.29) 
P 

The dispersion integral can then be computed by using a Hankel contour path of inte- 
gration C with left open end: 

1 f-°° , , -liyff 



r — OQ 

/ v- 

Jo P 



2vri J p'{-p' - z) 

= ^Tdp-^-, (2.30) 
vr Jo p(p-z) 

where in the last step the variable transformation p' — > —p has been introduced. Finally, 
the function / can be denoted as 



p— >x 2 



vr Jo p(p-z) 

2iz f°° , 1 
— / dx k 

IT Jo Z - X Z 

U/Z f°° , / 1 1 

dx — = h 



7T J n \ V z + X \J z — X 

i\fz~ 



7T 

i\fz 



ln(v^ + x)|jj° - ki(*/z - x)|jj° 



(-»7T - 0) 



7T 

rz . (2.31) 
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3 Quantum theory of unstable particles 



3.1 Non-relativistic Lee model 
3.1.1 Lee Hamiltonian 

In 1954, Le^] worked out a general technique to handle, with solvable interacting quan- 
tum field theoretical models, the study of the renormalization problem [53] . In this 
context he introduced a so-called Lee Hamiltonian which can be used to deal with decay 
processes in quantum mechanics. To this end, we need to couple an initial state^] \ S) to 
a continuum of two-particle states |k). In the following, we assume the decay process 
to be a transition of the form IS") —¥ |k) where the unstable particle S decays into two 



particles with momenta k and — k. Following Facchi and Pascazio 55 , we decompose 
the full Hilbert space H of all normalized states of the system in such a way that the 
state 1 5) lives in a subspace Tis- 

'H = 'Hs@'H^- (3.1) 

In order to be more precise, we can write explicitly the set of states living in each 
subspace: 

H S = {\S}}, (3.2) 

7i k = span({|k)}) . (3.3) 

For the moment, the whole system is to be considered in a finite volume V = L 3 and as 
a consequence of applying periodic boundary conditions the momenta k have the form 
k = *y n with n = (m, ri2, ns) T and rtj G Z, i = 1, 2, 3. The full Hamiltonian H shall 
be a sum of a free and an interaction part, namely 

H = H + H- mt , (3.4) 

and we assume the states {|k)} to be the eigenbasis of Hq in the subspace with 
eigenvalues w(k), while \S) shall be the eigenstate in Hs with eigenvalue wg. It is 



obvious that because of our decomposition (3.1 ) of the Hilbert space the states belonging 



1 In 1957, T. D. Lee and C. N. Yang were honoured with the Nobel Prize in Physics for their work on the 
violation of the parity conservation law in weak interaction, which C. S. Wu verified experimentally. 
Lee was just 31 years old. 

2 As was pointed out by Fonda, Ghirardi and Rimini [54], it is a controversial problem to assign a state 
vector to an unstable system, so we will call it simply 'initial' state without further commenting on 
this point. 



27 



to either subspace are ortho-normal to each other. Hence, we can define two operators 
Ps and Pk that project arbitrary states of % onto Us and 7i k : 

P S = \S)(S\, P k = t-\S)(S\. (3.5) 

Although we do not know the eigenbasis of the full Hamiltonian we nevertheless can 
work with the eigenbasis of the free part and write down a completeness relation 

l = |S><S| + 5>)<k| . (3.6) 

k 

The Lee Hamiltonian can now be constructed for the purpose to deal with our decay 
transition problem. The interaction part will mediate between states of the two sub- 
spaces Us and Hk, respectively, while the free part of course must deal with states from 
both: 

H = u, 5 |S)<S| + J>(k)|k)(k| , (3.7) 

k 

^ = E^(|k>(S| + |S)<k|) . (3.8) 
These expressions arise from 

H = PsHPs + P^HP^, (3.9) 
flint = P s Hp k + P k Hp s . (3.10) 

Thus, the interaction part operates only between the initial state \S) and the final states 
|k) (it is completely off-diagonal in the given basis), while the free part simply gives 
their eigenenergies and leaves a considered state untouched. Technically speaking, the 
main structure of our Hilbert space and the Hamiltonian is in full analogy with that of a 
two-state quantum system, commonly covered in an undergraduate quantum mechanics 
course. The only difference is that there is a set of final states |k) instead of just a 
single one and the interaction term to be time- independent, governed by the form factor 
function /(k) with explicit dependence on the momentum k of the final state. 



3.1.2 Schrodinger propagator 

As will be emphasized soon, the main starting point for our analysis in general is the 
propagator of a specific physical system, for now the non-relativistic Schrodinger prop- 
agator of the decaying particle S: 

G S (E) = (S\ I . \S) . (3.11) 
h — H + te 

It should be mentioned that there may be some confusion by looking into traditional 
literature due to a subtle difference in terminology. In quantum mechanics, referring to 
the term 'propagator' in the context of Feynman's path integral formalism, we address 
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the (probability) amplitude for a particle to travel from one spatial point at one time to 
another spatial point at a later time. One usually derives the propagator, often denoted 
as K(x, t, x', t'), and further the energy-dependent Green's function G(x, x', E), which is 
the Fourier transform of the propagator with respect to timej^] In this work, as done in 
most of the literature concerning (relativistic) decay processes, we will name both objects 
in the opposite way. In particular, Gs(E) shall be the propagator for the particle S. 
To calculate the propagator we first take a look at the survival amplitude in the initial 
state | S), 

a(t) = (S\U(t)\S) = (S\e- iAt \S) , (3.12) 
~ lHt is the time evolution operator. In Ref. 58 , Moshinsky et al. demon- 



where U(t) = e" 

strated for time-independent interaction terms in one dimension that the Schrodinger 
propagator (the Green's function) in position space can be obtained from the above 
mentioned amplitude K(x,t, x' ,t') by a Laplace transformation. In fact, this makes 
sense because the continuous Fourier transform is equivalent to evaluating the bilateral 
Laplace transform with imaginary arguments. Applying a similiar procedure to our 



problem, we take the Laplace transform of Eq. (3.12) and simply get the resolvent of 
U(t) by applying the geometric series: 



F(s) 



(S\ 



(S\ 



dt e 



dt e 



- st e- iHt \S) 



n=0 



(-*)' 



n! 



l \S) 



(SI^HrW die 



-at ' 



nl 



:\S) 



(S\ 



n=0 

s + iH 



H r 



>-iE 



(S\ 



E-H 



\S) . 



(3.13) 



This in principle gives the expression in Eq. (3.11), yet we do not stop here, especially 



without commenting on the variable transformation in the last step. Let us first recover 
our survival amplitude by applying the inverse Laplace transformation: 

f 7+ioo 



a(t) 



-iE 



(S\ 
i 

2^ 



1 
1 



( dse st 
7— zoo s + iH 

17—00 



2vri 



dE e 



-iEt 



1 



17+00 



+00 



dE e 



-iEt 



(S\ 



E-H 
1 



I S) 



E-H + ie 



\S) . 



3 A good introduction into path integrals in quantum mechanics can be found in Ref. 

detailed study one should consult classic literature on this topic, e.g. the textbook by Feynman and 



(3.14) 



For a more 



Hibbs 57. 
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The integration in the complex s-plane has to be done along a vertical line with Re s = 7 
such that 7 is greater than the real part of all singularities of F(s). More general, 7 is 
chosen in a way that the contour path of integration lies in the region of convergence 
of F(s). In order to guarantee real values for the energy variable E, the propagator 
from Eq. (3.11) will have either simple poles or a branch cut (starting and ending at 



some branch points) on the imaginary axis in the s-plane, depending on the coupling 



g in the interaction part (3.8) of the Hamiltonian. The integration path needs to be 
shifted by an infinitesimal number 7 = e > to the right. The s-plane is then rotated 
by the transformation s — > —iE, see Fig. 3.1 Thus, the propagator becomes analytic 



in the whole new physical i?-plane, except at the simple pole positions and the branch 
cut, respectively. These last statements are however not obvious and so the full analytic 
structure of the propagator will be studied in the next subsection. 



Li 



VWWWWSAA * 



Figure 3.1: Complex plane rotation s 
path displayed. 



— iE with branch cut and contour integration 



We are now able to determine the Schrodinger propagator explicitly. By following 
Giacosa 
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one can rewrite the propagator using the decomposition H = Hq + H u 



and performing a sequential insertion of completeness relations (3.6): 

1 



G S (E) 



(S\ 



(S\ 



(S\ 



(S\ 



E — 


* 

H + ie 
1 


E — 


Ho + ie 
1 


E — 


Ho + ie 
1 


E — 


H Q + ie 
1 



IS) 



El* 



n=0 



E — Hn + ie 



\S) 



1 + l# int l 



1 

E- Ho + ie 



+ ... )\S) 



E- Ho + ie 



1 + (\S)(S\ + £ |q)(q|)iI int (jS)(S| + £ |k)(k|) 

\S) . (3.15) 



+ . . . 
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We only have to evaluate matrix elements of the following form: 



(S\ 1 1 5) = 1 , (3.16) 

(q| 1 |k> = = ^ , (3.17) 

E-H + ie E-oj(k) + ie V ' 

(S\—!r Hk> = 0, (3.18) 

E — Hq + te 



and since the interaction part is completely off-diagonal we find 



(S\H int \S) = <q|ffint|k> = 0, (3.19) 
(S\H int \k) = (k\H- mt \S)* = g 1 ^ . (3.20) 



All terms in the infinite sum over n can be collected in a resummation after taking all 
the Kronecker deltas into account. We are left with: 

c (E) 1 1 V 1 g2f2{k) 



E — ojq + ie \E — ujq + if- V E — u)(k) + ie 

n=0 k 

-g 2 nE) 

E - us + ie \E - ojs + ie 

n=0 

1 



E - u s + g 2 ^(E) + ie ' 
where in the last steps we have defined 



(3.21) 



S (£) = _yl / 2 ( k ) . (3.22) 



V E-Lu(k) + ie ' 



Finally, performing the continuum limit we get 



(3 - 23) 



In the same way as was shown in section |2.4[ one immediately obtains an integral 
representation for the Schrodinger propagator by applying Cauchy's integral formula. 
This spectral representation can be also derived from a physical point of view in the 
context of many-body physics, similar as we did in the first chapter for relativistic 
quantum field theory (see for example Ref. [59] ). In the end the propagator reads: 



G S (E) = I ^ du -fM— , (3.25) 

with the spectral function 

p(u) = -lmG s (u + ie) . (3.26) 



31 



3.1.3 Analytic structure of the propagator 



We will now study a decay process within a specified Lee model. We restrict ourselves 
to final states |k) with eigenvalue w(k) > with respect to the free part Ho of the full 
Hamiltonian H. The decaying particle S will be prepared in its rest frame with energy 
us = Mq and the form factor function / (k) shall be a simple product of two Heaviside 
step functions of the form 



/(k) = 0(2k - E )Q(A - 2k) 



(3.27) 



where k = |k|. The parameter Eo is nothing else than the energy threshold; we thus 
assume Eq < Md < A, where A is the upper bound of the spectrum, acting as a cutoff 
(one could understand 6 (A — 2k) as a cutoff- function) . We chose the final energies as 
w(k) = 2k. Although this model represents an unrealistic simplification, the following 
section will provide us with the basic mathematical framework for our quantum field 
theoretical investigation. 

there exists a remarkable connection between the 



As was discussed in section 2.4 



branch cut structure of some complex functions and dispersion relation integrals of the 
generic form 

p(x) 



dx 



Z £ 



(3.28) 



We have also figured out in chapter 1 that whenever the (analytic) spectral function 
p(x) gives support to the integral, there is a branch cut arising in this region on the real 
axisj^] Taking a look at the relation between the non-relativistic propagator Gs(E) of the 
particle S and the spectral function p{u) in Eq. (3.25 ), we conclude that there is a branch 
cut if the spectral function is not a delta distribution function - i.e., if the coupling g in 



our Lee Hamiltonian is non-zero. Then the self-energy (3.24) will contribute to the full 



propagator and influence the position of the mass-pole in the complex E'-plane, where 
E —> z = x + iy. In one dimension, the self-energy takes the simple form 



E(E) 



dk 



/ 2 (k) 



2ttE - w(k) + ie 
J_ r°° dk Q{2k - E )O(A - 2k) 
2vr L E-2k + ie 



A 



2k^k 1 

1 

47T 



dk 



1 



k-E-ie 
ln(A — E — ie) — ln(E 



E-ie) 



(3.29) 



There is a very important aspect here to point out: In the last step we actually performed 
the limit e — > + . From a conservative point of view we should make use of the Sokhotski- 
Plemelj theorem, according to which an additional discontinuity of in would come into 



4 To be more precise: The function p(x) gives us two branch points on the real axis. The branch cut can 
be any arbitrary line connecting those two points, but in order to keep our definition of the complex 



logarithm in Eq. (2.9 1 and (2.121 we have to draw the cut along the real axis. 
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play. However, this term does not appear in our expression. In the complex Emplane the 
self-energy has a branch cut on the positive real axis starting from z = Eq and ending at 
z = A, so as long as Imz is finite, T,(z) is well-defined. Whereas, if we try to interpret 
the limit e — > + as a falling of the pole z = E + ie onto the real fc-axis, we will have 
a problem with our formalism. But this problem is actually not a real one: except for 
the spectral function, we will never evaluate the self-energy at such points, since for 
g 7^ the pole will be a complex number with non-zero imaginary part and for g = 



the self-energy will vanish in the propagator, see Eq. (3.21). The infinitesimal number 
ie is only necessary if one strikes real z = E. All this makes the expression (3.29) to live 



on the first Riemann sheet. The ^-integration should be taken with complex energies 
z = x + iy right from the beginning: 



f dk 

Je 



k-E 
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Je k 



ln(A 



ln(E 



(3.30) 



The last result is the final one because common multiplication rules for the logarithm 
are in general lost for complex arguments. A plot of the self-energy for the parameters 



(3.35) can be found in Fig. 3.2 In order to obtain the correct plot one needs to evaluate 



the function at z = E + ie, including explicitly a small imaginary number ie (or use the 
Sokhotski-Plemelj theorem and calculate the principal value). 




Figure 3.2: Real (dark blue) and imaginary part (red) of the self-energy T,(E) on the real 
axix with the choice A = IO-Eq, Mq = 2Eq and Eq = 1. 



We have already indicated that the self-energy has a branch cut on part of the positive 
real axis. This follows directly from our general discussion in the second chapter, because 
we are dealing with a complex logarithmic function with branch points at z = Eq and 
z = A. The point at complex infinity is not a branch point, which can be easily seen by 
taking a path z = a + pe^ with <fi = 0...2-7T that encircles both points z = Eq and z = A 
(so a needs to be in between both points on the real axis and p must be greater than the 
distance to the farthest one) : after a single full turn of 2tt we arrive again at the starting 
point. One may argue that both single logarithms in £(z) provide us with two branch 
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cuts, one starting at z = Eq and the other at z = A, each heading to infinity along the 
positive real axis. But the discontinuities collected for z > A cancel each other due to 
the minus sign. Note that this is the reason why the situation at the remaining branch 



cut is inverse to what we sketched in Fig. 2.2 This is also set up by the computational 



software Mathematica, though it will be shown in subsection 3.2.5 that the self-energy is 
actually continuous with the first quadrant (i.e., the real axis is situated right above the 
cut). Besides that, one can also investigate the behaviour of the point at infinity when 
making the variable transformation z — )■ 1/z and circling around the new origin. Here, 
too, a full circle does not give any change. 

Let us concentrate on the Schrodinger propagator which can be constructed with the 
specified self-energy. For the simplest case g = we have a free particle with mass 
us = Mq in its rest frame: 

G S (E) = . (3.31) 

^ v ' E-Mo + ie v ; 

We apply here our definition of the (real) mass of a free particle as the pole of the 
propagator in the complex E'-plane (E — > z = x + iy). The pole position of course gives 
the mass slightly shifted into the lower half plane to ensure convergence: 



E, 



pole 



(3.32) 



After turning on the coupling, g > 0, there is a contribution from the (complex-valued) 
self-energy function to the pole position. One could rephrase this statement in the 
following way: The appearence of an additional function in the denominator of the 
propagator forces the pole to move away from the real axis down into the lower half 
plane, which makes the self-energy per se to be evaluated at complex values. The new 
pole position as the zero of the denominator is determined by a system of two equations: 



x- M + g 2 Re^(x + iy) = 
y + g 2 Im Y,(x + iy) = 



We could decompose the self-energy for Eq < x < A in the following way: 

1 



T,(z = x + iy) 



ln(A — x — iy) — Ib.(Eq — x — iy) 



1 ln ( V(A-x)* + y^ 
a [^(E -x) 2 + y\ 

-y 



4tt 

— arctan 



+ 



47T 



arctan 



A 



Eq 



(3.33) 



(3.34) 



where we have used general conversion rules for complex numbers with argz = (p G 
(— 7r, 7r]. But the system of equations, combined with the decomposition of the self- 
energy, cannot be solved with analytic methods. We have to find the solution numerically 
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by performing a direct search for the pole of the propagator in the complex Emplane. 
The parameters are chosen to be 



A = 10£ n 



M = 2Et 



o 



En 



1 



(3.35) 



One should not be surprised to find the pole at z = Mn — ie for vanishing coupling, yet 
for small g there is no pole in the whole complex plane. We may get some numerical 
results near the previous pole position, but inserting those points inside the propagator 
reveals that they in fact are ordinary points. This is in full agreement with our discussion 
in the first two chapters: the complex mass-pole leaves the first Riemann sheet when 
we turn on the coupling! Our Lee model describes by construction the decay of the 
unstable particle S and consequently the particle should be represented as a pole with 
non- vanishing imaginary part on the second Riemann sheet j^] We have to look for the 
pole there. 

For g = the propagator Gs(z) is a single- valued function, while for non- vanishing 
coupling the multi-valued self-energy makes the propagator also not well-defined along 
the branch cut starting from z = Eq and ending at z = A. The cut structure is 'trans- 
mitted' to Gs(z). Although we can guess that the pure discontinuity of the logarith- 
mic self-energy across the branch cut will be i, we can also calculate it directly for 
E < E < A: 

DiscS(E) = £(£ + ie) - £(£ - ie) 

' f ln(A- {E + ie)) -ln(E - (E + ie)) 



4tt 



In (A - (E - ie)) + In (E - (E - ie)) 



1 

47T 



In 



V(A - E) 2 + e 2 
y/(E -Ef + e\ 

In ( V(*-E) 2 + t 2 \ 
\^{E -EY + e 2 ) 



l<p! - l(p 2 



i(p 3 + i(j) A 



1 

47T 

i 

47T 



H01 - ' 
arctan 



+ ( arctan ; , 

V Eq — E 



—-( 

A-E V 

e 



arctan 



En — E 



7T 



arctan 



A-E 



47T 



(vr + 7r) 



i 
2 



(3.36) 



3 It was actually proven by Aramaki and Osawa in Ref. 35 that in the general Lee model the propagator 



of the unstable particle must have poles in the second Riemann sheet. 
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Figure 3.3: Discontinuity of the self-energy (in units of i): numerical difference (dark blue) 
and analytic result (red). 



where some conversion relations for the polar angles have been used. We have plotted 



the discontinuity and the numerical difference of S(z) across the real axis in Fig. 3.3 



Surely, this result can be obtained quicker by making use of the integral representation 
of the self-energy and the Lorentzian form of the delta distribution function: 



DiscE(£) = H(E + ie) - E(£ - ie 

dk 



l r A „ / : i 



4-7T J Eq \k — E — ie k — E + ie 
1 ^d, 2 - 



inJ Eo (k-EY + e* 
rA 

i 

2 ' 



I dk iTi5{k - E) 



IE 

E < E < A , (3.37) 



from which we can observe immediately a discontinuous behaviour only at the branch 
cut. The analytic continuation down into (at least the fourth quadrant of) the second 
Riemann sheet therefore reads: 

£ n (z) = £(*) + \ , (3.38) 
and as a consequence the propagator can be written there as 



G\\(z, , f 



1 

M + g 2 Z u (z) 



(3.39) 



Af + 5 2 (S(z) + |) 
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To make sure we have found the correct expression in the first sheet, we proceed in the 
same way as we did in section 2.4 when the complex root function was computed from 
its cut structure and the requirement f(z) = f*(z*). First, we define the functions 

AC*) = -^MA-*), 9i(z) = A ^l 1 , (3.40) 
f 2 (z) = L H E Q -z), g 2 (z) = hiz) , (3.41) 

47T tjQ — Z — 1 

so that S(z) = fi(z) — fz(z). Cauchy's integral formula for each of the functions gi(z) 
and g2(z) then leads us to the correct result, for example: 



2m Jc £ 
A-z-l 



oo 

d.T 



Am J A ' (A - x - l)(x 

A-z-i r°° , i ( \ 

ax 



47T Ik -A + z + 1 U+l-A 



= ^(ln(x + l-A)|^-ln(x-z)|' 

= i-ln(A-z), (3.42) 

where we have used the partial fraction decomposition method. Beware of the fact that 
the discontinuity of gi(z) is the negative imaginary unit —i/2, because of the minus sign 
in the argument of the logarithm. 

Now, to find poles on the second sheet the system of equations (3.33) is modified to 

x - M + g 2 ReY,(x + iy) = 0, 

a 2 i 

y + 5 2 ImS(x + iy) + y = 0. (3.43) 

In addition to solving this system numerically, we will also state results from an analytic 
investigation of the relevant expressions in order to support all results. Before that, we 
want to point out a subtle difficulty. In literature concerning the search for propagator 
poles in the second Riemann sheet, almost always the analytic structure of the inverse 
propagator is studied (e.g. in Refs. 15,60l). A remarkably good and probably the 



most helpful presentation of the topic can be found in Ref. [15] (but beware of different 
sign conventions!). It is said that the difference of the inverse propagator at the branch 
cut gives the discontinuity to be used for constructing the inverse propagator on the 
second Riemann sheet. This is of course true since by that procedure one just gets the 
discontinuity of the self-energy: 

Disc Gg 1 (E) = Gs l {E + ie)-G s 1 {E-ie) 

= E-M + g 2 Y,(E + ie) + ie - (E - M + g 2 T,(E - ie) - ie) 

= lie + g 2 ! (£(E + ie) - E(E - ie)) 

= g 2 DiscT l (E) (3.44) 
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So instead of calculating the discontinuity of the self-energy, one can use the propagator 
as the multi-valued function. We know that for g ^ the propagator gains the same 
branch cut in the complex .E-plane but if we want to continue it down into the second 
Riemann sheet we need to find its discontinuity which of course will look different: 



DiscGs(-E) 



G s {E + ie)-G s (E-ie 

7T 



1 



-oo 
oo 



p{u) 



E 



uj + ie E — oj 
-2ie 



do; p(u>, . , , ,, , 
2i r °° 

IT 

-2ip(E) 



doj p{oj)i:5(E — oj) 
= 2iImG s (E + ie) . 



(3.45) 



The last step is due to the definition of the spectral function as the negative imaginary 



part of the propagator. We recognize the useful identity from Eq. (2.20) and remember 



the imaginary part to be evaluated right above the branch cut. The (preliminary) 



spectral function of our Lee model is shown in Fig. 3.4 for an arbitrarily chosen value of 
the coupling g. One can now carry out the analytic continuation of the spectral function 



P(E) 
25 - 



2 4 6 8 10 12 14 

Figure 3.4: Spectral function p(E) of the considered Lee model for g = 0.4. 



into the lower half plane, i.e., by replacing it by the imaginary part of Gs(E + ie): 

Gn(z) = G s (z) - 2ip(z) . (3.46) 

But care is needed here: Computational software in principle does not know how to 
perform the correct extension of the spectral function into the complex plane, because 
we either possess only a numerical expression (and the extension of such is often non- 
trivial) or use the imaginary part of the propagator. Nevertheless, we cannot simply take 
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the latter and substitute E —> z. If we first continue the propagator into the complex 
and afterwards take the imaginary part, we clearly see that Eq. (3.46) will not give us 
any pole since the resulting denominator tells us just to solve the old system of equations 
Lalytic 

G s (z) 



(3.33). The full analytic information never enters the final result: 

1 



M + 5 2 S(z) 

1 



M + g 2 ReS(z) + i(y + g 2 ImE(z)) 

-M + g 2 ReZ(z)-i(y + g 2 ImZ(z)) 
- M + g 2 ReZ(z)) 2 + (y + g 2 ImZ(z)) 2 ' 

y + g 2 ImS(z) 
(x-M + g 2 ReS(z)) 2 + (y + g 2 ImS(z))' 



(3.47) 



In fact, we first need to take the imaginary part of the propagator right above the cut 
and thereafter extend into the complex plane (this is the appropriate way to gain an 
expression for the spectral function). All this should at least make us very careful when 
dealing with analytic continuations at the level of the propagator, but mostly it tells us 
to proceed in another way. If we really want to work with the discontinuity of Gs(E) 
we can make use of the self-energy: 

AG(z) = G u (z)-G s (z) 

1 1 



M + 5 2 (S(z) + i) z-M +g 2 Z(z) 
-ig 2 



(3.49) 



[z-M +g 2 (Z(z) + ^)}(z-M + g 2 Z(z)) ' 

The propagator in the second sheet then reads: 

Gu(z) = G s (z) + AG(z) , (3.50) 

which is fully consistent, because on the one hand we get the same poles (from the 
denominator of AG(z)) and moreover we preserve a continuous transition from the first 
Riemann sheet to the second one, 

AG(E) = Gn{E -ie) -G S {E + ie) 

1 1 



(E - ie) -M + g 2 (Z(E - ie) + (E + ie) - M + g 2 Z(E + ie) 

1 1 
E-M + g 2 Y,{E + ie) - ie ~ E — Mq + g 2 Y>(E + ie) + ie 

= , (3.51) 
by taking e — > + in the last step. 
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3.2 Lee model poles 



3.2.1 Couplings g G [0.1,0.4] 



In the following subsections we search the complex E'-plane, using the parameters (3.35 ), 
for numerical solutions of the equation 

z - M + g^niz) = (3.52) 

by varying the coupling g in steps of Ag. In our case the numerics are very sensitive to 
variations in g, so we need to give a condition for deciding whether a found solution z po \ c 
marks a singularity or not (on a numerical level we take it as a pole if Abs G^ 1 (z po \ c ) < 
10~ 10 ). We may also examine the spectral function p(E) as the negative imaginary part 
of the propagator. In fact, it seems normalized only for couplings in the approximate 



interval [0.1, 1.2], see Fig. |3.5| (We will show in section 3.3 that it is indeed normalized 



over the full range of g.) Note that if g tends to zero the spectral function becomes a 




Figure 3.5: Incomplete normalization of the spectral function p(u) for different values of 
the coupling g. 



delta distribution function due to the fact that the pole approaches the real axis: 

P (E) = e + g 2 In^(E + i e) 

(E - M + g 2 Re E(E + ie)) + (e + g 2 ImS(E + ie)) 
cm>o e 

(E - M ) 2 + e 2 

= tt8(E - M ) . (3.53) 

The numerical integration i J^L, dw p(uS) is then not possible anymore (so pure poles 
in the vicinity of the integration contour cannot be considered). In general, p{oS) gets 
broader for larger couplings, and this effect is caused by the complex mass-pole moving 
deeper into the complex plane beyond the branch cut, i.e., the second Riemann sheet: 
the spectral function 'feels' this pole less and so shrinks very quickly the more the pole 
descends. 
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1 


.3634 


10 


,0085 


-0.3445 


9.9914 


10.0085 


2.3522 


9.9914 


2 


.1243 


1 


.0000 


4 


-0.4870 


1 


.3902 


10 


,0166 


-0.4870 


9.9830 


10.0166 


2.5350 


9.9830 


2 


.1514 


1 


.0000 


4.2 


-0.6312 


1 


.4202 


10 


,0296 


-0.6312 


9.9693 


10.0296 


2.7402 


9.9693 


2. 


.1573 


1 


.0000 


4.4 


-0.7769 


1 


.4519 


10. 


,0487 


-0.7769 


9.9485 


10.0487 


2.9732 


9.9485 


2. 


.1410 


1 


.0000 


4.6 


-0.9240 


1 


.4844 


10. 


,0750 


-0.9240 


9.9191 


10.0750 


3.2427 


9.9191 


2. 


.1037 


1 


.0000 


4.8 


-1.0724 


1 


.5167 


10. 


,1093 


-1.0724 


9.8792 


10.1093 


3.5638 


9.8792 


2 


.0485 


1 


.0000 


5 


-1.2219 


1 


.5480 


10. 


,1520 


-1.2219 


9.8273 


10.1520 


3.9686 


9.8273 


1 


.9791 


1 


.0000 


5.5 


-1.5999 


1 


.6196 


10. 


,2962 


-1.5999 


9.6382 


10.2962 


_ 


9.6382 


1 


.7701 


1 


.0000 


6 


-1.9833 


1 


.6795 


10. 


,4900 


-1.9833 


9.3595 


10.4900 


_ 


9.3595 


1 


.5494 


1 


.0000 


6.5 


-2.3711 


1 


.7283 


10. 


,7258 


-2.3711 


9.0030 


10.7258 


_ 


9.0030 


1 


.3453 


1 


.0000 


7 


-2.7626 


1 


.7677 


10. 


,9952 


-2.7626 


8.6015 


10.9952 


_ 


8.6015 


1 


.1678 


1 


.0000 


7.5 


-3.1571 


1 


.7996 


11. 


,2912 


-3.1571 


8.1974 


11.2912 




8.1974 


1 


.0174 


1 


.0000 


8 


-3.5544 


1 


.8255 


11. 


,6081 


-3.5544 


7.8254 


11.6081 




7.8254 





.8913 


1 


.0000 


8.5 


-3.9539 


1 


.8468 


11. 


,9415 


-3.9539 


7.5028 


11.9415 




7.5028 





.7858 


1 


.0000 


9 


-4.3555 


1 


.8645 


12. 


,2880 


-4.3555 


7.2323 


12.2880 




7.2323 





.6971 


1 


.0000 


9.5 


-4.7588 


1 


.8793 


12. 


,6450 


-4.7588 


7.0084 


12.6450 




7.0084 





.6221 


1 


.0000 


10 


-5.1637 


1 


.8918 


13. 


.0106 


-5.1637 


6.8235 


13.0106 




6.8235 





.5584 


1 


.0000 


10.5 


-5.5699 


1 


.9024 


13 


.3834 


-5.5699 


6.6702 


13.3834 




6.6702 





.5039 


1 


.0000 


11 


-5.9774 


1 


.9116 


13 


.7620 


-5.9774 


6.5421 


13.7620 




6.5421 





,4569 


1 


.0000 


11.5 


-6.3860 


1 


.9195 


14. 


,1455 


-6.3860 


6.4342 


14.1455 




6.4342 





.4162 


1 


.0000 


12 


-6.7955 


1 


.9264 


14. 


.5332 


-6.7955 


6.3426 


14.5332 




6.3426 





.3807 


1 


.0000 



Table 3.1: Selection of masses for the considered Lee model with parameters A = 
10i?o, M = 2E and E = 1. All numbers after the fourth digit are dropped. 
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name 


label 


determining equation 


description 




resonance pole mass 


(lIl)polG 


z - M Q - 


f g 2 Zu{z) = 


Real part of solution z. 




left /right pole mass 


•^(l/rjpole 


z-M 


+ g 2 T.{z) =0 


Real part of solution z. 










Nevertheless, z is purely 


real. 


Breit-Wigner resonance mass 


K (m)BW 


x - Mo + 


.9 2 ReS(i) = 


The solution x is purely 


real. 


Breit-Wigner left/right mass 


x (l/r)BW 


x - Mo + 


. 9 2 ReE(x) = 


The solution x is purely 


real. 


average mass 


^average 


* = h 

7T - 


f£ d P (u) 


Numerically evaluated. 




maxima of spectral function 


*^(l/r)max 


max p(E) 


Solution is x = E mSlX . 












Found numerically. 





Table 3.2: Description of masses for the considered Lee model. The general complex solu- 
tion of an equation is denoted as z = x + iy. 



Im z 



1.97 1.98 1.99 J»2.00 



Re c 



- 0.05 L 



Figure 3.6: Lee model pole (dark blue dots) for g £ [0.1,0.4] in the second Riemann sheet 
with Ag = 0.02, where the gray line marks the branch cut. 



The pole trajectory within the regarded interval [0.1,0.4] for the coupling is shown in 
Fig. |3.6| A summary of pole positions and additional information for different intervals 
can be found in Tab. |3.1[ To understand the adopted notation we have to define 
different masses of the considered Lee model (see also Tab. 3.2). The pole mass £( m ) pole 
of the main resonance pole is the real part of the complex solution found by solving 



Eq. ( |3.52 ), while a?(i) po ie an d £( r ) P oie are the real parts of the complex solutions in the 
first Riemann sheet, where we have replaced the self-energy Sn(z) in Eq. (3.52) by its 
equivalent on the first Riemann sheet. The corresponding Breit-Wigner masses are 
the solutions of the latter equation using just the real part of the self-energy on the first 
sheet (they are all real by construction) . These are simply the zeros in the denominator 
of the propagator on the real axis. The numerical maxima of the continuous part of the 
spectral function are denoted as xm mgx and £( r ) maX ) whereas x ave rage is the result of the 
numerically evaluated integral 



1 

^average 



A 

du up(u) . (3.54) 

E 



7T 

The last column of Tab. |3.1| simlpy gives the normalization of the spectral function 
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3.2.2 Couplings g 6 [0.4, 1.2] 



For higher values of the coupling the resonance pole descends deeper into the lower half 
plane of the second sheet, but around g ~ 0.9 there is a new pole arising in the first 
sheet left from the first branch point at z = Eq. The new pole seems to continuously 
enter the spectrum of found solutions]^] Although we could clearly fix g of this pole in 
the context of our above condition Abs G^ 1 (z po \ c ) < 1CT 10 , it is reasonable, and will 
be shown analytically, not to speak of the pole to suddenly 'pop up' on the real axis. 
Compare this situation to the behaviour of the resonance pole in the second sheet for 
vanishing couplings: these are obviously numerical problems we need to keep in mind. 



Im z 

o.o 
-o.i - 

-0.2 
-0.3 



-0.4 



1.75 1.80 1.85 1.90 1.95 _ 2.00 



Re z 



Figure 3.7: Lee model pole for g 6 [0.4, 1.2] in the second Riemann sheet with Ag = 0.05. 



Im z 



0.999 



Re z 



Figure 3.8: Left Lee model pole for g £ [0.4, 1.2] on the real axis of the first Riemann sheet 
with Ag = 0.05. 



This subtle point will be the main difference to the quantum field theoretical model in chapter 4. 
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3.2.3 Couplings g 6 [1.2,4.0] 

The presence of a new simple pole on the first sheet can also be noticed by the abruptly 



falling spectral function around g ~ 1.2, see Fig. 3.5 If we further raise g the pole 
passes the origin at g = 3.303, while the pole on the second sheet descends deeper into 
the lower half plane. Here, the real part decreases more slowly and at g = 3.069 it 
starts increasing. Additionally, there is another pole emerging around g ~ 2.7 on the 



Im z 



Re z 



( 



Figure 3.9: Lee model pole for g £ [1.2,4.0] in the second Riemann sheet with Ag — 0.1. 



Im z 







0.4 












0.2 








- 1.0 


-0.5 






0.5 


1.0 






-0.2 












-0.4 









Re c 



Figure 3.10: Left Lee model pole for g e [1.2,4.0] on the real axis of the first Riemann 
sheet with Ag = 0.1. 



first sheet to the right of the second branch point at z = A. We shall call the former 
the left (1) one and the latter the right (r) one, respectively, as was already explained 
during the discussion of Tab. 3.1 It should be clear that the pole masses for the left 



and the right poles are identical with the corresponding Breit-Wigner masses, because 
the imaginary part of the self-energy is only non-zero for Eq < z < A. Therefore, the 
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Im z 
0.4 - 

0.2 - 



-0.2 - 



-0.4 - 



10.00 10.01 



10.02 10.03 



Re z 



Figure 3.11: Right Lee model pole for g g [1.2,4.0] on the real axis of the first Riemann 
sheet with Ag = 0.1. 



purely real solutions of the equation 



= x - M + g 2 ReT,(x) 



are also solutions of 



(3.55) 
(3.56) 



= x - M + g 2 Y,(x) , 

which is the zero of the full denominator of the propagator. 

When looking at Fig 3.10, it seems that the left pole reaches a constant g-velocity 
in the parameter space spanned by Rez = x and the coupling g. Although the last 
equation cannot be solved for x analytically, it can be solved for g: 



g{x) 



4tt(M - x) 



In 



A— x 
E -x 



where we only take the positive branch from now on. The limit x 
derivative with respect to x in fact turns out to be a constant: 



(3.57) 
-oo of the first 



lim ^ 

x^-co dx 



\fn lim 



+ 



(E -A)(M -x) 



(*-£ )(*-A)ln 2 |^ 



(E - x)(x - A)ln( 



Mq-x 



111 



E -x 



A-x 
E —x 



x-Eq)(x- A) In ^ Eq _ x 

I / in 



Mp-x 



In 



E -x 



A -E n 



A-Ec 



4-7T 

A -En 



(3.58) 



45 




Figure 3.12: Analytic function g(x) (dark blue), numerical inverse x(g) (red) and asymp- 



totic lines from Eq. (3.621 (green, dashed) 



We obtain the same expression in the limit x — > +00 but with opposite sign: 



m+oo 



dg(x) 



X-S-+00 dx 



= lim 



47T 



A -En 



(3.59) 



The left and right pole, respectively, can then be taken to be located on a line g{x) = 
m± QO x + 6±oo with the above coefficients TTi-j-oo . The interception points (0|6-j-oo) are 
calculated again by taking the limit x — > ±00: 



lim g{x) 

x— >±oo 



m± 00 x + 6-too , 
=F 



1 



A-E, 



(A + E + 2M ) 



(3.60) 
(3.61) 



The resulting asymptotic line for the left pole gives deviations smaller than 5% in the 
coupling with x < —5.5 for our choice of parameters. For the right pole this is the case 
if x > 15.5. By inverting the linear equation (3.60) we finally get an expression for 
the function x(g) in the limit of huge couplings (here, deviations less than 5% can be 
obtained for g > 13 and g > 11.5, respectively): 



x(g) 



9 ~ b± c 



± 



A - E A + E + 2M 
9 + * 



4-7T 



(3.62) 



Both of these lines intercept at the same point on the x-axis. The limit of constant 



g-velocity is also observed for the right pole. In Fig. 3.12 we give the numerical inverse 



of Eq. (3.57) together with the asymptotic lines from Eq. (3.62). The most important 



result of that plot is the fact that, although we could not find poles for some intervals 
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of the coupling (no right pole for g < 2.7 and neither left nor right poles for g < 1.2), 
they really do exist near the branch points at z = Eq and z = A. A numerical analysis 
cannot fully reveal the full pole structure since one needs a high starting point precision 
beyond 10 -10 . Nevertheless, we can state now that there are poles on the real axis as 
long as g is non-zero. This is because the following limits indeed exist, 



lim g(E - rj) 



lim 

»7->-0+ 



lim 

r/->-0+ 



4tt(M - (Eq - n)) 



In 



A-(E - V ) 
Eo-XEq-v) 



4vr(M - E + rj) 



V 



(3.63) 



lim g(A + rj) 

r/->0 + 



lim 

r/->-0+ 



lim 

r/->-0+ 



4tt(Mo - (A + rj)) 







In 



A-(A+r;) 
Bo-(A-Hj) 



4vr(M - A - 77) 



In 



-v 

E -A-ri 



(3.64) 



and the root functions are both positive for z < Eq and z > A. It is satisfying to see the 
few available analytic expressions not only supporting the numerical analysis but indeed 
completing it. We may also show that the spectral function p(E) has delta peaks for the 



left (x = E < E ) and right (x 



p(E) 



E > A) poles in the case of non- vanishing couplings: 
e + g 2 lmT,(E + ie) 

2 



M + g 2 ReJ:(E + ie)) + (e + ff 2 Im£(£ 
e 



re 



))" 



— M + 5 2 Re E(,E + ie)) + e 2 
vr5(E-M + 5 2 ReE(S)) 
7r5(E-M + 5 2 S( J E)) . 



(3.65) 



3.2.4 Couplings # G [4.0, 6.0] 

The pole in the second Riemann sheet falls deeper into the lower half plane. The real 
part decreases but the pole seems to approach the constant line Hez = 2. This can be 
investigated in a similar way to the previous cases, yet now the full complex denominator 



of the propagator needs to be studied. After solving Eq. (3.52) for the positive coupling 
g we obtain the master solution 



g(z = x + iy) 



4tt(M - (x + iy)) 



In (A-(x + iyj) - In [Eq - (x + iyj) + 2ni 



(3.66) 
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from which one immediately observes the multi- valued character of the propgator due to 
an imaginary shift of 2ni, coming from the logarithmic structure. The numerical results 
already indicated that the imaginary part of the resonance pole goes to minus infinity, 
while the real part approaches the value Rez = 2. If we assume the latter to be true, 
then the coupling g must grow beyond all bounds when the imaginary part y goes to 
minus infinity. And this is in fact true since 



g(z = M + iy) 



4tt(M - (M + iy)) 



In (A - (M + iy)) - In (E - (M + iy)) + 2m 



lim g(z = Mo + iy) 

y— >— oo 



-4iriy 



ln(A - M - iy) - ln(E - Mo - iy) + 2iri 



(3.67) 
(3.68) 



The two poles in the first sheet just keep 'walking' on the real axis, the left one going 
to minus and the right one to plus infinity. We have just proved this behaviour in the 
previous subsection (the limit of constant ^-velocity in parameter space also holds for 
the right pole, see lower figure). 



0.2 



-0.2 



-0.4 - 



10.0 10.1 10.2 10.3 10.4 10.5 10.6 



Re z 



Figure 3.13: Right Lee model pole for g £ [4.0,6.0] on the real axis of the first Riemann 
sheet with Ag = 0.1. 
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3.2.5 Couplings g 6 [6.0, 12.0] 



The real part of the resonance pole approaches the constant line Rez = 2, see Fig. 3.14 



while the two poles in the first sheet further walk away on the real axis. When looking 



back at Tab. 3.1 we see that the values for the Breit-Wigner mass X( m )BW are very 
different for small and high g, seperated by a gap, while for g — > oo it seems to equal the 
value of X( r ) max . Numerically, one obtains for our choice of parameters £( m )BW ~~ * 5-5 in 
the limit of very high couplings g 3> 10. A natural suggestion therefore would be 



A + £ q 



lim_ x (r)max = lim x (m)B w = — in? — , (3-69) 



since this involves the only free parameters of our model and gives the right numerical 
value. Now, by considering the spectral function p{E) right above the branch cut, it is 



Im z 



Re z 



-20 



-40 - 



-60 - 



-80 - 




Figure 3.14: Lee model pole for g G [6.0, 12.0] in the second Riemann sheet with Ag = 0.2. 



possible to show explicitly that this is not the right limit. Remember that for Eq < z < A 
we cannot use the self-energy as calculated at the beginning of our investigation in Eq. 



(3.29) because the /c-integration would hit the pole on the real fc-axis. For that reason 



we need to perform the integration by using the Sokhotski-Plemelj theorem: 

A 1 f A 1 

dk = V dk + in 



E k-E-ie J Eo k-E 



O—77 1 rk—E [ 

dk - + dk - I + iir 



E ~E k J 0+7] k 

= ln(-T]) - ln(E - E) + ln(A - E) - In 7? + in 

= In T) + in — \n(E — Eq) — iir + ln(A — E) — In r/ + m 

where we have used the well-known definition of the logarithm for negative real numbers. 



This result can also be obtained by applying the latter to the self-energy (3.29) itself. 
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The spectral function now reads: 
p(E) -- 



e + g 2 Im£(£ + ie) 



(E-M + g 2 Re E(E + ie)f + (e + g 2 Im £(£ + ie)f 



1 



1 



E-M + ^ln 



A-E 
E-E 



(3.71) 



+ 



16 



A propagator pole on the real axis would correspond to a zero in the denominator of the 
last expression. But for the regarded interval there is no such zero except for E = Mq 
and g = 0, since the general solution for positive couplings, 



g( X = E) = 2v% 



(Mo - E) In - iry/-(M -E)* 



(3.72) 



is only real (with value g = 0) for i£ = Mq (of course it makes no sense to deal with the 
self-energy at all for vanishing couplings). The value x max is consequently the minimum 



of the denominator in Eq. (3.71): 



_d_ 
dE 



g 



l - 2 ( E -M + 9 -\n 



A-E 

E — En 



+ 



(E -A) + p(E - E)(E - A)] [g 2 In (j=g ) + 4vr(£ - M 



8it 2 (E - E )(A- E) 



(3.73) 



The denominator 8ir 2 (E — Eq)(A — E) vanishes for E = Eq and E = A which are 
independent of g, so there is no problem in discarding this term. On the other hand, 
the first part of the numerator has two solutions 



E± 



E + A ± J 9 - (E - A) + (E - Ay 



(3.74) 



where only the negative solution belongs to a local minimum (and to a local maximum 
of the spectral function) for g > 0. We cannot prove this in general, for our equations 
are transcendental, nevertheless the mentioned aspect becomes clear when plotting the 



second derivative of the denominator in the spectral function (3.71) and inserting the 



negative minumum solution from Eq. (3.74): the graph (now with dependence of g) lies 



in the first quadrant and there is no change until both solutions become complex for a 
numerically computed value of g > 5.317. This is in full agreement with the analytic 



condition for complex solutions coming from the root in Eq. (3.74), 



g_ 

7T 



(Eq - A) + (Eq - A) 2 = , 



9 



^-ir(E - A) « 5.317 , 



(3.75) 
(3.76) 
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and can be also seen in Tab. 3.1 where the maximum xn) mgx vanishes for g > 5.0. The 



other maximum x 



(r)max 



(l)max 

comes from the logarithmic term in Eq. (3.73). Its zeros are 



obviously the Breit-Wigner masses »( m )BW> too: 

g(x = E) 



4tt(M - E) 



In 



A-E 
E-E 



(3.77) 



The relevant behaviour for us to study now is the limit of huge couplings. Because 
the logarithmic term in the first derivative from Eq. (3.73) also gives the remaining 
minimum for g > 5.317, we only need to consider the limit of the last above expression. 
Our first suggestion 



lim x 

g— >+oo 



(r)max 



lim x 

g— >+oo 



(m)BW 



A + E 

M 



yields the correct numerical value of 5.5 (this is true if we just consider the right-sided 
limit): 



lim 

?7->-0+ 



4vr(2- (5.5 + t/)) 



lim 



lim 

J7->0+ 

oo . 



4tt(2 - 5.5 - rj) 



-14vr — 47T77 



4.5— t? 
4.5+J? 



(3.78) 



Nevertheless, this does not hold in general. In fact, our first suggestion has to be modified 
in the following way: 



lim x 

g— y+oo 



(r)max 



lim x 

g— >+oo 



A + E t 



(m)BW 







(3.79) 



which reveals the limit to be independent of the mass Mq. A real right-sided limit does 
not exist for arbitrary choices of parameters. We find in general: 



lim g(x = (A + E )/2 + rj) 

r;->-0+ 



lim 

r;->-0+ 



47T (M 



' A+Ep 
, 2 



+ r/)) 



In 



lim 



-2vr(A + E - 2M + 2rj) 

ln ( A-E -2r! \ 
111 \A-E +2 V J 



ooJ- sgn ((Eq - A) (A + E - 2M )) . (3.80) 
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3.3 Spectral function and its normalization 



After having investigated the pole trajectories and the behaviour of the spectral function 



p(E), let us now focus on the normalization of the latter. As was shown in Fig. 3.5 the 
spectral function is not normalized when using the numerical result of the integral 

1 f°° 

doj p{oj) , (3.81) 



7T 

simply because the numerical calculation does not consider the two poles arising on 
the real axis in the first Riemann sheet Q It is possible to take them into account by 
splitting the Schrodinger propagator Gs(E) into the contributions of two single-particle 
propagators and the remaining continuum part above the branch cut, similar as it was 
done in the first chapter: 

J-oo E-u + ie 



G S (E) 



Z 



(l)polc 



E 



"(l)pole 



+ ie 



+ 



^(r)pole 1 
E - Z( r )pole + ^ 7T 



do; 



p(u>) 



E 



E — uj + ie 



(3.82) 



where ^(i)p ie an d ^(r)poie are some renormalization constants. We calculate them by 
expanding the inverse propagator in a Taylor series around the two poles on the real 



axis at x 



"(l/r)pole 



in first order: 



Mo + 5 2 E(x) 



M + g 2 Z(x))\ x=x 



(l/r)polc 



d 

+ — I x 
ax 



M + 5 2 E(x) 



(x - X (1 / r ) pole ) 



d 



1 + 5 Z — ReS (x 

ax 



x a '(l/r)polc 

• (x - ^(l/r)pole) 



c (l/r)pol C 



1 + 



A-£ n 



4tt (A-x)(Eq 



(x - £(l/ r )pole) 



c (l/r)polc 



(3.83) 



Since the two simple poles and the branch cut are separated on the real axis, the renor- 
malization constants for each single-particle pole can be extracted from their residues: 

Z, 

(x - X(i/ r ) po ie) • G s {x) = ^ lim (x - x ( i /r)polc ) 



lim 

x-s>x (1( / r)polc 



y (l/r)pole 



X ~ X (l/r)pole + «C 



Res(G s (^), z=iE(i/ r ) P oic) 
| other pole and branch cut j 



Z 



(l/r)pole • 



(3.84) 



7 The normalization integral from Eq. (3.811, if calculated numerically, is computed only across the 



branch cut. Any extension to arbitrary bounds of integration beyond the cut does not change the 
result except for irrelevant fluctuations from the poles on the real axis. Hence, we write the bounds 
as plus and minus infinity. 
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If we insert our expansion of the inverse propagator from Eq. (3.83) into the left side, 
we finally arrive at: 



^(l/r)pole 



1 + 



9 



4-7T (A - aJ(l/ r ) p ole)(-So ~ £(l/r)pole 

Hence, the normalization condition can be written as 

,A 

1 = ^(l)pole + ^(r)pole + / dw p(uj) , 

Je 



(3.85) 



(3.86) 



which makes clear why the integral along the branch cut gives nearly the correct num- 
ber for small couplings (where the renormalization constants are very small). If we 



0.6 



Figure 3.15: Numerical verification that ^ J_ do; p(ui) — 1 (red) by including the single 
poles of the first Riemann sheet, and numerical contribution of the continu- 
ous part ^ J Eo duj p{uj) (dark blue). The necessity of the delta distribution 
functions is evident. 



perform our calculation including the renormalization constants, the spectral function 
is then normalized over the full range of g. In particular, sufficiently large numerical 
contributions from Z^ po i e are obtained for g > 1.2 and from Z( r ) po i e for g > 2.7^] The 
plot in Fig. |3.15| shows the completely fulfilled normalization condition for our choice of 
parameters (the points used are those from Tab. 3.1). For the sake of completeness, we 



provide a compilation of selected plots of the continuous part of the spectral function in 
Appendix E| from where one can review the behaviour described during this last section. 
An example for high couplings - including the two poles in the first Riemann sheet - 
can be found in Fig. |3.16| 



8 This is also true if we perform the average mass integral in Eq. (3.541 including the renormalization 
constants. Then, the value of :r av0 ragc 



is Mo for the whole range of the coupling. 
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Figure 3.16: Spectral function p{E) of the considered Lee model for g = 6.0. 



3.4 Time evolution and non-exponential decay properties 



We are also interested in the survival amplitude a(t) as denoted in Eq. ( 3.14| ) to make 
some statements about the decay of the unstable particle S in view of the analytic 
structure of its propagator. We have found an explicit result for this propagator in the 
previous section. The survival amplitude is obtained by using a Hankel contour-like 
path of integration in the complex E'-plane. We first extend the integration onto the 
whole complex plane and then close it counterclockwise in the lower half plane, while 



splitting the propagator in its single parts, see Fig. 3.17 The only contributions to the 
amplitude come from the residues of the two simple poles and from the branch cut. In 
a more pedagogical way, one may also express the spectral function as the difference 
between the propagator across the cut (in the limit e — > + ) and integrate along the real 
axis in a first step from — oo to +oo, in a second step from +oo to — oo. Either way, the 
survival amplitude becomes 



a(t) 



2tt 

i 

2^ 



J +0 ° dE e- iEt G s (E) = ±- jfdz e- izt G s (z 



Z(Y)pole <*> dz 



-izt 



d z ~ Z(l)pole + le 



(r)polc 



-izt 



dz 

C r z ~ £(r)pole + *e 



A 



-izt 



H / duo p(uj) <b dz — 

7T Je Jc z-u + ie 



i 

2^ 



27r ^(i) P oie Res 



-izt 



Z = X 



27riZ( r ) pole Res 



z — X 



z - x (1 ) pole + te 

g—izt 



(l)polo 



(r)pole ' ^ 



2i 



/ dw p( 

JEo 



oj) Res 



-izt 



z — u + ie 



z = UJ 



1 



^(i) P oiee" i:C(1)polct + ^(r) P oiee- M 'M^ + - ^ du e^* p{u>) 



(3.87) 
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Figure 3.17: Contour paths of integration C\, C r , and C (red, dashed) enclosing the two 
simple poles and the finite branch cut when blowing up the full contour C 
(black, dashed). 



and thus the survival probability p(t) = \a(t)\ 2 consists of different terms: First of all 
we get the squared renormalization constants with g-dependence belonging to stable 
particles, because the imaginary part of the two poles on the first Riemann sheet vanish. 
If we had instead just one simple pole with finite negative imaginary part at z po \ c = 
M — iT/2 in the lower half plane, the well-known exponential decay law would have 
been reproduced: 

o(t) 



Note that in general all resonance poles lie on the second Riemann sheet. Pure expo- 
nential decays therefore cannot be realized exactly. This fact is due to the mathematical 
structure of the propagator: the self-energy Ti(E) is a multi-valued complex function, 
creating a branch cut on the real axis whenever the spectral function gives continuous 
support in the integral representation of the propagator (as was already explained in 
the first chapter). All poles on the cut slip through it into the second sheet - as a con- 
sequence, the survival amplitude will have a cut contribution (as well as contributions 
from other poles in the first sheet) which changes the familiar decay law. It is well-known 
that the exponential behaviour can be preserved by approximating the self-energy at its 
value on the real axis [55] . In this case the resonance pole is assumed to be situated 
near and below the real axis, i.e., the imaginary part y po ie = —T/2 to be much smaller 
than the real part x po ie = M. The self-energy in the second Riemann sheet around the 
pole can then be taken at z = Mo — 5M + ie = Mbw + above the branch cut in the 
first sheet. This is obviously the Breit-Wigner approximation! Going over to general 



i 
i 

2^ 



oo —iEt 

dE 

- oo -^pole 
p — izt 

dz — 

C z ~ z po\e 



-rt 



(3.88) 
(3.89) 
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expressions, including a normalized form factor function / 2 (k)/(27r) = / 2 (k), this means 
En(* = M-iT/2) » £(z = M BW -He) 

/ 2 (k) 



7> / dk 



.-m bw +ot/(Mbw) 

and thus one can write for the Schrodinger propagator around the pole 



(3.90) 



Gs(z) 



1 



z-M + c/ 2 £(M BW + tc) 



M + 5 2 Re E(M BW + *e) + » {y + 9 2 Im £(M BW + ie)) 



.7; 



<5M 



5 2 ReE(M BW + ie) , 
2 5 2 ImE(M BW + ie) 



2^ 2 / 2 (M BW ) , 



(3.91) 
(3.92) 
(3.93) 



where the last relation turns out to be Fermi 's golden rule for the decay width (sometimes 
M B w is replaced by Mo). If we make the unphysical assumption that the energy spec- 
trum of the full Hamiltonian has no lower bound (and no upper bound), a non-relativistic 
Breit-Wigner distribution is immediately obtained for the corresponding spectral func- 
tion |54| . On the other hand, since the spectral function displays a pole in the second 
sheet right below the real axis as a relatively sharp peak, the branch cut contribution 
can be replaced by a simple pole in the lower half plane of the first sheet, yielding the 
exponential decay law (3.89) for the survival probability p(t). 

Besides the branch cut, we need to take the real poles on the first sheet with real 
renormalization constants into account. The complete survival proability reads: 

^{ x (r)po\c~ x (\)pole)t _j_ g — ^( x '(r)polc — x (l)polc)^ 



P(t) 



z 



+ z, 



(r)pole 



(l)pole •^(r)pole 
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+ -f^mpoi e e" (1)polci + Z Mpole e ix (^ 
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7T 
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7T 
1 

V 2 



y (l)polc 



e -W!(l)polet _|_ 2 



(r)polc 



^^'(r)pole^ 



du> e 

A 



^p(u) 
doj e _iwt p(w) 



E 



du: e- iut p{uj 



Eo 



dco e'^piuj) 



Eo 



After introducing 



Eo 



duj e lu) p(u>) , w r i = x (r ) pole 



X 



(l)polc 



(3.94) 



(3.95) 



this becomes the more compact expression 
p(t) - ^ 1 ^ 2 



^(l)pole + ^(r)pole + 2 ^(l)pole^(r)pole COs(ui T \t) 



+ -Re«S( Z, 

7T 

2 



7T 



(l)pole ^(^(^pole*) + ^(r)pole COs(X( r ) 

Im5^Z (1)pole sin(x (1)polc t) + Z( r)pole sin(x (l .) po i c t)^ + ^|5| 2 . (3.96) 
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In the following Fig. 3.18 we show the full survival amplitude p(t) from Eq. (3.96) (red), 
the branch cut contribution S (dark blue) and the single pole approximation as denoted 
in Eq. ( 3.89| ) (green, dashed) for some values of the coupling g. We observe: 



• Small couplings g < 1: 

The renormalization constants and all terms containing at least one of them are 
negligible since the self-energy has sharp cusps at the branch points, making the 



expression (3.85) zero. The survival probability consists mainly of the branch cut 
contribution S, which can be approximated as an exponential function. In the 
limit of very small times we realize p(t) to flatten with vanishing first derivative 
- a feature observable in the full range of the coupling and closely related to the 



famous quantum Zeno effect 61-65 



• Sizable contribution from left pole, g = 1.4: 

The pure exponential law does not hold any more, and the branch cut contribution 
also does not give the full survival probability. This is clear because the left pole on 
the first Riemann sheet brings a relevant renormalization constant into play and 
hence the probability (though very small) stays nearly constant for large times as 



a direct manifestation of the first term in Eq. (3.96). Oscillations given by the 



cosinus of the fourth and the sinus of the fifth term are very weak and finally vanish 
due to the real part of the branch cut contribution tending to zero. 



• Sizable contributions from both poles, g > 2.8: 

For large times the survival probability p(t) oscillates with the frequency uj t \ be- 
tween the amplitude values ^ ( 2 1)pole +^ ( 2 r)po i e +2Z (1)pole Z (r)pole and ^ )pole +Z ( 2 r)pole - 
2Z(i) po i e .Z( r )p i e . The intermediate oscillating behaviour with a much lower fre- 
quency is due to the interplay of the real and imaginary parts of the branch cut 
contribution in the last three terms. This transient overall oscillation is not im- 
portant for the evolution of the system and tends to zero, because the spectral 
function p(u) goes like 0(g~ 2 ) in the highest order for huge couplings. 
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3.5 Concluding remarks of the chapter 



The Lee model presented here provides us with a variety of interesting deviations from 
the general transition properties of unstable quantum systems as induced by the purely 
exponential decay law: 

1. Although the Breit-Wigner parameterization reproduces very well the mass (and 
decay width) of the resonance pole in the regime of small couplings, there exists 



no such mass (and width) for intermediate values of g, compare Tab. 3.1 Above 



this regime, the mass (and width) is completely detached from the pole values. 

2. The spectral function describes the pole motion adequately through its maximum 
x (r)max = ^(m)BW an d full width at half maximum as long as the two newly arising 
poles can be neglected. Otherwise, the maximum does not correlate with the real 
part of the resonance pole in the second Riemann sheet. 

3. The left and right poles have their origin in the form of the self-energy. While the 
right pole disappears after removing the right branch point by sending A — > oo (and 
performing a subtraction), the other pole would remain. We interpret the latter as 
a bound state emerging from the strong coupling to the continuum of states. One 
should stress that the whole model is academic, but serves to demonstrate how to 
perform the analytic continuation into the second Riemann sheet. At least smooth 
cutoff functions have to be used since the right threshold (and consequently the 
right simple pole in the first sheet) is unphysical. 

4. Deviations from the decay law can be obtained for small and large times due to the 



branch cut contribution. However, we do not find the theoretical power law 54 
for very large times. If present in that regime, it is covered by the dominant 
oscillations coming from the left and right pole. 

The two poles in the first Riemann sheet are dynamically generated by the in- 
teraction term appearing in the full Hamiltonian. There is no evidence for other 
resonance poles created in the second sheet. 
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4 Quantum field theory of resonances 



4.1 Construction of the £</></>- model 
4.1.1 Lagrangian 

In this chapter we study a quantum field theoretical model with two scalar fields, S and 
(f>, in d = 3 + 1 dimensions. Such a model can be applied to many problems in the realm 
of hadron physics, most notably it is possible to investigate the resonance pole of the 
cr-meson. The same model was already studied by Veltman [66] as well as by Giacosa 
and Pagliara [6~7]|68| regarding its spectral function(s): it contains an interaction term 
for a one-channel decay process S — > <j><j> (see Fig. |4.1[ ), described by a Lagrangian of the 
form 

C = l 2 {d^S){d»S) + \{d^)(d»4>) - \M*S 2 - ^m^+gScp 2 . (4.1) 

Here, the (pseudo) scalar </>-fields represent pions. The equation of motion for the a- 
meson alias the field S, 

(□ + M 2 )S = gtf , (4.2) 



is of course the solution of the Euler-Lagrange equation 

dC d dC 



dS dx» d(d^S) 



. (4.3) 



The corresponding Feynman propagator in position space for vanishing coupling g 
is defined as in the first chapter, 

Af%x-y) = -i(0\T{S(x)S(y)}\0) 

r „—ip-(x—y) 



-ip-(x-y) - /aa\ 

(2tt) 4 P 2 _ M 2 + ie ' V'*) 

where T is the time-ordering operator. Hence, the propagator in momentum space 
becomes simply: 

A * C V) = 2 L - ■ • (4-5) 

One should mention that the theory described by the above equations is known to be 
super-renormalizable. In fact, the superficial divergence D in terms of the number of 
vertices n and the number of external lines is proportianal to — n, so the number of 
divergent diagrams is finite. 
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Figure 4.1: Schematic decay process S —> 4>4>. 



4.1.2 Tree-level decay width 

The decay width of the resonance S can be calculated by evaluating the invariant S- 
matrix element of the process S — > (jxfi in first order in the coupling. In the case of our 
simple interaction term C m t = gScf) 2 the S'-matrix can be immediately written down as 

S = t + ig J d 4 x T{:S{x)cp 2 {x):} 

= 1 + , (4.6) 

where the dots mark the normal ordering prescription. The crossed-out terms in the 
resulting transition matrix element, 

(final|5|initial) = (piP 2 |£ (1) |p) 

= (P1P2I i9 J d 4 x:gj^^ 

+ gT¥^2:|p), (4.7) 

give no contribution because the creation and annihilation operators combine in such a 
way that their scalar product vanishes. The superscript at the S- and (f)- fields denotes 
those parts of the field that contain a creation or annihilation operator, e.g. 



</>(+) = <^+)( 



x 



1 bU ip - x . (4.8) 



(2vr)3 ^2E V 



Writing out the Fourier expansion of the fields in full detail and using general commu- 
tation relations for the bosonic creation and annihilation operators, the remaining steps 
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are straightforward: 



■ [a* f a3 ' f a3 ' f a3 ' ^2E P 2E P 2E P 
ig I d x / d pi I a p 2 / d p — 

J J J J J2E„,2E„,2Ew(2Tr) 9 



-Pi— p 2 — P '(2vr) E 
x e i W+^^')- a; {0|6 P2 6 Pl 6j ) ,6 t p ,a p ,at ) |0) 

v / d 4 x / d 3 P ; / dV 2 / dV ^ 2E ^ 2E ^ 2E » 

J J J J . 2E„> 2E„< 2E„> 



2Ep>2Epi2Ep. 

x e iW+^- P ')-| (5 (3) (p / _ p)5 (3) (pi _ p' 2)(5 (3) ( p 2 _ p ; } 
+ <5( 3 )(p'-p)^)( Pl -pi)5( 3 )(p 2 -p' 2 )} 

= 2ig j d 4 x e<(w+M-p)-* 

= 2ir7(27r)V 4 )(pi+p 2 -p) 

= -iM (2vr) 4 5 (4) - pi - p 2 ) , (4.9) 
=> -iM = 2i# . (4.10) 

The decay width is now obtained by performing the phase space integral over the invari- 
ant amplitude —iM. The infinitesimal expressions read: 

'-E»)(n^). <^> 

where 5 is a symmetry factor and n is the number of particles created. For a general 
discussion of these formulas see for example Ref. [19] . In our case the symmetry factor is 
just one half - since the directions of the outgoing momenta are determined by conser- 
vation laws, only a half sphere in position space needs to be taken into account. There 
is no angular dependence as a consequence of Lorentz invariance: spinless particles at 
rest have no preferred direction in which they decay into two other particles. We then 
quickly arrive at the phase space integral for the tree-level decay width, 

r tree (V^ = M) = I d 3 Pl / d 3 P2 oEl x — s^( P - Pl - P2 ) 



(2n) 2 M J "J ^ 2E Pl 2E P2 

=5( 3 )( Pl +p 2 )5(M-B P1 -Bp 2 ) 



9 2 I A3, 



^g 2 r°° du 

(27r) 2 M J U 16(u2 + m2)v /Mi_ m 2 



4vrc? 2 f°° , Mu 2 nrj^r t ^ 

1 du 5(u - yM 2 /4 - m 2 ) , 



(4.13) 
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where u = \pi\- A real- valued expression requires a threshold value 2m, so we add a 
Heaviside step function to the preliminary result 



r -(M) = ^^e(M-2m) 



(4.14) 



Here, we have redefined the coupling as gs<f><f> = V%9 and introduced the magnitude of the 
three-momentum |p 1 | = |p 2 | = Ps4>4> = \/M 2 /4 — m 2 of the two outgoing (^-particles. 



The latter is derived from ordinary kinematics, see Appendix C for more details. 



4.1.3 Self-energy and interacting propagator 



In a free field theory, as was explained in the first chapter, the propagator (4.4) has the 



simple interpretation of being the amplitude for a stable particle S to propagate from 
y to x. This is different if we turn on the interaction: the field S cannot be considered 
anymore as an asymptotic state of the Lagrangian since loop contributions given by 
strongly coupled hadronic intermediate states dominate its dynamics [3l] . They appear 




Figure 4.2: Self-energy £(p 2 ) as a mesonic loop. 



in the full interacting propagator of the field S, so that the free propagator is modified: 



A s (p 2 ) 



1 



(4.15) 



The emerging self-energy function T,(p 2 ) in the denominator can be represented by a 
one-particle-irreducible (1PI) diagram in momentum space with two (/>-particles inside 
a mesonic loop where the incoming and outgoing momentum is denoted by p, see Fig. 



4.2. The full propagator (4.15) is obtained by a resummation of the one-loop diagrams, 
as shown on the next page, and separating its inverse by multiplying with (Ae ee (p 2 )) 1 
from the left and (A^(p 2 )) 1 from the right. From general Feynman rule^ for scalar 
fields we get the following expression for the self-energy: 



5V) 



d 4 q 



/!(<?) 



( 2? 4 (I + qf - m? + ie (f - qf - m 2 + ie 



(4.16) 



1 We adopt the Feynman rules used by Peskin and Schroeder 13 , yet the minus si gn in front of the 
self-energy contribution II(p 2 ) appearing in the denominator of the propagator (1.321 is absorbed 
into £(p ). It is therefore II(p 2 ) = -g|^ E(p 2 ). 
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Figure 4.3: Modification of the propagator by resummation of 1PI diagrams. The free 
propagator is represented by a solid line with an arrow and the interacting one 
by a double line with an arrow. 



The regularization function f\(q), which depends on a UV cutoff scale A, was introduced 
to make otherwise logarithmic divergent integral finite. Since we deal with an effective 
Lagrangian in the low-energy regime to study light mesons, the mass scale of our model 
is determined: it is reasonable to set A between 1 and 2 GeV, yet it should be possible to 
perform a full detailed study and eliminate any dependence on the cutoff. The integral is 
then evaluated by assuming the regularization function to depend only on the magnitude 
of the three-momentum q. Actually, this violates Lorentz symmetrjj^J but we accept this 
drawback in order to arrive at simple analytic expressions for Y,(p 2 ). First we split the 
g-integration into its spatial and time part, continue the latter into the complex and 
perform the integration in the rest frame of the particle S with mass y/p* = M by 
making use of the residue theorem: 

d 3 g 
(2^ 
dq° 



m) = -i 1 7^x3 /i(lql) 



27T (f + q ) 2 - |q| 2 - m 2 + ie (f - q ) 2 - |q| 2 - m 2 + ie 



(1) (2) 



(4.17) 



from this we have simple poles 

(1) q° 1/2 = -M-± x /\ (l \2 + m2-i e = -M.± y /\ fl \2 + m 2^ie , (4.18) 

(2) q° 1/2 = — ± v / |q| 2 +m 2 -ie = — ± \Z\q\ 2 + m 2 ie . (4.19) 

Because the final result does not depend on how we close the contour, only two poles 
need to be considered, namely those with negative root, situated above the real g°-axis. 
They give a path of integration circling counterclockwise in the upper half plane. The 
residues 

Res(£(M), s = o2 m ) = : 1 : , (4.20) 

V V ' 2{1) ' 2Mv/|q| 2 + m 2 (M + 2 v / |q| 2 + m 2 ) V ' 

Res (E(M), z = qi (2 \) = ; 1 ; , (4.21) 

V V ; 1(2); 2M v / |q| 2 + m 2 (Af - 2V|q| 2 + m 2 ) V ; 



2 See also Ref. 691 for more details on this aspect. 
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can be combined and, after performing the polar and azimuthal integration of the spatial 
part in spherical coordinates, the self-energy reads: 



S(M) 



2vr 2 



du 



« 2 /X(«) 



o Vu 2 + m 2 (4(u 2 + m 2 ) - M 2 - ie) 



(4.22) 



where we have shifted the integration path slightly into the lower half plane of the 
complex u-plane (u = | q| ) . 



Our preliminary result for the tree-level decay width in Eq. (4.14) is surely correct, 



but it is crucial to realize that the cutoff does not exist in our model on the Lagrangian 
level. It can be implemented by taking a non-local interaction term into account, 



£mt = gS(x)4> 2 (x) 



gS(x) / d 4 y 4>{x + y/2)(j>(x - y/2)$(y) 



(4.23) 



where the regularization function appears as the Fourier transform of <£(?/) in the loop 
integral and 'heals' the divergence in the case of Ja{q) = /a(|q|) [67] • This changes also 
the tree-level result for the width such that 



r tree (M) -> r tree (M) • fl(ps4>4> = VM 2 /2 - m 2 ) , 
which is demonstrated explicitly in Appendix D 



(4.24) 



4.2 Sharp cutoff 

4.2.1 Analytic structure of the propagator 

In order to proceed we need to determine /a(|<i|)- Two different choices for the regular- 
ization functions will be studied in this work and for the sake of simplicity we start with 
a sharp cutoff, 

/ A (|q|) = 9(A 2 - |q| 2 ) . (4.25) 



Inserting this function into Eq. (4.22), we are left with a non-trivial integral, which 
nevertheless has an analytic solution^ (that is found by Mathematica if we set e = 0): 



, , v^m 2 - M 2 ( KM \ 1 / m 

= ^977 arctan /,„ . o rr o , ,o " ^ ln 



8vr 2 M Vv / A 2 + m 2 v / 4m 2 -M 2 / 8ir 2 \A + VA 2 + m 2 , 

(4.26) 

Since this function has only a physical meaning for ^/p 2 = M > 0, we will not be con- 
cerned in what is happening for negative values of the mass. This is really an important 



3 One should never forget that in order to look for poles in the second Riemann sheet in the same way 
as was done for the non-relativistic Lee model in the previous chapter, we really need an analytic 
expression for the self-energy. Otherwise, the continuation into the complex plane has to be performed 
by a (finite) series expansion, while the continuation into the second sheet would be fully non- 
trivial. This way can be avoided by computing the loop integral solely numerically (while numerical 
distortions in the vicinity of the branch cut have to be managed). 
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point to mention: in the literature, people usually work in the complex p 2 -plane and 
indeed, after all calculations in this work were done, we found that it is much easier to 
do so as well. 

A first hint of the multi-valued character of the above self-energy function is given 
by the root. Clearly, there should be a branch cut in the complex M-plane (s/p^ = 
M —> z = x + iy) starting from z = 2m and heading to plus infinity on the real axis. 
Furthermore, since the complex inverse trigonometric function arctan has in general an 
ordinary series representation (with convergence radius r = 1, which already indicates 
the presence of a singularity, e.g. a branch cut), there is no problem at all to make the 
entire self-energy T,(M) a complex multi- valued function. However, we immediately face 
a strange property: plotting the self-energy wit h A = 1.5 GeV shows no kind of dis- 



continuous behaviour above 2i/A 2 + m 2 , see Fig. 4.4 This curiosity can be clarified by 
investigating the quite complicated inverse tangent inside the self-energy. The common 
arctan has two branch points at z = ±« (excluded from its domain) and conventionally 
the branch cut does not connect them both by a straight line on the imaginary axis 
through the origin, but rather by two straight lines. One part of the cut is along the 
positive imaginary axis above i, and one along the negative imaginary axis below —i, 
both heading to complex infinity: 

branch cut of z h-> arctan z = w : {(— ioo, —i], [i, ioo)} . (4.27) 

This cut structure determines the range of the arctan for real z = x as (— f,f) and 
usually the function is well-defined right from the upper and left from the lower part. 
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Looking for branch points in our case means 



(1) 
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(2) 



m 2 \/4m 2 
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m 2 \/4m 2 
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lim (2m + 77) 



(4.28) 
(4.29) 
(4.30) 
(4.31) 



so obviously the situation looks quite complicated. It is therefore much easier to find the 
discontinuity of S(z) on the real axis in exact the same manner as we did for the non- 
relativistic Lee model. There, we tried two ways and we may begin with the simpler one 
using the integral representation of the self-energy from Eq. (4.22). The discontinuity 
turns out to be: 
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u 



\/u 2 + m 2 (4(u 2 + m 2 ) - M 2 - ie) 



M 2 + ie) 
2ie 



A 



du 



Vu 2 + m 2 (4(n 2 + m 2 ) — M 2 ) + 
5{A{u 2 + m 2 ) 
1 



u 



\Ju 2 + m 2 
A „,2 

du 



u 







V« 2 + m 2 g /M£ _ m2 



M^j 



Af 2 
4 



4vrM 



2m < M < 2y / A 2 T 



(4.32) 



where the restrictions are predefined by the cutoff and the necessity of purely real values 
in the argument of the delta distribution function. It is also possible to argue that in 
order to avoid missing a real pole in the complex u-plane during the integration from 
to A, the mass M needs to lie inside the upper interval. We shall also verify the useful 
identity from Eq. (2.20) between the discontinuity and the imaginary part right above 
the branch cut, 

Disc S(M) = 2* Im E(M + ie) , (4.33) 
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by applying the Sokhotski-Plemelj theorem to the self-energy's integral representation: 



1 f A . n 2 



\Ai 2 ~+m 2 " (4(u 2 + m 2 ) - M 2 - ie) 



1 /- A , u 2 
du 



2vr 2 7 4V« 2 + m 2 u 2 _ p( m 2 _ Ml ] _ /( 



1 /- A J u 2 
dn 



2^7o ^2 + m 2 ^ + ^ _ m 2} u _^_ m 2_ ie 



(4.34) 



M 2 2 

^ — m 2 



ImS(M) = V , 2m < M < 2\J A 2 + m 2 , (4.35) 



which indeed yields the same result as in Eq. (4.32) after multiplying with 2i. Note 
that this result is valid on the positive real axis, thus it gives us information with which 
quadrant the branch cut is continuous (it is continuous with the first quadrant). 



We now would like to use the final result of the self-energy from Eq. (4.26) to arrive 
at the same result. For that purpose one has to rewrite the inverse tangent with the 
well-known decomposition 



ln(l - iz) - ln(l + iz) , (4.36) 



i 

arctan z = - 
2 

which is especially used in computational software | 70| . The two logarithmic terms 
determine the range and branch cut structure of the inverse tangent fully correctlyj^] 
This form is actually equivalent to 

arctantuz 

arctan z = , (4.37) 

i 

recommended by Kahan (Tl] . Let us guess what will happen along the positive real axis: 
The ratio of negative complex root function and single variable z in the expression for 
the self-energy will give a branch cut for z > 2m with an easily obtained value for the 
discontinuity, while we expext another conribution from the inverse tangent for values 
larger than z = 2-\/ A 2 + to 2 . If this part of the overall cut is somehow totally cancelled 
out it is reasonable to assume the inverse tangent to be responsible for this effect. An 
explicit calculation illustrates this aspect (the logarithmic constants are not important 



4 This statement is not true if one uses 

arctan z = — i In I (1 + iz ', 



1 + z< 



With this representation the inverse tangent would be well-defined left from the upper part of the 
cut and right from the lower part. 
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and were already subtracted): 

DiscS(M) = £(M + ie) — E(M — ie) 
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8vr 2 (M 
In ( 1 + i 



IE 



It I 
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8vr 2 M 

-In | 1 
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In 1 
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ie 



AM 



+ 



\/A 2 + m 2 VM 2 - Am 2 
WM 2 - Am 2 i 



8ir 2 M 
In I 1 + 



In 

AM 



+ ie 
AM 



VA 2 + m 2 VM 2 - 4m 2 



ie 



VA 2 + m 2 VM 2 - Am 2 



+ ie 



M 2 



— m A 



2m< M < 2\/A 2 + m 2 . 



(4.38) 



47rM 

Here, the discontinuity vanishes for M < 2m because the real parts of all radicands after 
the second equality stay positive, so we are beyond the branch cuts of the roots and 
logarithms. The self-energy above and below the positive real axis consequently equal 
each other in the limit of e — > + and subtract to zero. In order to get a discontinuity 
from the roots outside the arguments, we also need the logarithmic cuts to contribute - 
the real part of their arguments has to be negative or zero. Only the inequality 

AM 

1 ^= < (4.39) 

VA 2 + m 2 VM 2 -Am 2 ~ 

has real positive solutions, namely 2m < M < 2\/ A 2 + m 2 . We have plotted the discon- 



tinuity and the numerical difference of T,(E) across the real axis in Fig. 4.5 



In the end, we can perform the analytic continuation of the self-energy down into the 
second Riemann sheet (this is true for the fourth quadrant): 



S n (z) = E(z) + 



m^ 



Airz 



(4.40) 
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Disc E( M) 
0.04 h 



Figure 4.5: Discontinuity of the self-energy (in units of i): numerical difference (dark blue) 
and analytic result (red). 



and as a consequence there the propagator can be written as 

1 



An(z) 



1 



M 2 + 



E(*) + 



4ttz 



(4.41) 



which gains the same branch cut structure in the complex M-plane for gs^ / as 
the self-energy. Although we will work with the inverse propagator and are thus only 
concerned with the self-energy, one can write down, for the sake of completeness, the 
analytic continuation of the propagator itself by using the spectral function. Here, one 
should keep in mind that the starting point in the literature is the complex p 2 -plane: 



Disc As(p 2 



A s (p 2 + ie)-A s {p 2 
k Jo 



ie) 



P(s 2 



J> 

„2 / 2\ 



ds p{s ) T - i 

roo 



s z + ie p z 
-2ie 



)€ 



a 2\2 



+ e 2 



2i 

n Jo 
-2ip(p 2 



ds 2 p(s 2 )tt5(p 2 - s 2 ) 
i = 2ilmA s {p 2 + ie) . 



(4.42) 



We of course recognize the useful identity (2.20) and remember the imaginary part to 



be evaluated right above the branch cut. Because we work in the complex M-plane, the 
spectral function has to be changed by a variable transformation, 



P(s 2 



ImAs(s 2 + ie) -)• d(s) = 2sp(s) = -2s Im A s (s + ie) 



(4.43) 
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where the spectral function of the resonance S is from now on denoted as d(s). Finally, 
the propagator on the second sheet is 

A u (z) = A s (z)-i^- . (4.44) 

The full interacting propagator belongs for the simplest case gs<f><t> = to a stable 
particle with mass Mq in its rest frame: 

A S (M) = — —= . (4.45) 

We have applied our general definition of a stable particles (real) mass as its propagator 
pole. The pole position of course gives a mass which is slightly shifted into the lower 
half plane of the complex M-plane (-y/p 2 = M — >■ z = x + iy): 

Mp i e = M -ie. (4.46) 

By switching on the coupling there is a contribution from the (complex-valued) self- 
energy to the pole position in exactly the same way as for the non-relativistic Lee model. 
The new pole position, as the zero of the denominator, is determined by a system of two 
equations 

x 2 -y 2 -M 2 + <? 2 ^Re£n(x + iy) = 0, 

2ixy + g 2 SH Im S n {x + iy) = , (4.47) 
and we will search for numerical solutions of 

z 2 -M 2 + 5 2 ^S„(z) = (4.48) 

in the complex M-plane. 



4.2.2 Couplings g s<M> G [0.1, 1.0] 

We need to add some words concerning the parameters. The choice of the cutoff A = 
1.5 GeV was already motivated by the energy and mass scales of the effective theory 
presented here. Nevertheless, there is still one free parameter in our model, namely the 
bare mass Mq. In a fundamental theory (where A — > oo) such a quantity can neither 
be obtained from field theory, nor from any experiment and has to be treated in the 
framework of mass renormalization. This is different from our point of view: the bare 
mass is model dependent and was for example taken by Tornqvist |29| to be a free 
parameter (that is fixed after fitting to experimental data), as it will be in this work. If 
we interpret this number as somehow related to the physical mass of the resonance S, 
namely its mass for a vanishing interaction, then it is modified by quantum fluctuations 
of the mesonic loops. The renormalized Breit-Wigner mass Mbw of the resonance is 
defined like in the first chapter as the real zero of the propagator's real part: 

M 2 W - Mq + gl^ Re S(M BW ) = . (4.49) 
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We fix the number by Mbw = 0.6 GeV and thus every resonance pole in the second 
Riemann sheet will describe a particle with this value of the Breit-Wigner mass (this 
choice was made before the PDG update on the c-meson occurred). Such a constraint 
was already introduced by Achasov and Kiselev in order to take into account finite- 



width corrections when studying the propagators of light scalar mesons 72 . Since 



ReX(0.6 GeV) > 0, the renormalized mass is smaller than Mq. In the following subsec- 



tions we search the complex M-plane for numerical solutions of Eq. (4.48) by varying 
the coupling g s ^ in steps of Agg^. The condit ion for deciding whether a solution Zp Q \ e 
marks a singularity or not shall be again Abs (Aj^ 1 (z po i e )) < 1CT 10 (GeV) 2 . 

Similar as in the case of the non-relativistic Lee model, it seems the spectral function 
d(M), defined as the negative imaginary part of the propagator right above the branch 
cut (see first chapter), is normalized only for couplings in the approximate interval 



[0.2 GeV, 4.5 GeV], see Fig. 4.6 For smaller gs6<f> there would be a steep descent at 



0.0 £ S JI GeV II 

5 10 15 20 25 

Figure 4.6: Incomplete normalization of the spectral function d(M) for different couplings. 

the beginning which we ignore since, if the coupling tends to zero, the spectral function 
becomes a delta distribution function due to the fact that the pole approaches the real 
axis: 

d{M) = 2M( e + ^ImS(M + i e )) 

( M 2 - Ml + g\ H Re E(M + ie)) 2 + (e + g S(f>(f> Im £(M + ie)) 2 
gs^o 2Me 

(M 2 - M 2 ) 2 + e 2 

2Mtt5(M 2 - M 2 ) . (4.50) 

As we know from the previous chapter, the numerical integration ^ J °° ds d(s) is not 
possible in the vicinity of a simple pole. This was also true for poles on the real axis in 
the first sheet: we expect additional poles in the first Riemann sheet once the spectral 
function is decreasing for sufficiently large couplings. This time though, one of the new 
poles will not emerge right from the beginning. 

In general, d(M) becomes broader for larger couplings, because the complex mass- 
pole moves deeper down into the complex plane (i.e., the second Riemann sheet). It 
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Im z 



0.2 0.4 



T 



0.8 1.0 1.2 1.4 



Re z 



-0.015 - 



Figure 4.7: S4>4>-mode\ poles (dark blue dots) for gs^ £ [0.1 GeV, 1.0 GeV] in the second 
Riemann sheet with Agg^ — 0.05 GeV, where the gray line marks the branch 
cut. 



shrinks very quickly the more the pole descends. The pole trajectory within the regarded 



coupling interval [0.1 GeV, 1.0 GeV] is shown in Fig. 4.7 In contrast to the results in 



chapter 3, one clearly sees another pole coming from the origin and heading to the first 
branch point at z = 2m J^] Such a phenomenon was not observed in the non-relativistic 
Lee model. Note that there is no possible way of finding an analytic solution of Eq. 



(4.48), not even for purely real z-values. Yet, we can solve for gs<j><j> and, after taking 



just positive solutions into account, perform the limit 



gs<jxf>(z = 0-ie) 




M 2 



it 



Re S(M BW ) - S(M - ie) - Disc S(M - ie) 



(4.51) 



This result is due to the divergent behaviour of the discontinuity at the origin in the 
second sheet. 

A summary of pole positions and other relevant informations for different intervals 



can be found in Tab. 4.2. The adopted notation is very similar to the one for the Lee 
model (for a better understanding see Tab. 4.1): The pole mass es X( m ) P oi e and X(o) po ie 



are the real parts of the complex solutions found by solving Eq. 
the real part of the complex solution in the first Riemann sheet, w 



4.48), while x 



(l)pole 



IS 



rere we have replaced 

the self-energy T,u(z) in Eq. (4.48) by 'S(z). The corresponding Breit-Wigner masses 
are the solutions of the latter equation using only the real part of the self-energy on 
the first sheet (they are all real by construction). These are of course the zeros of the 
propagator's denominator on the real axis. The numerical maxima of the continuous 
part of the spectral function are denoted as xn-\ max and X( r ) max , whereas ^average is the 



5 Compared to our strategy in the non-relativistic Lee model, we have chosen this time the accuracy of 
the starting values on the lattice in such a way that poles on the real axis can be found the moment 
they really exist. 
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name label determining equation description 



resonance pole mass 


^(mjpole 


z 2 - 




Real part of solution z. 


gap pole mass 


^(OJpole 


z 2 - 


-M$ + g 2 s ^ 11 (z) = 


Real part of solution z. 


left pole mass 


*£(l)pole 


z 2 


-Mi + g 2 s ^(z)=0 


Real part of solution z. 








Nevertheless, z is purely real. 


Breit-Wigner resonance mass 


I (m)BW 


x 2 - 


M 2 + ff | # ReE(x) = 


The solution x is purely real. 


Breit-Wigner left /right mass 


^(1)BW 


x 2 - 




The solution x is purely real. 


average mass 


^average 


X = 


1 r 2\/A 2 +m 2 , „ v 


Numerically evaluated. 


maxima of spectral function 


^(l)max 




max d(s') 


Solution is x = M max . 
Found numerically. 



Table 4.1: Description of masses for the considered S4>4>- model. The general complex so- 
lution of an equation is denoted as z = x + iy. 



result of the numerically evaluated integral 



^average 



7T 



2 v / A 2 +m 2 



ds sd(s) 



(4.52) 



2m 



The last column of Tab. 4.2 simply gives the normalization value of the spectral function 



(which indeed is normalized to one, as will be shown in subsection 4.2.7). All numbers 
after the fourth digit are dropped. 



4.2.3 Couplings g SH> 6 [1.0,2.8] 

The resonance pole descends deeper into the lower half plane and turns right, while the 
pole near the origin, from now on denoted as the gap pole, keeps approaching the first 
branch point with rising parameter speed. No poles are found in the first sheet. 



Im z 



- 0.05 



-0.10 



Re z 



Figure 4.8: 5(/>0-model poles for gs^ € [1.0 GeV, 2.8 GeV] in the second Riemann sheet 
with Ag s<p< p = 0.05 GeV. 
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9S<jj<t> ^(OJpole ^(ljpole a; (m)pole I (0)BW ^(lJBW I (m)BW ^max ^average Md Norm 



0.1 


- 
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0.6000 


- 


- 


0.6000 


0.6000 




- 





.6001 




- 


0.2 


0.0012 




- 


0.6000 


- 


- 


0.6000 


0.6000 




- 





.6007 




- 


0.3 


0.0027 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6009 





.6017 





.9999 


0.4 


0.0049 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6016 





.6030 


1 


.0000 


0.5 


0.0076 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6025 





.6047 


1 


.0000 


0.6 


0.0110 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6036 





.6068 


1 


.0000 


0.7 


0.0151 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6050 





.6093 


1 


.0000 


0.8 


0.0197 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6065 





.6121 


1 


.0000 


0.9 


0.0255 




- 


0.6000 


- 


- 


0.6000 


0.6000 





.6082 





.6153 


1 


.0000 


1.0 


0.0309 




- 


0.6000 


- 


- 


0.6000 


0.6001 





.6102 





.6188 


1 


.0000 


1.2 


0.0445 




- 


0.6002 


- 


- 


0.6000 


0.6002 





.6147 





.6269 


1 


.0000 


1.4 


0.0605 




- 


0.6004 


- 


- 


0.6000 


0.6005 





.6199 





.6364 


1 


.0000 


1.6 


0.0788 




- 


0.6008 


- 


- 


0.6000 


0.6009 





.6260 





.6471 


1 


.0000 


1.8 


0.0991 




- 


0.6014 


- 


- 


0.6000 


0.6016 





.6329 





.6591 


1 


.0000 


2.0 


0.1212 




- 


0.6025 


- 


- 


0.6000 


0.6025 





.6405 





.6722 


1 


.0000 


2.2 


0.1444 




- 


0.6042 


- 


- 


0.6000 


0.6039 





.6489 





.6864 


1 


.0000 


2.4 


0.1680 




- 


0.6068 


- 


- 


0.6000 


0.6058 





.6580 





.7016 


1 


.0000 


2.6 


0.1911 




- 


0.6108 


- 


- 


0.6000 


0.6084 





.6679 





.7178 


1 


.0000 


2.8 


0.2128 




- 


0.6166 


- 


- 


0.6000 


0.6119 





.6785 





.7349 


1 


.0000 


3.0 


0.2318 




- 


0.6247 


- 


- 


0.6000 


0.6166 





.6898 





.7528 


1 


.0000 


3.2 


0.2476 




- 


0.6354 


- 


- 


0.6000 


0.6228 





.7018 





.7715 


1 


.0000 


3.4 


0.2597 




- 


0.6487 


- 


- 


0.6000 


0.6308 





.7145 





.7909 


1 


.0000 


3.6 


0.2684 




- 


0.6646 


- 


- 


0.6000 


0.6408 





.7279 





.8110 


1 


.0000 


3.8 


0.2740 




- 


0.6828 


- 


- 


0.6000 


0.6531 





.7418 





.8317 


1 


.0000 


4.0 


0.2770 




- 


0.7031 


- 


- 


0.6000 


0.6678 





.7565 





.8529 


1 


.0000 


4.2 


- 




- 


0.7251 


- 


- 


0.6000 


0.6848 





,7717 





.8748 


1 


.0000 


4.4 


- 





.2773 


0.7485 


0.3024 


0.2773 


0.6000 


0.7040 





,7793 





.8971 


1 


.0000 


4.6 


- 





.2753 


0.7732 


0.3288 


0.2753 


0.6000 


0.7252 





,7884 





.9198 


1 


.0000 


4.8 


- 





.2723 


0.7991 


0.3561 


0.2723 


0.6000 


0.7481 





.7995 





.9430 


1 


.0000 


5.0 







.2685 


0.8259 


0.3840 


0.2685 


0.6000 


0.7725 





.8123 





.9666 


1 


.0000 


5.5 







.2563 


0.8966 


0.4568 


0.2563 


0.6000 


0.8388 





,8510 


1 


.0271 


1 


.0000 


6.0 







.2416 


0.9718 


0.5335 


0.2416 


0.6000 


0.9109 





,8973 


1 


.0895 


1 


.0000 


6.5 







.2251 


1.0507 


0.6138 


0.2251 


0.6000 


0.9875 


0. 


,9491 


1 


.1535 


1 


.0000 


7.0 







.2075 


1.1329 


0.6972 


0.2075 


0.6000 


1.0679 


1 


.0052 


1 


.2189 


1 


.0000 


7.5 







.1887 


1.2182 


0.7838 


0.1887 


0.6000 


1.1518 


1. 


.0646 


1 


.2854 


1 


.0000 


8.0 







.1688 


1.3064 


0.8732 


0.1688 


0.6000 


1.2390 


1 


,1267 


1 


.3529 


1 


.0000 


8.5 







.1475 


1.3973 


0.9656 


0.1475 


0.6000 


1.3297 


1 


,1909 


1 


.4212 





.9999 


9.0 







.1241 


1.4909 


1.0608 


0.1241 


0.6000 


1.4238 


1 


,2569 


1 


.4902 





.9999 


9.5 







.0972 


1.5870 


1.1590 


0.0972 


0.6000 


1.5218 


1 


,3241 


1 


.5599 





.9998 


10.0 







.0623 


1.6857 


1.2603 


0.0623 


0.6000 


1.6239 


1 


,3917 


1 


.6301 





.9994 


10.2 







,0420 


1.7258 


1.3018 


0.0420 


0.6000 


1.6660 


1 


,4186 


1 


.6584 





.9991 



Table 4.2: Selection of masses (in units of GeV) for the S '^-model with A = 
1.5 GeV, M BW = 0.6 GeV and m = 0.139 GeV. 
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Im 
0.10 r 



Re z 



Figure 4.9: S (jxp-model gap pole for gs<p<j> S [1-0 GeV, 2.8 GeV] on the real axis of the 
second Riemann sheet with Ags^ = 0.05 GeV. 



4.2.4 Couplings g s<M> G [2.8, 5.0] 

For this range of the coupling we are faced with an interesting behaviour. As we turn 
up gscfxf)-, for some value larger than 4 GeV the gap pole hits the first branch point at 
z = 2m and vanishes, in fact it slips into the first Riemann sheet. Although this strange 
pole has no physical interpretation, its existence below threshold lasts as long as gs<t><t> 
is non-zero, because then the self-energy, its real part at the point z = Mbw an d the 
discontinuity are all three purely real and positive right below the real axis. In this case 
both the numerator and the denominator under the root of the master solution 



gs<t><t>(M 



M 2 



Kw 



ie 



Re S(M BW ) - E(M - ie) - Disc E(M - ie) 



(4.53) 



yield negative values and give a positive one for gs<j><f> ■ It is very important to realize 
that the pole must lie below the real axis, so the infinitesimal number ie is crucial. 
Otherwise, the root inside the expression for the discontinuity would be evaluated in the 
wrong quadrant and there would be an imaginary solution for the coupling. Indeed 
the gap pole leaves the second sheet after arriving at the first branch point for 



5S^(2m - ie) 



V ReS(M BW ) 



(2m — r\Y 



It 



S(2m — rj 



Disc S(2m — rj — ie) 



4m 2 



ReS(M BW ) -S(2m) 



(4.54) 



which numerically gives a value of gs<jxj> = 4.209 GeV. In order to make this clear, one 



has to evaluate the limit analytically by extracting the real part of the self-energy (4.26) 
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Figure 4.10: Scfxfi-mode\ poles for gs<f><f> S [2.8 GeV, 5.0 GeV] in the second Riemann sheet 
with Ag s<j)l p = 0.05 GeV. 
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Figure 4.11: Left Scfxfi-model pole for 
first Riemann sheet with 



£ [2.8 GeV, 5.0 GeV] on the real axis of the 
= 0.05 GeV. 



by performing the principal value integral on the integral representation from Eq. (4.22 ) 



We start with Eq. fl4.34[ ): 

Re S(M) : 



2tt 2 
1 

2^2 



V / du 



u 



Ml_ m 2 



\/m 2 + m 2 (4(ii 2 + m 2 ) - M 2 ) 

m z —T) „,2 



du 



ij- 



+ 



Vu 2 + m 2 (4(n 2 + m 2 ) - M 2 ) 
„,2 



du 



^M2_ m2+r] x/n^Tm 2 (4(u 2 + m 2 ) - M 2 ) 



(4.55) 
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Only for 2m < M < 2\/A 2 + m 2 the whole integral gives another result than the already 
known self-energy. By using the same indefinite integral as during the calculation of the 
self-energy in the beginning of our investigation and performing the indicated limits, we 
arrive at a hardly suprising result: 



ReS(M) 



M 2 



8vr 2 M 



'J^-m 2 



arctan 



AM 



1 



4vr 2 



+ m 2 \/Am 2 - M 2 J 8?r 2 
arctan(— i + irj) — arctan(— % — irj) 



In 



m 



A + VA 2 + m 2 



E(M) - ilmS(M) 



(4.56) 



Obviously, the imaginary part (which comes from the known two terms already found 
on the previous pages) is subtracted and leaves the real part we were looking for. We 
emphasize that although this last expression might look trivial, it is in fact an important 
result. We are now able to give explicitly the real value of the self-energy right above the 
real axis and in principle could also do the analytic continuation down into the second 
sheet with the help of the spectral function. However, with this result the denominator 
in the upper limit 



4m 2 



Re£(M BW ) -S(2m) 



(4.57) 



can be written explicitly, but it again yields a transcendental equation, so we have to 
perform the numerical result for the coupling. The unphysical gap pole finally emerges 
on the first sheet left from the first branch point, as can be seen in Fig 4.11 This is 
surely true since the master solution in the first sheet leads to 



'2m — ie) 



(2m - V ) 2 - M 2 W 



ie 



v™+ V Re S(M BW ) - S(2m -rj-ie) 



4m 2 - M 2 W 



ReS(M BW ) - S(2m) 



(4.58) 



thus the same result as in the second sheet. Furthermore, this equation has no real 
solution for g < 4.209 GeV, yet they are purely real for larger values. We may show the 
pole appearing in the first sheet simply by considering the spectral function for M < 2m: 



d{M) 



2M(e + g 2 ImS(M + ze)) 



( M 2 _ M 2 + g 2^ Re E(M + ie) y + ( £ + g 2^ Im E(M + ie) y 

2Me 



[M 2 -M 2 W 
2Mtt5(m~ 



ReS(M BW ) -ReS(M + ie))] 2 + e 2 



ReS(M BW ) - ReS(M) 



(4.59) 
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where the argument of the delta distribution function naturally becomes the master 
solution in the first sheet (the self-energy is real below threshold). From here, we also 
see that the Breit-Wigner mass £(i)bw f° r the left pole equals the corresponding pole 
mass. 



4.2.5 Couplings g SH> G [5.0, 10.0] 



The resonance pole in the second sheet is now the only one there. It descends deeper 
down in the lower half plane and its imaginary part becomes comparable to its real 
part. In the first sheet, the new pole speeds up and approaches the origin. Its existence 
has dramatic effects on the normalization of the spectral function. It seems we can 
understand the lack of normalization by finding this new pole, but for completeness 
reasons one should pay attention to the second branch point at z = 2\/A 2 + to 2 . We do 
not expect a third pole slipping from the branch cut into the complex plane, because 
the analytic continuation of the discontinuity into the second Riemann sheet provides 
us with a branch cut heading there to infinity. In spite of that, the limit 



(o / A2^ 2 ■ ^ v (2^K 2 +m 2 +r 1 ) 2 -M 2 w -te 

QsssV^V A 2 + m z — le) = hm \ 

^ ' r,^o+ V Re S(M BW ) - S(2VA 2 + to 2 + rj - ie) 



4A 2 + 4m 2 - Mr 



bw 



V ReS(M BW ) - X^^/A 2 " 
= , (4.60) 

is obtained, where it is positive for every z > 2\/A 2 + m 2 (since the inverse tangent in the 
self-energy is evaluated at a branch point, i.e., — i arctan(— i) = — oo, and stays imaginary 
for M — > oo, while the prefactor root function is imaginary, too). So, we just made clear 
that there is another pole coming from the second branch point at z = 2\/A 2 + m 2 right 
from the beginning, similar to the non-relativistic Lee model. This follows also from the 
spectral function for M > 2a/ A 2 + to 2 : 

d{M) = 2M( e + 4lmS(M + i£ )) 

(M 2 - M 2 + g 2 H Re E(M + ie)) 2 + (e + ff 2 ^ Im S(M + ze)) 2 

2Me 

[M 2 - M 2 W - g 2 ( Re E(M BW ) - Re E(M + ie))] 2 + e 2 



= 2Mvr ( 5(M 2 -M2 w -« 7 2 ^(ReS(MBw)-ReS(M))) . (4.61) 

Again, the argument of the delta distribution function requires the master solution in 
the first sheet. The important point here is that this additional pole is due to the finite 
cutoff, thus it is not a physical effect - one either needs to perform the limit A — > oo 
(and subtract the divergence) , or use smooth regularization functions to make the second 
branch point disappear at complex infinity (where it has to be, see first chapter). That 
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is why we completely disregard this polejand do not list its masses in Tab. 4.2 
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Figure 4.12: S^-model resonance pole for gs<f><f> € [5.0 GeV, 10.0 GeV] in the second 
Riemann sheet with Ags^ = 0.05 GeV. 
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Figure 4.13: Left 500- model pole for gsd>6 S [5.0 GeV, 10.0 GeV] on the real axis of the 
first Riemann sheet with Aqsss = 0.05 GeV. 



4.2.6 Couplings g s<M> 6 [10.0, 10.5] 

On top of that, we notice another remarkable behaviour of the left pole. It speeds up 
and approaches the origin but seems to vanish for some specific value of the coupling. 



6 The numerical consequences are not worth mentioning: On the one hand the pole can only be found 
on the lattice simultaneously with tachyonic excitations (for huge couplings), on the other hand its 
renormalization constant is so small that it can be neglected for the presented range of the coupling. 
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Then, there is another pole emerging on the imaginary axis of the complex M-plane. 
We find the critical value by 



9SH>(° ~ ie ) 




M 2 - M 2 W - %e 
Re£(M BW ) -S(M-ie) 
Mbw 



A /S(0-ie)-ReE(Af B w) 
A /ReE(M BW )-S(0-ie) 



Re£(M BW ) < S(O-ie) 
Re£(M BW ) > S(0-ie) 



(4.62) 



where in our case the first is true with gs<f><f> = 10.371 GeV. The left pole will consequently 
arrive at the origin for a positive value of the coupling and the question now is whether 
the tachyonic pole on the imaginary axis in the lower half plane manifests at exactly the 
same gsss- The following limit clarifies this: 



gsssiv ~ iQ ) = lim 



(r? + iyf - M 2 W 
y^o- V Re S(Mbw) - £(77 + iy) 



- v /S(j ? -iO)-ReE(M B w) 
^/RoE(Af B w)-S(»?-iO) 



Re E (M BW ) < Efa - i0) 
Re£(M BW ) > Efa-iO) 



(4.63) 



The inverse tangent inside the expression for the self-energy is purely imaginary on the 
negative imaginary axis and has absolute magnitude smaller than one, because the ar- 
gument is smaller than one. Together with the prefactor, which is also purely imaginary, 
this combines to a negative real number subtracted by a positive logarithmic term. In 
the end, the self-energy is real and positive on the negative imaginary axis with the 



same limit as in Eq. (4.62). It declines along the negative imaginary axis and so the 
radicand in the first of the two above equations stays positive while increasing. This 
means that the purely imaginary pole keeps heading to complex infinity for increasing 
coupling. Hence, the tachyonic pole shows up when the left pole vanishes. Note that all 



these statements are strictly valid only for the first case in Eq. (4.63). It clearly depends 
on the choice of the Breit-Wigner mass Mbw an d although we possess an expression for 
the real part of £(M) we cannot perform further studies in that direction in virtue of 
the transcendental character of our equations. The numerical bound in order to obtain 
tachyonic excitations is Mbw = 0.4984 GeV. Since the whole model has then become 
unphysical due to the tachyonic pole on the imaginary axis, it of course makes no sense 
to study the pole's behaviour for larger couplings. 

4.2.7 Spectral function and its normalization 

Let us now focus on the spectral function d(s) and especially on its normalization. As 



was shown in Fig. |4.6| the spectral function is not normalized when using the numerical 

ds d(s) , (4.64) 



result of the integral 

1 



7T 
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Figure 4.14: Scfxfi-mode\ poles for gs<p<p € [10.0 GeV, 10.5 GeV] in the first Ricmarm sheet 
with Ag S4><j) = 0.01 GeV. 



simply because the additional simple poles are not seen by numerics. It is possible to 
take them into account by splitting the full propagator A$(M) into the contributions of 
two single-particle propagators and the remaining continuum part above the branch cut 
(in particular, the right pole will be ignored completely; its renormalization constant is 
very small and gives no significant contribution): 



A S (M) 



vr J M 2 — s 2 + ie 

-^(l)pole 1 



2VA 2 +m 2 



ds 



d(s) 



M 2 — s 2 + ie 



(4.65) 



The renormalization constant ^(i) po i e is calculated by expanding the inverse of the full 
propagator in a Taylor series around the pole at x = £n/ r ) P oie i n first order: 



x 2 -M%+g%^nx) 



(x 2 - Ml + g 2 s ^(x))\ x=x 



(l)polc 



+ — [x 2 -M 2 +q%^Y,(x. 
dx 



2 dReS(x) 
2aJ(i)pole + 9su fa. 



C (l)polG 



( x ~ Z(l)pole) 



(x - XQpoie) 



[ '(l)polc/ 



(4.66) 



The same argument as for the non-relativistic Lee model can be adopted here: The 
simple poles and the branch cut lie separated on the real axis at different masses, so 
the renormalization constants for each single-particle pole can be extracted from their 
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residues: 



lim (x 

»->^(l)pole 



"(l)pole. 



•A 



s(x 



lim (x 

^^(ljpolo 



x (l)poleJ ' 



J (l)polc 



'(l)pole 



+ ie 



Res(A s (a;), x=x (1)polc ) 

+ {other pole and branch cut} 
= ^(l)pole • 

After inserting the above expansion into the left side, we finally get: 



(4.67) 



y (l)polc 



1 + 



9sh> dReE(x) 



2x 



(l)polc 



dx 



c (l)polc. 



(4.68) 



The normalization condition thus can be written as 



l = Z, 



(l)polc 



+ 



ds d(s) . 



2 m 



(4.69) 



If we perform our calculation including the renormalization constants, the spectral func- 



tion is normalized over the full range of gs4><p ■ The plot in Fig. 4.15 shows a completely 
fulfilled normalization condition for our choice of parameters (the points used are those 



from Tab. 4.2). For the sake of completeness, we provide a compilation of selected 



spectral function plots in Appendix F No analytic study is done since the expressions 
for the real and imaginary part of the self-energy are far too complicated to handle. 
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Figure 4.15: Numerical verification that ^ J °° ds d{s) = 1 (red) by including the single 
poles of the first Riemann sheet, and numerical contribution of the continuous 

part ~ f 2 J^ +m ds d(s) (dark blue). The necessity of the delta distribution 
functions is evident. 
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4.2.8 Decay width(s) 

It is reasonable to have a look at the full decay width of the resonance as two times the 
negative imaginary part y( m ) po io = — r/2 of the complex pole descending on the lower 
half plane of the second Riemann sheet. As was argued long time ago by Matthews 
and Salam |73| , the spectral function can be interpreted as a mass distribution of the 
unstable particle S. Remember the first chapter, where we have introduced the spectral 
function: given an interacting single-particle state, the spectral function p(s 2 ) describes 
the probability of finding a free state with invariant mass s 2 . This in mind we already 
applied such a view when we calculated the average mass x average in the non-relativistic 



Lee model and the Scjxfi-model, too. In Ref. 67 , Giacosa and Pagliara have given an 
intuitive argument for the correctness of the upper interpretation. In agreement with 
the authors point of view, we can define an average decay width by the expression 

1 f°° 

Taverage = ~ / ds d(s)T tTee (s) . (4.70) 
n Jim 

The first two panels of the next figure show the plot of the full width T (red) in com- 



parison with the tree-level result (4.24) for the constant Breit-Wigner mass Mbw (blue, 
dashed) and the average width (blue) as a function of the coupling. For small couplings 
9s^><f> < 1 GeV there is actually no mentionable difference between those three, and they 
all vanish for gs<jxj> — > 0. The situation changes completely for higher values: the full 
width becomes larger than the tree- level result, which again is larger than the average 
width. While this is true only for the latter two over the whole range of gs(j><j>, the tree- 
level result passes the value of the full width at gs^ = 5.622 GeV and is nearly two times 
the full width when arriving at the end of the interval. It therefore cannot be regarded 
as a good approximation for large couplings, though it gives 'the best fit' for all values 
below the intersection point. On the other hand, the average value of the decay width 
is far too small in the intermediate range and large-coupling limit, but nevertheless the 
difference to the full width is 'only' 30%. 

As an interesting consideration, one can study other decay widths as a combination of 
spectral function properties and the tree- level result. In the second part of Fig. 4.16 we 
provide plots of the full (red) and average width (blue) , and plots of the tree- level results 
for each pole mass £( m ) pole (green, dashed), each bare mass Mq (magenta, dashed), each 
maximum value x max of the spectral function (brown, dashed) and each average mass 
^average (black, dashed). All of them may differ in several ways, but overall they are 
comparable and still yield too small values. The black dashed curve has the smallest 
error for large gs</xf,, but not for intermediate ones. 
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Figure 4.16: Comparing decay widths for the with sharp cutoff (and varying 

scales): the first two panels show the full width T (red) in comparison with 
the tree-level result ( 4.24 1 for the constant Breit-Wigner mass Mbw (blue, 
dashed) and the average width (blue) as a function of the coupling. The lower 
two panels are plots of the full (red) and average width (blue), and of the 
tree-level results for each pole mass X( m ) po ic (green, dashed), each bare mass 
Mo (magenta, dashed), each maximum value x max of the spectral function 
(brown, dashed) and each average mass average (black, dashed). 
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4.2.9 Concluding remarks 

The present analysis is based on a fixed Breit-Wigner mass of Mbw = 0.6 GeV. This 
is why the pole trajectory on the second Riemann sheet cannot be taken as the path of 
the resonance pole for different values of the coupling, but it outlines the pole position 
when requiring a constant Breit-Wigner mass. Some relevant points before we continue: 

1. We have shown that a single seed pole for the u-meson entails two other emerging 
poles, once the interaction term in the Lagrangian is turned on. One pole appears 
in the first sheet and can be discarded as non-physical, induced by the finite branch 
cut, while the other pole is dynamically generated in the vicinity of the origin in 
the second sheet. This pole moves to the branch point at threshold and slips 
through the cut onto the first sheet, where it has to be taken into account in the 
normalization condition of the spectral function. It can be interpreted as a bound 
state. The whole model becomes unstable due to the occurrence of tachyonic 
excitations for large couplings. 

2. Although the Breit-Wigner parameterization reproduces very well the mass (and 
decay width) of the resonance pole in the small-coupling regime, this is not true 



anymore for intermediate values of gs<f><f>, compare Tab. 4.2 Above this regime, 
the mass is completely detached from the pole mass, while the average and bare 
mass, and the maximum value of the spectral function are sometimes in acceptable 
agreement with the pole mass. 

3. The tree- level result for the decay width cannot be regarded as a good approxima- 
tion for large couplings, yet it gives 'the best fit' in the intermediate regime (with 
deviations up to 100 MeV). 

4.3 Smooth cutoff 

4.3.1 Analytic structure of the propagator 

A regularization function with sharp cutoff behaviour can only be motivated for reasons 
of simplicity. Fortunately, as long as /a(|<i|) depends only on the magnitude of the 
three-momentum q, any smooth function modifies the self-energy integral 

E(M) = ifd» , " 2/ » f (4.71) 

V ' 2vr 2 y Vu 2 + m 2 (4(u 2 + m 2 ) — M 2 — it) 



as a multiplication factor that has to be included during the it-integration. In Ref. 67 
a regularization function 

/A(|q|) = (4.72) 
was chosen, yielding an exponentially decreasing interaction strength for increasing dis- 



tance of the two ^-particles (see non-local Lagrangian (4.23)). Here, the analytic result 
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for S(M) reads: 

S(M) 



A 3 



4(4A 2 - 4m 2 + M 2 ) 2 7r 2 



A(4A 2 - 4m 2 + M 2 



A 2 



8A 2 - 4A 2 m 2 - 4m 2 + m 2 M 2 



(m 2 



A 2 ) I 



arctan 



A 2 



A 



8A\/4m 2 



M 2 



M 



arctan 



M 



\/4m 2 



M 2 



(4.73) 



Considering only positive values = M > 0, the imaginary part of E(M) will be 

non-zero above threshold and falls rapidly for M —> oo, see Fig. |4. IT This gives an 
infinite branch cut. In comparison to the self-energy for the sharp cutoff, all additional 



terms in Eq. (4.73) are single- valued, so it is clearly the root function and the inverse 
tangent in the last line that provide us with the branch cut starting at z = 2m. The 
overall prefactor in the first line then makes the imaginary part of the self-energy fall for 
high M as an effect of the incorporated regularization function. Note that the inverse 
tangent in the last term has two branch points in the complex M-plane, which are the 
same as for the root function: 



(1) 



M 



m c 



M 2 



M 



(4.74) 



(4.75) 



(2) 



M 



V 4m 2 



M 2 
M 



-too 



lim (2m + rj) 



(4.76) 
(4.77) 
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Let us look for the discontinuity by using the integral representation of the self-energy. 
As seen after the calculation carried out in the last section, we also here find the regu- 
larization function to be only a multiplicative factoiQ evaluated at the three- momentum 
u = pstjxfi = \/ M 2 /A — to 2 , namely: 

DiscS(M) = E(Af + ie) - S(M - ie) 
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2^2 



u 



du . , 

" (1 + *) 
1 



Vu 2 + to 2 (4(u 2 + to 2 ) - M 2 - ie) 



Vu 2 + m 2 (4(u 2 + to 2 ) - M 2 + ie) 
„ 2 

du 



7T 
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Vn 2 +TO 2 (l + ^) 2 8 



Ai 2 



on — 



TO 



m/^-to 2 

M 2 _ w 2\2 

l7T.\A[ 1 + 



M > 2to . 



(4.78) 



The useful identity (2.20) between the discontinuity and the imaginary part right above 



the branch cut still holds true, since 



£(M) 
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du 



27r 2 7o Vu 2 + to 2 (1 + g) 2 (4(u 2 + to 2 ) — M 2 — ie) 
1 f A . n 2 1 



2vr 2 
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v / n 2 + m 2 (l + g) 2 (n + 



ImS(M) 



M 2 
4 



TO^ 



87rM^l + ^^- 



f - m 2 ii/- 



M > 2m , 



M 2 o 

to — ze 



(4.79) 



(4.80) 



which indeed yields the same as in Eq. (4.78) after multiplying by Ii. This result is valid 



on the positive real axis, so that the branch cut is continuous with the first quadrant. 



The same expression can be obtained directly from the self-energy function (4.73) 



by using the logarithmic decomposition (4.36) for the inverse tangent. The ratio of 



7 In particular, this would make further studies with more complicated regularization functions easier: 
Instead of working with loop integrals, one can use the general expression 



DiscE(Vp 2 = M) = 



AttM 



/a(ps^) 



and perform a dispersion integral in the complex p 2 -plane. This should simplify the computation in 
most cases. By adopting such a procedure, used by authors when fitting the scattering amplitude 
to experimental data, we automatically satisfy physically correct analytic properties, i.e., we get no 
spurious poles or cuts and have the correct asymptotic behaviour. 
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negative complex root and single variable z will give a branch cut for z > 2m along the 
positive real axis, while the inverse tangent once more just acts as a 'doorman' allowing 
the discontinuity of the root function to appear. Since this part gives support for all 
M > 2m, the overall cut goes to infinity - here lies the difference to the sharp cutoff. 
Explicitly this means (the single- valued terms were already subtracted): 
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(4.81) 



The discontinuity vanishes for M < 2m because the real parts of all radicands after 
the second equality stay positive, so we are beyond the branch cuts of the roots and 
logarithms. The self-energy above and below the positive real axis consequently equal 
each other in the limit of e — > + and subtract to zero. In order to get a discontinuity 
from the roots outside the arguments, we also need the logarithmic cuts to contribute - 
the real part of their arguments has to be negative or zero. Only the inequality 

M 

1- = < (4.82) 

VM 2 - Am 2 ~ V 7 

has real positive solutions, namely M > 2m. We have plotted the above discontinuity 



and the numerical difference of across the real axis in Fig. 4.18 
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Figure 4.18: Discontinuity of the self-energy (in units of i): numerical difference (dark 
blue) and analytic result (red). 



The analytic continuation of the self-energy down into (at least the fourth quadrant 
of) the second Riemann sheet can now be written as 



£n(z) = S(z) + V 4 ^_ 2x2 , (4.83) 



and as a consequence there the propagator reads: 
An(z) = 



1 



1 



z 2 -M% + gl H E(z) + 



^(i + i^) 2 



(4.84) 



4.3.2 Couplings g s<M> G [0.1, 1.0] 



In the following, there will be one difference to the case of a sharp cutoff concerning 
the choice of parameters. The values for m and A will stay unchanged, whereas it 
is interesting to study the poles in the second Riemann sheet under the view of two 
different mass definitions. We have already introduced the (renormalized) Breit-Wigner 
mass Mbw oi the resonance as the real zero of the inverse propagator's real part, 

Ml w - Ml + 5 2 Re S(Mbw ) = o , (4.85) 



which was fixed by setting it to 0.6 GeV, a rough estimate of the u-meson mass before 
the PDG update. In addition to that, we will now investigate the pole structure by 



fixing the bare mass Mq = 0.6 GeV. We will solve Eq. (4.48) numerically by varying 
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the coupling gs4>4> in steps of Ags^, using these parameters. The condition for deciding 
whether a solution 2: po i e marks a singularity is still Abs(Ajj (z po i c )) < 10- 10 (GeV) 2 . 
Again, the spectral function d(M) is normalized for both cases only for couplings 



belonging to an interval more or less [0.2 GeV, 4.5 GeV], see Fig. 4.19 All distortions 
for small gs<jxj> are neglected, since if the coupling tends to zero, the spectral function 
becomes a delta distribution function due to the fact that the pole approaches the real 
axis. Numerical integration of - f °° ds d(s) is not possible when poles on the real axis 




0.0 1 gs^IGeV] 

5 10 15 20 25 

Figure 4.19: Incomplete normalization of the spectral function d(M) for different cou- 
plings. The dark blue curve belongs to the case of fixed Mbw, the red curve 
to the case of fixed M . 



in the first sheet are present. We therefore expect extra poles in the first Riemann sheet 
for sufficiently large gs<j><j>- 

In general, both spectral functions become broader for larger couplings: the resonance 
poles move deeper down into the complex plane (i.e., the second Riemann sheet). The 
spectral functions 'shrink' very quickly the more the poles descend. The pole trajectories 



within the regarded coupling interval [0.1 GeV, 1.0 GeV] are shown in Fig. 4.20 It seems 
that the pole in the case of fixed Mq is faster (concerning gs ^-yelocity) and its path is 
twisted to the left, while the real part of the other one (with fixed Mbw) never falls below 
its starting value. Nevertheless, in both cases one clearly notices a gap pole heading to 
the first branch point at z = 2m. The one for fixed bare mass takes precedence. To catch 
up with the analysis of a sharp cutoff, we have chosen the accuracy of the starting values 
on the lattice in such a way that poles on the real axis can be found the moment they 
really exist and indeed, the whole situation is very similar to the previous investigation. 
Equation (4.48) can be solved for gsss and the usual limits can be studied, whereas we 



will only consider the solutions for fixed bare mass Mq, since all analytic expressions for 
the other case were already found in the last section (see especially Eqs. (4.51), (4.54), 



(4.58), (4.62) and (4.63)). Clearly, just the numerical values will be different. We first 
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look at 



lim 

Af-K)+ 

. 



Ml 



M 2 



re 



S(M - ie) + Disc E(M - ie) 



(4.86) 



This limit is due to the divergent behaviour of the discontinuity at the origin in the second 
sheet. It should be stressed that the infinitesimal number ie is crucial. The coupling 
increases if the mass increases (as an issue of the interplay between discontinuity and 
self-energy) . 

A summary of pole positions and other relevant information for different intervals can 
be found in Tab. 4.3 for fixed Mbw an d in Tab. 4.4 for fixed Mq. As a reminder (see 
also Tab. |4.1|): The pole masses X(- m ) po lc and x 



(0)polc 



are the real parts of the complex 



solutions found by solving Eq. (4.48), while X(i) po i e is the real part of the complex 
solution in the first Riemann sheet, where we have replaced the self-energy £n(z) in Eq. 
(4.48) by S(z). The corresponding Breit-Wigner masses are the solutions of the latter 



equation using only the real part of the self-energy on the first sheet (they are all real by 
construction). These are of course the zeros of the propagator's denominator on the real 
axis. The numerical maxima of the continuous part of the spectral function are denoted 



as X ma x and x. 



average 



is the result of the numerically evaluated integral 



x. 



average 



1 



7T 



ds sd(s) . 



(4.87) 



2ni 



The last column of both tables gives the normalization value of the spectral function. 
All numbers after the fourth digit are dropped. 
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Figure 4.20: Scfxfi-model poles for gs<p<p G [0.1 GeV, 1.0 GeV] in the second Riemann sheet 
with Agscf>cf> = 0.05 GeV, where the gray line marks the branch cut. The dark 
blue dots correspond to the case of fixed Mbw and the red dots to the case 
of fixed Mq . 
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9S(t>4> ^(OJpole ^(l)polfi x (m)pole ^(OJBW ^(1)BW ^(m)BW ^max ^average M Norm 
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Table 4.3: Selection of masses (in units of GcV) for the Scfxfi-model with A = 
1.5 GeV, M BW = 0.6 GcV and m = 0.139 GeV. Here, M BW is fixed. 



94 



9Sct>4> ^(OJpole x (l)po\e ^(mjpole x (0)BW X (1)BW x [m)WW x max ^average Mq 
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Table 4.4: Selection of masses (in units of GeV) for the S^x^-model with A = 
1.5 GeV, M = 0.6 GeV and m = 0.139 GeV. Here, M is fixed. 



4.3.3 Couplings g SH> G [1.0, 2.8] 

No poles are found in the first sheet. The resonance poles descend deeper into the lower 
half plane, displaying a huge difference between the two mass definitions: the dark blue 
pole turns right (as in the case of a sharp cutoff), while the red one is twisted to the left 
(as in the non-relativistic case), yet one can see that this curve becomes flatter (see Fig. 



4.21 ). We also realize that the red poles leave the blue one behind, e.g. the red gap pole 
keeps approaching the first branch point with a much higher parameter velocity. This 
makes sense, because the limit 



ffS^(2m - ie) 




Ml - (2m - r]) 2 + ie 



r] — ie) + Disc S(2m — rj — ie) 

(4.88) 
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is larger than the result in Eq. (4.54) for every point z < 2m on the real axis, namely 
for Mq = Mbw the overall ratio is 



S(M) — Re S(Mbw) 
£(M) 



(4.89) 



as long as S(M) > ReX(MBw)- The latter, though, depends on the choice of the 
Breit-Wigner mass. The numerical bound in order to make the upper condition hold is 
Mbw > 0.4721 GeV for the smooth cutoff (this will turn out to be the same number so 
as to make tachyonic excitations appear for sufficiently large values of the coupling). 
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Figure 4.21: Scfxfi-model poles for gs<p<p G [1-0 GeV, 2.8 GeV] in the second Riemann sheet 
with Z\(?s00 = 0.1 GeV. The dark blue dots correspond to the case of fixed 
Mbw an d the red dots to the case of fixed Mq. 
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Figure 4.22: 5</>0-model gap poles for gs<f><f> G [1.0 GeV, 2.8 GeV] on the real axis of the 
second Riemann sheet with Ags^ = 0.1 GeV. The dark blue dots correspond 
to the case of fixed Mbw and the red dots to the case of fixed M . 
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4.3.4 Couplings g s <p<p G [2.8, 5.0] 

For this range of the coupling we are once again faced with an interesting behaviour. As 
we increase gsss-i f° r some values the gap poles hit the first branch point at z = 2m and 
vanish, possibly they slip into the first Riemann sheet. Although these strange poles 
have no physical interpretation, their existence below threshold lasts as long as gs<t><t> is 
non-zero. For fixed Mq the master solution 



9s<p<t>( M ~ ie ) 



Mi 



M 2 + ie 



S(M - ie) + Disc S(M - ie) 



(4.90) 



yields a positive and real-valued result for the coupling as long as the bare mass is larger 
than 2m, because for M lying below threshold the self-energy T,(M) and the discontinuity 
are both purely real and positive right below the real axis. It is very important to realize 
that the pole must lie below the real axis. Otherwise, the root inside the expression for 
the discontinuity would be evaluated in the wrong quadrant giving an imaginary solution 
for gs<jxj> • Indeed, we find the gap poles leaving the second sheet after they arrived at 
the first branch point. In the case of fixed Mq this happens for 



gs<p<p(2m - ie) 



lim i / — 

V £(2m 



'M 2 - 4m 2 



M 2 



(2m — rj) 2 + ie 



rj — ie) + Disc £(2m — r\ 



S(2m) 



167r 2 (m 2 - A 2 )i(M 2 - 4m 2 ) 
\| A( v / m 2 - A 2 + (m 2 + 2A 2 ) arctan ^' 2 ~ a:2 ' 



(4.91) 



or numerically gscj>6 = 2.935 GeV for our choice of parameters. For fixed Mbw we find 
9s<f><f> = 4.083 GeV. It is evident that there consists a huge difference inside the model 
when we either hold the bare mass constant or the Breit-Wigner mass. As one can see 



from Fig. 4.23 not only the resonance pole lies deeper in the lower half plane, also the 
gap pole on the real axis of the second sheet slips much earlier through the branch cut. 
Both gap poles finally emerge in the first sheet left from the first branch point due to 
the outcome of the previous section and the limit 



gs</><t>(2m - ie) 



lim 

r?->-0+ 



I M 2 



(2m — r]) 2 + ie 



£(2m — rj — ie) 



'Mi 



Am 2 



S(2m 



(4.92) 



This equation has no real solution for gs^ < 2.935 GeV and is purely real for larger 
Note that the new left pole for fixed bare mass tends to the origin in the first 
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Figure 4.23: S^x^-model poles for gs<f>4> G [2.8 GeV, 5.0 GeV] in the second Riemann sheet 
with Agsffxf, = 0.1 GeV. The dark blue dots correspond to the case of fixed 
Mbw and the red dots to the case of fixed M . 



sheet and reaches it for 



lim , 

M->o+ V £(M 



M 2 + 



te 



ie) 



Ml 



S(0 



if 



16M 2 vr 2 (m 2 



A 2 ): 



\ A 3 ( - 3AVm 2 - A 2 + (m 2 + 2A 2 ) arctan 



v/7 



-A 2 



) 



(4.93) 



or numerically gs<p<f> = 4.184 GeV. 

In the end we may illustrate the new left pole showing up in the first sheet simply by 
considering the spectral function for M < 2m: 



d(M) 



2M(e + g 2 s ^ImZ(M + ie)) 



( M 2 _ M 2 + g 2^ Re S(M + ie) y + ( £ + g 2^ ImS(M + i£) y 
2Me 



( M 2 _ M 2 + Re S(M + ie) y + e 2 
2Mvn5(M 2 -M 2 + 5 |^Re£(M)) . 



(4.94) 



From here we also see that the Breit-Wigner mass a?n)BW equals the corresponding 
pole mass (the argument of the delta distribution function is the master solution in 



the first sheet). But, compare Tab. 4.4, once the coupling exceeds a critical value of 
gs4>4> = 2.969 GeV, there is actually no Breit-Wigner mass for the resonance pole in the 
second sheet! This is not true in the case of fixed Mbw> both for sharp and smooth 
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Figure 4.24: Left Scfxp-model poles for gs </></> £ [2.8 GeV, 5.0 GeV] on the real axis of the 
second Riemann sheet with Ags<f,<f, = 0.1 GeV. The dark blue dots correspond 
to the case of fixed Mbw an d the red dots to the case of fixed Mq. 
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Figure 4.25: Scfxfi-model poles for gs<p<p G [2.8 GeV, 5.0 GeV] in the first Riemann sheet 
with Z\(?s00 = 0.1 GeV. The dark blue dots correspond to the case of fixed 
Mbw an d the red dots to the case of fixed Mq. 



cutoff, since we forced all equations to satisfy 



Af| w - M 2 + gl^ Re S(M BW ) = 



(4.95) 



Figure 4.25 shows the left pole arriving at the origin for a positive value of the coupling 
and after that it becomes tachyonic in the regarded range of the coupling for 




E(7? - iO) 



(4.96) 
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4.3.5 Couplings g SH> G [5.0,9.0] 

The only thing that is left is to find the coupling for which the left pole on the real axis 
becomes tachyonic in the case of fixed M B \y. We already know the limit 



9S<M>{0 ~ *e) 



lim 

M-S.0+ 



' M 2 
Re E (M BW 



A /S(0-ie)-ReE(Af B w) 
—iMbw 



ie 

S(M-ie) 

Re£(M BW ) < S(0-ie) 
Re£(M BW ) > S(0-*e) 



(4.97) 



A /ReS(A/ B w)-S(0-ie) 

where the first is true with value gs6<t> = 8.748 GeV. Of course, this is the same as 



lim 



1 (r? + iyf - Mg w 
ReS(M BW )-S(r ? + %) 



Mr, 



^/ S (7?-i0)-ReS(M B w) 

-iM B w 
^/RcS(Af BW )-S(»?-iO) 



Re£(M BW ) < Efa-iO) 
Re E(Mbw) > Efa-iO) 



(4.98) 



Hence, the tachyonic pole appears when the left pole vanishes. This statement is only 



valid for the first case in Eq. (4.98). It depends on the choice of the Breit-Wigner mass 



M B w and although we possess an expression for the real part of S(M) we cannot perform 
further studies in that direction in virtue of the transcendental character of our equations. 
The numerical bound in order to obtain tachyonic excitations is M B w = 0.4721 GeV. 
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Figure 4.26: poles for gs<j>c/> & [5.0 GeV, 9.0 GeV] in the first Riemann sheet 

with Ags^ = 0.1 GeV. The dark blue dots correspond to the case of fixed 
Mbw and the red dots to the case of fixed M . 
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4.3.6 Spectral function and its normalization 



Let us now have a quick look at the spectral functions d(s) for both cases and especially 
on their normalization. As was shown in Fig. 4.19[ the spectral functions are not 
normalized when using the numerical results for the integral 



7T 



ds d(s) 



(4.99) 



simply because the additional simple poles are not seen by numerics. As usual, we need 
to compute renormalization constants: 



A S (M) 



1 



7T 



ds 



d(s) 



M 2 



s 2 + ie 



(l)polc 



M 2 



X 



(l)polc 



+ ~ 

+ ie 7T 



oo 

2m d " 



d(s) 



s 2 + ie 



(4.100) 



The renormalization constant ^(n po i e is calculated by expanding the inverse of the full 
propagator in a Taylor series around the pole at x = £(i) po ie i n first order and since the 
simple pole and the branch cut lie separated on the real axis at different masses, the 
renormalization constant for the single pole can be extracted from its residue. The final 
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Figure 4.27: Numerical verification that ^ J °° ds d(s) — 1 (blue/red dashed for fixed 
A/bw/Mo) by including the single pole of the first Riemann sheet, and nu- 
merical contribution of the continuous part i f 2 ds d(s) (dark blue and dark 
blue dashed). The necessity of the delta distribution function is evident. 



result equals the one for the sharp cutoff: 

^2 



v (l)pole 



1 + 



dReS(x) 



2x( 1)po i es dx 



The normalization condition can thus be written as 



c (l)polc 



>(l)pole + / ds d(s) . 
1 2m 



(4.101) 



(4.102) 



101 



If we perform our calculation including the renormalization constants, the spectral func- 



tions are then normalized over the full range of gs<j>(j>- The plots in Fig. 4.27 show the 



completely fulfilled normalization condition for our choice of parameters (the used points 



are those in Tab. 4.3 and 4.4). For the sake of completeness, we provide a compilation 



of selected spectral function plots in Appendix G 



4.3.7 Decay width(s) 

It is reasonable to have a look at the full decay width of the resonance in both cases, as 
two times the negative imaginary part y( m ) po i e = — r/2 of the complex pole descending 
on the lower half plane of the second Riemann sheet. In the next figure the first two 
panels show plots of the full width (red) in comparison with the tree-level result 



M 2 i 
^ 77T 



\ree{M) = - 2 — G(M - 2m) (4.103) 



AnM 2 ( 1 + ^_ 



for a constant Breit-Wigner mass (blue, dashed) and the average width (blue) from Eq. 



(4.70) as a function of the coupling. For small couplings gs<jxj> < 1 GeV there is not a 
large difference between those three, and they all vanish for gs^ — > 0. The situation 
changes completely for larger values: the full width becomes larger than the tree-level 
result, which again is larger than the average width. While this is true only for the latter 
two over the whole range of gs<f><f>, the tree-level result passes the value of the full width 
a t gs<j><j> = 5.044 GeV and is more than two times the full width when arriving at the 
end of the interval. It therefore cannot be regarded as a good approximation for large 
couplings, though it gives 'the best fit' for all values below the intersection point. On 
the other hand, the average value of the decay width is far too small in the intermediate 
range and high-coupling limit. 

We can also study other decay widths as a combination of spectral function properties 



and the tree-level result. In the lower panels of Fig. 4.28| we provide plots of the full 



(red) and average width (blue), and plots of the tree- level results for each pole mass 
a; (m)poie (green, dashed), each bare mass Mq (magenta, dashed), each maximum value 
^max of the spectral function (brown, dashed) and each average mass ^average (black, 
dashed). All of them may differ in several ways, but overall they are comparable and 
still yield too small values. The black dashed curve surprisinglyhas an intersection point 
with the full width at gscjxfr = 6.195 GeV. 

An analogous discussion is obtained for the full width (red) in comparison with the 
tree-level result for a constant bare mass (magenta, dashed) and the average width (blue) 



in Fig. 4.29. It is not necessary to repeat all similar properties again - the important 
point is the fact that up to intermediate values of the coupling the tree-level result is not 
a better approach than the average width, it does not even produce 'the best fit' , which 
is clearly given by the tree- level result with each average mass (black, dashed) . The plots 
of the tree-level results for each pole mass £( m ) po ie (green, dashed) and each maximum 
value x max of the spectral function (brown, dashed) are also showing comparable, yet 
too small, values for intermediate and large couplings. 
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Figure 4.28: Comparing decay widths for the 5<^0-model with smooth cutoff and varying 
scales (here, Mbw is fixed): the first two panels show plots of the full width 
(red) in comparison with the tree-level result (4.103) for a constant Breit- 
Wigner mass (blue, dashed) and the average width (blue) from Eq. (4.70) 
as a function of the coupling. The lower two panels are plots of the full 
(red) and average width (blue), and plots of the tree-level results for each 
pole mass a;( m ) po i e (green, dashed), each bare mass Mq (magenta, dashed), 
each maximum value x max of the spectral function (brown, dashed) and each 
average mass a; average (black, dashed). 
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Figure 4.29: Comparing decay widths for the S cjxfi-model with smooth cutoff (here, Mq is 
fixed): On the left, full width (red) in comparison with the tree- level result 
for a constant bare mass (magenta, dashed) and the average width (blue) 
as a function of the coupling. On the right, the tree- level result with each 
average mass (black, dashed) and plots of the tree-level results for each pole 
mass X( m ) po i c (green, dashed) and each maximum value cc max of the spectral 
function (brown, dashed). 



4.3.8 Concluding remarks 

The analysis presented has compared the two cases of a fixed Breit-Wigner and bare 
mass, in particular we specified Mbw = Mo = 0.6 GeV. Some relevant points: 

1. We have shown that in both cases a single seed pole for the u-meson entails two 
other emerging poles, once the interaction term in the Lagrangian is turned on. 
One pole appears in the first sheet and can be discarded as non-physical, induced 
by the finite branch cut, while the other pole is dynamically generated in the 
vicinity of the origin of the second sheet. This pole moves to the branch point 
at threshold and slips through the cut onto the first sheet, where it has to be 
taken into account in the normalization condition of the spectral function. It can 
be interpreted as a bound state. The whole model becomes unstable due to the 
occurence of tachyonic excitations for large couplings. 

2. Although the Breit-Wigner parameterization reproduces very well the mass (and 
decay width) of the resonance pole in the small-coupling regime, this is not true 



anymore for intermediate values of gs<j><j>, compare Tab. 4.3 and 4.4 Above this 
regime, the mass is completely detached from the pole mass, while the average 
mass and the maximum value of the spectral function are sometimes in acceptable 
agreement with the pole mass. 

One can try to connect our introductory discussion about the Breit-Wigner param- 
eterization of narrow resonances with the numerically obtained pole positions on the 
second Riemann sheet. If we start with a single seed pole for the u-meson, the trajec- 
tory for fixed bare mass Mq can be taken as the path of the resonance pole for different 
values of the coupling. The spectral function describes this motion due to the fact that 
poles slightly below the branch cut dominate its shape through the maximum and full 
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width at half maximum - these are the appropriate parameters to describe the unstable 
particle. Though we have not focussed on the full width at half maximum, we listed the 



maximum of the spectral function in Tab. 4.4 from which we see very good agreement 
with the pole mass and the Breit-Wigner mass in the case of small and up to interme- 
diate couplings. The whole situation changes when the pole in the second sheet starts 
turning right: both masses become too small and the Breit-Wigner mass has even no 
value for gscj>4> > 2.969 GeV. We plot the latter and the pole mass as a function of gs<t><t> 
in the next figure. The reason for the described difference lies in the real part of the 




ft* I GeV I 



Figure 4.30: Full resonance mass (red) as the real part a;( m ) po i c = M of the complex pole 
descending on the lower half plane of the second Riemann sheet compared to 
the Breit-Wigner mass £( m )BW (blue). 



self-energy. As can be seen in Fig. 4.17, the self-energy has a square root cusp near 
threshold that affects the solution of 



x (m)BW - M l + 9SH> Re S 0(m)B\v) = 



(4.104) 



dramatically. Such influence is expected near every threshold and is well-known, see for 
example Refs. |29||31|. The resonance is too broad to describe it with the help of the 



two equations (1.33) and (1.35), while one could still try to fit the spectral function by 



a Breit-Wigner distribution. 

On the other hand, the Breit-Wigner width is usually defined by Eq. 
include loop contributions, whereas in effective field theories one instead identifies the 



(1.35) if we 



(preliminary) tree-level result (4.14). Yet, since the quantum fluctuations from mesonic 
loops cannot be neglected in the intermediate and large-coupling regime, we may choose 



the modified expression (4.103), so that the Breit-Wigner width just gains an enhance- 
ment factor, 



namely the renormalization constant evaluated at 
significant effect (see Fig. 



'(m)BW- 



This has a 



4.31): the Breit-Wigner width maps the full width much 



better as long as £( m )BW is sufficiently far from threshold. Otherwise, it grows beyond 
all bounds, while the expression without renormalization stays finite but too small. 
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Figure 4.31: Full decay width (red) as two times the negative imaginary part y( m ) po ic = 
— r/2 of the complex pole descending on the lower half plane of the second 
Riemann sheet compared to the Breit-Wigner width (blue) and the tree-level 
result (blue, dashed) evaluated at the Breit-Wigner mass. 



4.3.9 Determination of the cr-meson or / (500) state 

It was stated in the first chapter that many works in recent years attempted to determine 
the o"-pole by using Roy equations with crossing symmetry, analyticity and unitarity. 
Similar results were found, which actually differ quite a lot from the values of Tornqvist 



and Roos 27 in 1996. We have quoted one result by Caprini et al. |28| according to 
whom the 7T7r-scattering amplitude would contain a pole with the quantum numbers of 
/o(500): 

^s~T e = (441^ 6 - i272± 9 12 5 ) MeV . (4.105) 

The Particle Data Group 2012 also just recently revised its values for the mass and width 
estimates on the a-meson |20] , giving a range 

= (400-550) - i(200-350) MeV . (4.106) 

Our very simple model can indeed reproduce this pole in the case of fixed bare mass by 
setting 

M = 0.5447 GeV , g SH> = 2.783 GeV . (4.107) 

The pole lies inside the PDG estimate. Besides that, we notice the left pole below 
threshold to exist in the first sheet with coordinates 

z pole = (0.2705 - it) GeV . (4.108) 

In other words, if we take the cr-pole of Caprini et al. as a correct determination of its 
mass and decay width, then we additionally would get another pole, interpreted as a 
bound state, slightly below threshold. Both are then one and the same manifestation of 
the seed state included in our Lagrangian. Note that we cannot exclude the possibility 
of having found a spurious pole not belonging to a correct description of nature since no 
experiment has yet indicated a hadronic particle with a mass lower than the two-pion 
threshold. 
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Figure 4.32: Trajectory for the cr-pole (white dot) in the second Riemann sheet when the 
result of Caprini et al. is used. The black dot marks the position of the 
dynamically generated pole (bound state) in the first sheet. 
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5 Conclusion and outlook 



In conclusion, we studied the propagator poles of a scalar hadronic resonance within a 
simple quantum field theoretical model in which this scalar state, namely the u-meson, 
decays into two pions. The initial state was settled as a seed pole above the corresponding 
two-pion threshold in the limit of vanishing coupling. We investigated how the motion 
of the propagator pole(s) is influenced by loop contributions of the two pions if the 
interaction is turned on. For this, the self-energy diagram of the two-pion loop was 
calculated analytically and taken into account by applying a Dyson resummation scheme 
- the scalar propagator was then continued into the second sheet of the Riemann energy 
surface in order to numerically look for poles there. This method was also demonstrated 
for a non-relativistic Lee model. 

The Lee model exposed some interesting deviations from the general transition prop- 
erties of unstable quantum systems as induced by the purely exponential decay law. 
Although the Breit-Wigner parameterization was fully adequate in the small-coupling 
regime, it completely failed for intermediate and large values. It seems quite obvious that 
strongly coupled unstable quantum systems lose this feature and should be described by 
including the full propagator and the complete spectral function, respectively. Moreover, 
two poles were dynamically generated in the first sheet within the Lee model, either ab- 
sent in the continuous spectrum. The one above threshold, discarded as non-physical, 
was assumed to be a direct consequence of a non-smooth cutoff function which therefore 
should be avoided in general. We nevertheless interpreted the pole below threshold as 
a bound state emerging from the strong coupling to the continuum of states. The pole 
trajectory in the second sheet revealed that the unstable particle has a constant mass 
for large couplings, thus, while the interaction with intermediate states first lowers its 
mass, it finally stabilizes it. When this complex seed pole is taken to mirror the mass 
and decay width by 

r 

£pole = M — i— , 

those values are completely detached from quantities like the maximum of the spectral 
function in the large-coupling regime. There was no evidence for other resonance poles 
created in the second sheet. 

In the case of the a-meson, we compared the two cases of a fixed Breit-Wigner and a 
fixed bare mass, as well as two different cutoff functions for the former. As a main result 
of the analysis based on a smooth cutoff function, a single seed pole for the cr-meson 
entailed another pole emerging in our model, once the interaction term in the Lagrangian 
was turned on. It was dynamically generated in the vicinity of the origin in the second 
sheet and moved to the branch point at threshold, finally slipping through the cut onto 
the first sheet where it had to be taken into account in the normalization condition of 
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the spectral function. We interpreted this pole as a putative bound state belonging to 
a complete description of one and the same object, i.e., the state /o(500). Note that 
by increasing the coupling the whole model became unstable due to the occurrence of 
tachyonic excitations. 

Thus, besides the expected resonance pole in the lower half plane of the second Rie- 
mann sheet, we have found - for a sufficiently large coupling constant - an additional 
pole on the first sheet below the two-pion threshold. It was clarified under which condi- 
tions such a stable state can emerge and it turned out that it factually can be generated 
simultaneously with the <r-pol^] found by Caprini et al. in Ref. |28|. The appearence 
of such a pole was observed for two different kinds of cutoff functions and seems to 
have its origin in the (resummed) two-pion loop. Further studies should indeed target 
the self-energy contribution in more detail and reconsider other form factors, as well as 
different values for the cutoff parameter A. Note that we cannot exclude the possiblity 
of having found a spurious pole not belonging to a correct description of nature (except 
for the claims in the footnote, no experiment has yet indicated a hadronic particle with 
a mass lower than the two-pion threshold). It is additionally worth to mention that al- 
though there exist phenomenological models in which the state /o(500) is interpreted as 
a tetraquark state and where the theory reduces to our Lagrangian (4.1) when only the 



pion-pion channel is studied 23-25 , more advanced investigations have to go beyond 
our simple model. One should especially include chiral symmetry. 

Despite all that, our work provided numerical bounds for the breakdown of the tree- 
level values for mass and decay width of a resonance below 1 GeV characterized by our 
model, as well as displayed the limits of the Breit-Wigner parameterization in the case 
of small and intermediate couplings. When the seed pole in the second Riemann sheet 
moved deep down into the complex plane, especially the Breit-Wigner mass became too 
small and had even no value for some critical coupling due to the square root cusp of 
the self-energy near threshold. Such influence is indeed expected near every threshold 



and is well-known, see for example Refs. [29 , 31 . Concerning the decay width we have 



found that the Breit-Wigner width maps the full width much better than the tree-level 
one as long as the corresponding mass is far enough from threshold. 

The main topic of our future work could be the study of hadronic resonances within a 
specified effective quantum field theory, e.g. the extended linear sigma model (eLSM) as 
was established in the group of Prof. Dr. Dirk H. Rischke, see for example Refs. [78[|79| . 
One should investigate the influence of loop contributions on, for instance, the mass and 
decay width of scalar resonances since these observables are so far calculated only at tree- 
level, leading to the Breit-Wigner parameterization. This seems at least questionable 
for the scalar mesonic sector. A reasonable starting point would be the lightest scalar 
mesons, i.e., the /o(500) state, the /o(980) and ao(980) states, which are not yet included 
explicitly in the mentioned eLSM. It is important to notice that the precise value of 
the (relatively small) decay width of the isovector state ao(980) is still under debate. 



of a novel scalar boson with mass 38 



x It should be mentioned that according to the recent claim 74 -76 

MeV (termed E(38)), this result was already put under further investigation, leading to a publication 
771. 
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Moreover, it can be interpreted as a glueball or kaonic bound state |80| . Bound states 
actually need to be studied by solving a Bethe-Salpeter equation, so this could be another 
interesting project for the future. 

In such a prespective work we would like to clarify if the scalar resonances not yet 
included in the eLSM can be found as propagator poles due to hadronic loop contribu- 
tions. Prom such an investigation one could possibly also learn more about the general 
dependence of the eLSM - and effective field theories in general - on strongly coupled 
hadronic intermediate states, possibly giving new insight into the low-energy regime, 
scalar resonances and both its theoretical description and physical interpretation. 
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Appendix A 



Fourier transformation: 



/d 3 A; 
_____ e « kx y(i_) ; 

/(k) = J d 3 x e- k x /(x) . 
f d 4 k 

four dimensions: f(x) = / - — -j e~ ik ' x f(k) , 

J (27T) 4 

f(k) = I d 4 x e ihx f(x) , 



Laplace transformation 

(• 7+ioo 



/(*) = 2^/ d-e-FW, 

poo 

F{s) = at e~ st f(t) , sec 

JO 



Here, the integration is done along a vertical line with Res = 7 such that 7 is 
larger than the real part of all singularities of F(s). To make a more general 
statement, 7 is chosen in a way that the contour path of integration is in the 
region of convergence of F(s). 

Representation of the delta distribution function: 

1 e 



6 (a — b) = lim 



e ^o+ 7r (a - b) 2 + e 2 



Sokhotski-Plemelj theorem: 

-b 



lim [ dx -^L = V f dx ^ T i7r/(0) , 



1 = pfi) T m5(x), 



a? =L ie \ x 



iv 



with the Cauchy principal value 

pc r pb—rj pc 

V I dx f(x) = lim / dx f(x) + dx f(x) 

J a V^0 + [Ja Jb+v 



• Polar form of a complex number z = x + iy = pe 1 ^ for <fi £ (— tt, tt]: 
p = \/x 2 + y 2 , 







arg z 



arctan - 

X 



arctan 



arctan ^ 

X 



x>0 

tt x < 0, y > 

tt x < 0, y < 

§ x = 0, y > 

-f x = 0, y <0 

undefined x = 0, y = 



Fourier transform of the non-relativistic Breit-Wigner distribution: 



/(*) 



/ 



du e 



-icot 



(5) 



vr (w _ M) 2 + (F)2 
e -iMt-§ |t| 5 t > Q 
e -i(M-q)t 
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Appendix B 



• Natural units: 

c = h=l 

• Minkowski metric: 

„2 ij.2 j„2 



ds z = dt z - dr z , ( gfiu ) = diag(l, -1, -1, -1) 
Normalization of states and commutator relations: 

Ip> = v^4i°) > 

(p|q) = 2^ p (2^) 3 5( 3 )(p-q) , 

[a p ,at] = (2^) 3 <5( 3 )(p-q) 
Spectral representation of Green's functions (with e — > + ): 



1 r c 

QFT: A(p 2 ) = - ds 
n Jo 



,2 



p 2 — s 2 + ie 



1 f°° , 2sp(s 2 ) 
ds 



7T ,/ p 2 — s 2 + ze 



7T Jo 



p2 _ g 2 _|_ ^ 



^(s 2 ) = -ImA(s 2 + ie) , 
d(s) = -2sImA(s + ze) 

QM: G(E) = — [°° duj P ^ , , 
K J~oo E -uj + ie 

p(u) = — Im.G(uj + ie) 
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Appendix C 



Consider a particle S in its rest frame, decaying into two particles A and B. Form 
kinematics we can find an expression for the three-momentum of the outgoing particles, 
namely psab = \Psab\ : 



[M S \ 


/ 












V o J 


\ 



Psab 





-PSAB 





(5.1) 



M\ + M B + 2p 2 SAB + 2J(M| +p| AB )(M2 +pl AB ) , (5.2) 



PSAB 



2M S 



Ml + {M\ - Ml) 2 - 2{M\ - M 2 B )M. 



(5.3) 



A 



B 



Figure 5.1: Schematic decay process S — > AB. 
And for M A = Mb this simplifies to 

PSAA 



1 



'Mi-AMlMl 
2M S v S as 



Ml 



(5.4) 
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Appendix D 



The tree-level decay width r tree (M), when including a non-local interaction term (4.23) 
on the Lagrangian level, is modified by the regularization function /a (I q| = Ps<t><f>) just 
as a multiplication factor. We start in the same way as was done in the case without 



cutoff (see subsection 4.1.2): 

S - i + hi I d 'x T{:S(x) J d A y <f>(x + y/2)<j>(x - y/2)$(y):} 



1 + ig j d A x 

1 + sw 



(5.5) 



where S is the 5-Matrix, T is the time-ordering operator and the dots mark the normal 
ordering prescription. We need again to evaluate the matrix element 



(P1P 2 I5 (1) I P ) 

(P1P2I W j d 4 x 














^g^JO/) + S ( - } J d 4 y 












3$^Hv) ■ Ip> • 





(5.6) 



Here, the crossed-out terms give no contribution because the creation and annihilation 
operators combine in such a way that their scalar product vanishes. Note that the 



s 



s 



<5-sr 



Figure 5.2: Schematic decay process 5 — > 4>(f> after including a non-local interaction. Every 
triangle corner in the new vertex marks a different point in space-time. 
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superscript at the S- and ^-fields denotes those parts of the field that contain a creation 
or annihilation operator, e.g. 



= <j) {+) {x + y/2) 
= f d " p 1 & t e «P-(*+|Q 



The remaining steps are straightforward: 

• /"j4 / i3 / f i3 / / ,3 / ^2£ p 2£ p 2£ p 
= ig d x d°p 1 / d J p 2 / 



(5.7) 



2E p ,2E p ,2E pl (2^f 
x y d 4 y e ^^(-+f) e ^ ^-f) e -V-^ (y)(0 | &p2&Pi6 t >/& t ) , ap , a t|o) 

/ d 4 * / dVi ! dV 2 / dV ^ gE5B / d 4 y e*t-M)e«-(-l)e-*'-*(») 
J J J J ^2E p ,2E p ,2E p , J 

x j>)(p' - p)<^ 3 )( Pl - p 2 )<^(p 2 - pi) + <^(p' " P)^ (3) ( Pl - Pi)<5 (3) (p 2 - p' 2 )} 

J d 4 x J d 4 y|e ip2 '( x+ i)e ipi '^"2)e" ip ' s + e ipi '( x+ 2) e ip2 '( x "2) e " ip ' x |$(y) 

J d 4 x e <(Pi+P»-P)-« | d 4 y jX^) + e ^-(^) | $(y) . (5.8) 



«5 



If the regularization function /a(q) is now simply the Fourier transform of <3?(y) and 
depends only on the magnitude of the three-momentum q, then the invariant amplitude 
—iM is changed only by the former: 

(final|5|mitial) = 2^/ A (|q|)(2vr)V 4 )(pi + p 2 - p) 

= -iM (2vr) V 4 ) (p - p x - p 2 ) , (5.9) 
=> -iM = 2i 5 / A (|q|) , (5.10) 

by remembering that Pi = p/2 + q and p 2 = p/2 — q. Leaving these two momenta inside 



each argument of the regularization function, the integral from Eq. (4.13) becomes 

WM) = f" du Mu2 fl^) 6{u _ y/WJl - m 2 ) , (5.11) 

(2ir)*Mj i 6(n 2 + m2) /M£_ m 2 



where u = \pi\. From this we obtain the real solution 

,2 

~8ttM 2 



r tree (M) = g TZT flM 6(M - 2m) , (5.12) 



with the modified coupling gs^ = \[2g and the magnitude |p 1 | = |p 2 | = ps^ 
\J M 2 /4 — m 2 of the two outgoing ^-particles three-momenta. 
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Appendix E 



For the sake of completeness, we provide a compilation of selected plots of the continuous 
part of the spectral function for the non-relativistic Lee model. From this one can review 
the behaviour described during chapter 3. 
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Appendix F 



For the sake of completeness, we provide a compilation of selected plots of the continuous 
part of the spectral function for the .S^-model with sharp cutoff. From this one can 
review the behaviour described during chapter 4. 
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Appendix G 



For the sake of completeness, we provide a compilation of selected plots of the continuous 
part of the spectral function for the S^-model with smooth cutoff. From this one can 
review the behaviour described during chapter 4. The dark blue curves belong to the 
case of fixed Mbw> the red curves to the case of fixed Mq. 
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